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UNITED STATES 
DEPARTMENT OF THE INTERICR 
FISH AND WILDLIFE SERVICE 





Near Colleaque: 


The Eastern Eneray and Land Use Team (EELUT) is pleased to provide you 
this report on the evaluation of the sensitivity and potential effects 
of acidic deposition in watersheds characteristic of the Rocky Mountain 
Region. This report is part of the series of technical reports on air 
pollution and acid rain developed at EELUT. Previous reports are listed 
on the inside front cover. 


Areas within Rocky Mountain National Park and Yellowstone National Park 
were selected as representative of geologic types in a large portion of 
that in the total Rocky Mountain Region. In addition to determining the 
sensitivity of characteristic watersheds, the study also evaluates the 
extent of current acidification, the impacts on fish populations, and 
recommendations for assessment of future trends in both changing water 
chemistry and fishery impacts. 


Please feel] free to send suggestions or comments on this report to EELUT. 


Sincerely, 


K Ket Schuster 


R. Kent Schreiber 
Acting Team Leader, EELUT 
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EXECUTIVE SUMMARY 


As a result of the growing concern for the potential effects of acid 
deposition in western mountain watersheds, this project was undertaken to: 


1. Determine the sensitivity of watersheds characteristic of the Rocky 
Mountain Region and the relationship of watershed sensitivity to 
geology and soils. 


2. Evaluate the extent of current acidification and the potential for 
increasing acidification with increasing deposition of nitrate and 
sulfate. 


3. Evaluate the results of the above in terms of impacts on fish popu- 
lations. 


4. Develop recommendations for assessment of future trends in both 
changing water chemistry and impacts on fish populations. 


Areas were selected for study which had minimal human impact and for 
which the maximum amount of data on soils, geology and water chemistry already 
existed. The Rocky Mountain National Park (RMNP) and Yellowstone National 
Park (YNP) areas selected exemplified two different geologic types. The 
geology of these areas is representative of a large portion of that in the 
total Rocky Mountain region. In addition, data on precipitation chemistry 
were available from the National Atmospheric Deposition Program (NADP) moni- 
toring studies in the two parks. In Yellowstone National Park up to 30 years 
of water chemistry data were available, along with information on geology and 
soils. In Rocky Mountain National Park a lake and stream sampling program was 
conducted to collect the water chemistry information. Throughout the project 
the relationship between water quality data and fish responses determined in 
research studies in Scandinavia and eastern North America were relied upon to 
assess potential impacts on fish populaticns. 


RESULTS 


The two parks represent a contrast in geslogic materials; Rocky Mountain 
National Park being primarily underlain by granite and Yellowstone National 
Park by volcanic materials, although the geologic material in one watershed in 
Rocky Mountain National Park (the Upper Colorado River Basin) is of volcanic 
origin. The examination of the geochemistry of Rocky Mountain National Park 
has shown that many areas in RMNP are sensitive to acidic deposition and that 
this sensitivity is primarily determined by bedrock geology. In addition, 
sensitivity varies inversely with elevation. The analyses show that water- 
sheds underlain by granite and biotite gneiss and schist are equally 
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sensitive to acidic deposition. The lakes and streams in these watersheds had 
alkalinities <200 weq/l. Waters at the higher elevations (>3300 m) were very 
sensitive (alkalinity < 100 peq/1). The Upper Colorado River Basin and the 
Upper Fall River Basin contains tertiary intrusive rocks in their drainage, 
resulting in low sensitivity (alkalinity >200 peq/!). Waters in Yellowstone 
National Park had alkalinities generally above 200 yq/!l, with a few < 200 yq/! 
on rhyolite or basaltic flows. 


The analysis of sensitive aquatic systems have been extrapolated to the 
Central Rocky Mountain Region by delineating areas underlain by granite 
biotite gneiss and schist and similar gneisses and schists. Areas underlain 
by these formations are classified as sensitive (alkalinity <200 peq/!), lakes 
and streams located at higher elevations (>3300 m) can be classified as very 
sensitive (alkalinity < 100 peq/!). Areas underlain by tertiary intrusive 
rocks were classified as nonsensitive (alkalinity >200 peq/!). 


The Central Rocky Mountain Region is currently receiving precipitation 
that is somewhat acidic (ave. annual pH = 5.0, NADP 1982). The analysis of 
the data collected in RMNP shows that little, if any, acidification of lakes 
or streams has occurred; however, areas that are subject to periodic 
deposition of pollutants during upslope air movements from population centers 
such as Denver may be experiencing some acidification. 


As stated, high-elevation lakes and streams in the Central Rocky Mountain 
Region are very sensitive to acidic deposition. Much of the region is 
underlain by rock with low buffering capacities that is covered by highly 
permeable soils with low ion-exchange capacities. As in RMNP, high-elevation 
,akes and streams in this poorly buffered region probably will become 
acidified if acidic deposition increases to the level currently experienced in 
the northeastern United States. An evaluation of the potential impact of 
increased acidic deposition on specific lakes and streams would require a 
drainage-by-drainace assessment, as local variability in bedrock, hydrologic 
flow path, and soil development may have an overwhelming influence on 
sensitivity. 


With respect to fish population there is currently no evidence of chronic 
acidification and thus no apparent impact on fisheries. However, the very low 
base cation conrentration observed in the headwater drainages of Rocky 
Mountain National Park suggests extreme sensitivity to acidification. Fish 
populations present in these low calcium waters may be particularly 
susceptible to osmoregulatcry stress from episodic acidification. The few 
remaining native trout (Salmo clarki) located in the interior regions of the 
Rocky Mountains persist only in small isolated headwater drainages. The 
displacement of these rare and endangered genotypes to headwater drainages 
also makes them most susceptible to potential acidification in these sensitive 
habitats. Waters in volcanic areas such as Yellowstone National Park are 
generally of high alkalinity and thus do not represent potentially sensitive 
habitats. 





RECOMMENDATIONS 


It is recommended that long-term watershed studies be established at 
several points in the Rocky Mountain region. Water samples should be analyzed 
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for major inorganic constituents, DOC, aluminum, alkalinity and total acidity. 
In addition, surveys of water chemistry should be performed in other 
mountainous areas, not only in the Rocky Mountain region but also in other 
areas of the western United States. In addition, an effort should be made to 
determine the rate of dry deposition of neutral salts, i.e., calcium sulfate. 


Behavioral responses and immigration tendencies of S$. ciarki populations 
exposed to episodes of acidification should be determined experimental ly. 
Potential for aluminur mobilization in the Rocky Mountain watersheds exposed 
to increased acid Geposition should be studied. Studies of winter-spring 
water chemistry in headwater catchments should also include evaluation of the 
movement of trout populations in response to chemical change. The potential 
sensitivity tc acidification of watersheds currently occupied by endangered or 
threatened populations of S. clarki or of watersheds considered as potential! 
candidate sites fur reintroduction should be determined and given special 
consideration. 
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INTRODUCTION 


Acid precipitation and its subsequent impact on aquatic and terrestrial 
ecosystems is a matter of extensive research and debate in much of the 
industrialized world in the northern hemisphere. Northern Europe, United 
States, Canada, and Japan, have al] experienced an increase in the acidity of 
precipitation and in many of these areas this increase has been associated 
with deterioration in aquatic and terrestrial ecosystems. The extent of this 
deterioration and its relationship to acidic precipitation has been the 
subject of considerabie controversy. However, it appears clear at this time 
that in the Scandanavian countries, and limited areas of northeastern U.S. and 
eastern Canada, lake and stream acidification as a result of acidic precipi- 
tation has been demonstrated beyond reasonable doubt. In addition, in many of 
these areas there have been documented declines in fish populations associated 
with the increase in acidity. Effects on terrestrial systems and their 
» sociation with acidic precipitation are less well established, but there do 
appear to be declines in forest productivity in areas of Europe and north- 
eastern United States that are correlated with atmospheric deposition. These 
observations have spurred initiation of extensive research programs to 
determine the relationship between acidic deposition and declines in the 
productivity of terrestrial and aquatic systems (National Research Council 
1981). 

The emphasis of research to date has been in eastern North America, areas 
of which now receive rainfall with an average pH of 4.1 or below along with 
elevated levels of sulfate and nitrate (NADP 1980). In addition, areas in 
which the most significant impacts of acidic deposition have been demonstrated 
in aquatic systems are underlain by geologic material of low buffering 
capacity. Such areas are found in eastern Canada and northeastern United 
States. The combination of limited buffering capacity of soils and parent 
materials and elevated levels of acid deposition have led to lake and stream 
acidification, with subsequent effects on the biotic community. The effects 








on aquatic systems, particularly the decline in fish populations, has been the 
most significant factor in encouraging governments and industry that we are 
dealing with a serious problem which requires not only increased levels of 
research, but also the development of policy to bring about reductions in 
anthropogenic emissions of sulfur and nitrogen oxides. 

While the emphasis in the eastern U.S. and Canada is understandable 
because of the documented problems, in terms of both levels of acidic 
deposition and impacts, it has been noted that the same conditions of 
sensitivity exist in areas of the western United States including the Rocky 
Mountains, the Sierras, Cascades, and other western ranges. These areas are 
primarily of granitic bedrock, and the limited data available has shown that 
the high elevation lakes and streams have sensitivities comparable to those in 
the eastern United States and Canada. On the other hand, data from the 
National Atmospheric Deposition Program (NADP) and other research has not 
demonstrated that these areas are receiving precipitation with acidity 
comparable to that of eastern North America. Preliminary data from NADP would 
indicate that the pH of rainfall] in the Rocky Mountain West is averaging 4.8 
to 5.0, with a possibility that in some areas at higher elevations the pH may 
be somewhat below these values. Precipitation with pH values below 5 would 
generally be considered to be more acidic than normal. This raises the 
question whether the most sensitive lakes and streams at the higher elevations 
are being acidified. Because of the importance of these lakes as trout 
fisheries and their aesthetic value to millions of visitors, it is important 
to gain a better understanding of the current status of these lakes with 
respect to buffering capacity and to assess the potential for acidification 
and subsequent impact on fish populations. A survey of the literature 
indicates however, that little water chemistry data are available for high 
lakes (>2500 m) in the central Rocky Mountains. 

The most extensive watershed study in this region has been at the Como 
Creek watershed in the northern front range of the Colorado Rockies (Lewis and 
Grant 1979). At least 150 weeks of consecutive stream and precipitation 
chemistry data are available (Lewis and Grant 1979; Lewis 1982), and the first 
report of acidic deposition in the Rockies was made at this site (Lewis and 
Grant 1980). The alkalinity in the waters of Como Creek averaged 192 yeq/! as 
determined by potentiometric titration to pH 4 (Lewis and Grant 1979). Other 
Studies of high mountain lakes include those by Dodson (1981) in the Elk 
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Mountains of westcentral Colorado. This area, located on the west slope of 
the Continental Divide, is underlain by quartzite, siltstone, sandstone, and 
shale (Harte et al. 1983). Dodson (1982) found these lakes to have low 
alkalinity, 8 to 350 peq/]. A recent study by Harte et al. (submitted) in the 
same area reveals alkalinities ranging from 8 to 250 yeqg/! with pH vaiues of 
5.9 to 7.88. Turk and Adams (1983) have carried out a study of the chemistry 
of high-elevation lakes in the Flat Tops Wilderness area of western Colorado. 
The bedrock in this area is predominantly basalt with some granitic outcrops. 
Alkalinities in this region range from 70 to 1400 peq/!. In addition, data 
from a lake survey of the Mount Zirkle Wilderness area by Turk (unpublished) 
shows alkalinities ranging from 12 to 315 peq/]. The bedrock in this area is 
primarily granite. In his study of the South St. Vrain Watershed in the Front 
Range of Colorado, Thurmon (in press) found that alkalinities in the head- 
waters averaged 82 yeq/!]. This area is underlain by silverplume granite and 
biotite, gneiss, and schist. 

These studies of water chemistry of headwater lakes and streams in the 
Central Rocky Mountain Region have only recently been conducted. The lack of 
historical data prevents the evaluation of any long-term trend in acidity. 
Lewis and Grant (1979) reported an increase in stream hydrogen ion concen- 
trations in the Como Creek watershed over a 150 week period. This record, 
however, is too short to be considered evidence of increasing acidification. 
Lewis (1982) reported an average decrease in alkalinity of 180 peq/! for 104 
lower elevation lakes and streams surveyed by Pennak in 1938-1942, and 
resurveyed by Lewis in 1979. Although Lewis attributed this loss to acidic 
deposition, the concurrent increase in total dissolved solids (TDS) suggests 
that hydrologic variability is responsible for this decrease in alkalinity. 
In addition, a 180 peq/! loss of alkalinity is much greater than any 
alkalinity loss observed in strongly acidified lakes in the northeastern 
United States (Hendrey et al. 1980). Thus, no convincing evidence of 
acidification in the Central Rockies has yet been presented, although several! 
areas have been shown to be potentially sensitive (alkalinities less than 200 
yeq/|) to acidic deposition. Also, none of the studies in the Rocky Mountain 
Region to date have associated the current alkalinity and pH of high-elevation 
lakes and streams with potential changes in the future and subsequent effects 
on fish populations. 











In general, the objective of this project was to provide an assessment of 
lake and stream sensitivity of selected areas in the Rocky Mountain region and 
to relate this to potential effects on fish populations. It was also felt 
that by coupling water chemistry data with data on the geology and soils in 
the sampled watersheds, that a relationship could be established between water 
chemistry and geology and soils types and that this relationship could be 
extrapolated to other areas of the Rocky Mountain Region. More specifically 
the objectives were to: 

1. Determine the sensitivity of watersheds characteristic of the Rocky 


Mountain region and the relationship of watershed sensitivity to 
geology and soils. 


2. Evaluate the extent of current acidification and the potential for 
increasing acidification with increasing deposition levels of 
nitrate and sulfate. 


3. Evaluate the results of the above in terms of effects on fish popu- 
lations. 


4. Develop recommendations for assessment of future trends in both 

changing water chemistry and effects on fish populations. 

The plan called for selecting areas for the study which had minimal human 
impact and for which the maximum amount of data on soils, geology, and water 
chemistry might already exist. The plan relied heavily on the established 
relationship between water quality data and fish responses as determined in 
research studies in Scandinavia and eastern North America. After considering 
potential study sites available, two areas of different geologic type were 
selected: Yellowstone National Park and Rocky Mountain National Park. The 
geology of these areas is representative of a large portion of that in the 
total Rocky Mountain Region. Another factor considered was that each of these 
parks participated in the National Atmospheric Deposition Program (NADP) 
precipitation chemistry monitoring network, and therefore had available the 
necessary precipitation chemistry data to be used in the analysis of the 
extent of acidification of the surface waters. The existence of information 
on geology and soils in both parks, while incomplete, proved of considerable 
value. In Yellowstone National Park, water chemistry data was available for 
up to thirty years in some lakes. A single season of lake and stream sampling 
was conducted in Rocky Mountain National Park to determine water chemistry. 











The first sections of the report begins with an evaluation of the 
sensitivity of lakes and streams in Rocky Mountain National Park and 
Yellowstone National Park. Sensitivity (alkalinity) is related to both soi! 
and bedrock influences, as well as elevational gradients. Tte following 
section discusses the relationship of the water chemistry data to fish 
populations based on currently available fisheries data. Conclusions and 
recommendations are presented in the final two sections of the report. 














ROCKY MOUNTAIN NATIONAL PARK 


INTRODUCTION 


A discussion of the studies performed in Rocky Mountain National Park 
(RMNP) is presented in this section of the Rocky Mountain Acidification Study 
(RMAS) report. Rocky Mountain National Park is located in northcentral 
Colorado along the Continental Divide. It has an area of 106,700 ha and 
ranges in elevation from 2,329 m to 4,345 m. Approximately one-third of the 
area and the 107 named mountains are over 3,353 m high. The mountain building 
was accomplished by a regional uplift of Precambrian igneous and metamorphic 
rocks during Late Mesozoic and Early Cenozoic time (Richmond 1974). Volcanic 
activity has more recently altered the western side of the park. Alpine 
glaciation during the Pleistocene formed characteristic U-shaped valleys, 
-teep-sided cliffs, and lateral moraines. Cirque lakes are common. These 
geologic characteristics made this park an ideal site in which to evaluate 
lake and stream sensitivity. 


DESCRIPTION OF STUDY WATERSHEDS 


Four watersheds were selected within the park boundaries in which to 
conduct the studies. The geologic criteria used in selecting watersheds were 
geologic control (bedrock and glacial erosion/deposition) as determined from 
existing information, with attention being given to the representativeness of 
selected watersheds to the Rocky Mountain region in general and RMNP in 
particular. Climatic factors suggested the selection of watersheds on both 
the east and west side of the Continental Divide. Access to these watersheds 
by foot trail was also a consideration. 

The specific watersheds selected for study were Fall River Basin, Glacier 
Gorge Basin, East Inlet Basin and Upper Colorado River Basin (Figure 1). 
Information on the bedrock geologic control of these watersheds is available 
with the bulk of the data compiled on a 1:48000 scale map (unpublished) by wW. 
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Figure 1. Location of study areas, Rocky Mountain National Park, Colorado. 
Subbasins are (A) Roaring River, (B) Ypsilon Creek, (C) Upper Fall River, (D) 
Tyndall Gorge, (E) Loch Vale, (F) Glacier Creek. 


A. Braddock, University of Colorado, Boulder. There are also two pertinent, 
unpublished theses (Abbott 1974; Shroba 1977). 

The geology in RMNP is similar to about 70% of the mountainous areas 
within the Rocky Mountain region. The most extensive geologic materia! 
present in the Colorado Rocky Mountain region is Precamorian granites and 
metamorphic rocks. These rock types are present in three of the four basins 
selected for this study: Fall River, Glacier Gorge, and East Inlet Basins. 
The fourth watershed, the Upper Colorado River Basin has as its major rock 
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types Tertiary volcanics with intrusives of andesitic to rhyolitic composi- 
tion. Glacial till also constitutes a major geologic material in al! four of 
the watersheds, but especially in the Fall River Basin. 

Fall River, Glacier Gorge, and East Inlet Basins have tributaries 
beginning at the Continental Divide at elevations above timberline. The 
tributaries start in material that is predominantly rock talus with some 
tundra vegetation present. Soil development above timberline is generally 
weak thus making geologic control] the most significant factor controlling 
water chemistry. Most of the tributaries have as their source some type of 
snow field, glacier or cirque lake. 

The Upper Colorado River Basin watershed, with its bedrock of inter- 
mediate to rhyolitic composition, intrusive and extrusive igneous rocks, has 
approximately analagous geology to that present in the San Juan mountains of 
southwest Colorado and the northwest section of Yellowstone National Park. 
This area is the only section of RMNP where these geologic conditions exist. 
These rock types contribute more in terms of buffering capacity to the soils 
because of mineralogy and more rapid weathering. 

The third broad type of geologic material, glacial deposits, is present 
in all four basins with Fall River having the largest percentage. Glacial 
deposits range from late moraines of up to 200 feet thick to thin veneers in 
upstream regions. These deposits have a major bearing on the soils encoun- 
tered as well as the genera] water quality. The morainal materia! in general 
is coarse grained with abundant boulders, gravels and sands. This material is 
wel] drained and thus water interacts not only with the soils developed on the 
glacial deposits but also with the deposits themse!ves. 

Another factor, climatic regimes, has profound influences on watershed 
characteristics. It is noted that two watersheds, Fall River and Glacier 
Garge basins, are located on the east side of the continental divide, while 
the other two watersheds, East Inlet and Upper Colorado River Basins are on 
the west side of the divide. This provides a contrast in precipitation levels 
and vegetation types, both major factors in soil development. 

A bedrock and surficial geology map using a U.S. Geological Survey 7.5 
minute quadrant topographic map was prepared for each watershed using 
primarily existing data (Braddock, unpublished; Cole 1977; Abbott 1974). 
Also, a map using the same base has been prep..ed showing the geomorphic 
provinces in each watershed. The map of the geomorphic provinces (talus 
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slopes, scoured bedrock, moraine veneer, etc.), shows the prime source of the 
soils (parent material) that have developed in those provinces along with the 
primary vegetation present (Locke, unpub! ished). 

The four watersheds were subdivided on the basis of tributary drainage as 
follows (Figure 1): 


Glacier Gorge: Tyndall! Gorge, Loch Vale, Glacier Creek 
subbasins 
Fall River: Roaring River, Ypsilon Lake, Upper Fall 


River subbasins 
East Inlet: Not subdivided 
Upper Colorado River: Box Canyon subbasin 
Lithologic Units 
Igneous and metamorphic rocks are present in the four watersheds and 
range in age from Precambrian to Tertiary; the metamorphic rocks are confined 
to the Precambrian. Unconsolidated materials of Quaternary age are abundant 


in all basins. The major bedrock units exposed in the four watersheds are 
(Abbott 1974; Cole 1977; O'Neil! 1981): 





Xqs: Biotite gneiss and schist 

Xam: Amphibolite 

Xgg: Quartz diorite gneiss 

bc: Boulder Creek granodiorite 

Ysp: Silver Plume granite 

PEa, PEab, PEap: Andesite flows, breccia flows, prophyry 
Ngd: Granodiorite of Mt. Richtofen stock 

Ngr: Granite of Mt. Cumulus stock 

Nvr: Ash flow tuff 


The mineralogical compositions of bedrock geologic units are summarized 
in Table l. Using the Goldich (1938) weathering sequence and the moda! 
percentages of the lithologic units, a suggested weathering stability 
classification is shown. The susceptibility to chemical weathering increases 
to the right. The classification is essentially based on the relative moda! 
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Table 1. Major bedrock units: Composition and relative 


susceptibility to chemical weathering. 





Lithologic units 
Minerals (by %) 








Nor Ysp bc Nur? Ngd Xgg Xqs Péa, PEab, PEap* Kam 
LOGE WORTROFERTE SO SS ee Re eRe ee HE Se eH eS Most weatherabiec 
Quartz 35 27 31 High 10 30-65 37 None 7-10 
Microcline 60 38 26 High 32 7-40 ll Minor 
Orthoc lase 3 Minor 
Plagioclase 4 24 25 Minor 40 10-30 22 High 20-40 
Plagioclase type” Olig. And. And. Olig And. Labr. 
Biotite l 7 15 Minor 18 15-35 18 High 7-20 
Hornblende 20-60 
Magnetite 2 5-30 
Cordierite 6 
Accessory 2 2 





pDeterminat ion of percentage by optical means was not possible. Rating based on known rock type. 


Olig. = Oligoclase, And. = Andesine, Labr. = Labradorite. 
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percentages of mafic minerals, hornblende and biotite, and plagioclase, as 
these minerals are more susceptible to chemical weathering (Goldich 1938). 
Some attention has been given to the physical stability of the rock, t.e.- 
schistosity in the biotite gneiss and schist (Ags). The location of the ash 
flow tuff (Nvr) in the classification is questionable because of possible 
influence of the moderately high temperature of formation and the glassy 
matrix on weathering. 

Of the four basins selected for study (two on each side of the 
Continental Divide), three are Gominated by metamorphic rocks and soils 
derived from those rocks. The third is in a region of volcanic rocks typical 
of a significant section of the Rocky Mountain region. The three basins 
representing metamorphic materials differ in the percentage of glacial! 
deposits in the basin vs. the steep slope areas with either thin or no soil. 


SAMPLING SCHEME AND METHODS 


Chemistry of Surface Waters 

Conceptual basis. The general objectives of the RMAS were to determine 
the sensitivity of waters in the Rocky Mountain region to acidification by 
acid precipitation and to determine whether this will have an impact on fish 
populations. An earlier section of this report outlined the strategy for 
selecting the sampling sites and descrited how existing geologic and soils 
information had been used to select systems and subsystems in the Park that 
are representative of the Rocky Mountain region. Lakes and streams in nine 
large watersheds (four basins) were sampled on an elevational gradient. for 
each subsystem, it was expected there would be an elevational trend in the 
concentrations of the base cations (calcium, magnesium, sodium, and potassium) 
because these ions are derived from terrestrial ecosystems by the process of 
soil or parent materia! weathering. An elevational gradient in the concen- 
trations of the acid anions (sulfate, chloride, and nitrate) was not expected. 
This concept is illustrated in Figure 2, a theoretical illustration in which 
concentration of the sum of base cations (Cg) and the sum of the acid anions 
(C,) is plotted as a function of elevation. At the point where they 
intersect, the system has lost all alkalinity and is about to develop strong 
acidity. In the lower elevation areas where Cp is larger than C,, there has 
to be another anion for charge-balance considerations. In these systems, the 
ar.ion is bicarbonate, which is equated to alkalinity because of the absence of 
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Figure 2. Conceptual behavior of Ce, Cy and alkalinity concentrations as a 
function of altitude. 


Al and organic compounds that could contribute to the alkalinity. For the 
higher elevations where C, are larger than Ce. there have to be other cations, 
again for reasons of charge balance. The cations are hydrogen and aluminum. 
By definition this region has already developed strong acidity. The effect of 
increased atmospheric deposition of the acid aniuns 04°, No; and C1? will 
cause the horizontal line of Cy to move upward resulting in a larger number of 
lakes at the lower elevations developing strong acidity. This suggests that 
the most sensitive (lowest alkalinity) lakes will be at the watershed head- 
waters. 

The dotted line in Figure 2 is worthy of mention. At lower elevations 
where soil development is more extensive, the process of sulfate adsorption 
may remove atmospherically derived sulfur from the waters of the watershed. 
One of the goals of this project was to determine the degree of sulfate 
adsorption as a function of elevation. Therefore, if the elevational gradient 
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for the acid anions, specifically sulfate, foilows the dotted line, we wil! be 
able to field test this by the determination of sulfate adsorption capacity. 
This concept, relating the concentrations of various species to eleva- 
tion, was the basis for devising the general water chemistry sampling plan. 
Sampling plan. Water chemistry sampling within each sub-basin was 





carried out in a one to four day period to reduce variation in hydrological 
conditions. Samples were collected under base flow conditions, i.e., sampling 
did not occur within 24 hours after rainstorms. Lake samples were collected 
at each lake inlet, outlet, and center location. Stream samples were taken 25 
meters below each confluence and at approximately 150 meter elevation 
intervals. (Sample location maps are presented in Appendix A.) Stream 
samples were collected in the middle of the stream under falling water, while 
lake samples were taken 0.3 meters below the surface at the center of the 
lake. 

Two 250 ml samples were collected at each site in clean polyethylene 
bottles. Each bottle had been washed with hot water and detergent, rinsed 
five times with hot tap water and five times with deionized water. The 
analysis of blank samples indicated no contamination from this cleaning 
procedure. Each bottle was rinrsea three times in situ with sample water 
before filling. One sample was immediately preserved with reagent grade 
chloroform for later base cation and acid anion analysis. The sample 
identification number, location, date, and time of collection were printed on 
the bottle and recorded in the field log book. 

Conductivity and pH were measured in a separate aliquot of sample at the 
site. Unpreserved 250 ml samples for alkalinity analysis were stored at 5°C 
for one to two weeks before analysis. The preserved samples were analyzed for 
chloride, nitrate + nitrite, sulfate, phosphate, silicate, calcium, magnesium, 
sodium, potassium, and ammonium within two to three months afte’ sample 
coliection. The pH of the stream waters was too high to include aluminum 
analysis. 

Field analyses. Conductivities were measured in the field using a 
Beckman mode] RC-16C meter and a Yellow Spring conductivity cell (YSC model 
3404). Corrections to 25°C were calculated according to Standard Methods 
(A.P.H.A. 1976). Field pH measurements were taken with a Corning digital pH 





meter and a Corning mode! 476182 pH electrode. The meter was periodically 
calibrated with cold dilute, strong acids in the laboratory. Each measurement 
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was preceded by a two point calibration using pH 7.00 and pH 4.00 buffer 
solutions. Corrections for the temperature difference between the sample and 
buffer were made using meter adjustments. Stream samples for the field pH 
analysis were collected in polyethylene beakers from the center of the stream. 
The pH was allowed to stabilize before a reading was taken in a quiecent 
sample. Lake samples for the field analysis were collected in polyethylene 
bottles and brought immediately to shore where pH and conductivity measure- 
ments were taken. pH was also measured in the laboratory as part of the 
alkalinity titration. 

Laboratory analyses. Calcium, magnesium, sodium, and potassium were 
analyzed by atomic adsorption spectroscopy using an instrumentation |abora- 
tories mode! 751 AA/AE spectrophotometer. Samples were spiked with a solution 
of lithium and lanthanum to suppress ionization in the magnesium and calcium 
analysis (E.M.S.L. 1978). A standard curve consisting of one blank and five 
Standards was prepared between every 34 samples to check for drift. If 
significant (>10%) drift had occurred, the samples were reanalyzed. The 
instrument was recalibrated after every standard curve (Emme! 1977). 

Sulfate, phosphate, nitrate, chloride, silicate and ammonia were analyzed 
by automated wet chemistry using a Technicon II Auto-Analyzer. Sulfate was 
measured using a modification of the Thorin technique developed by the 
Norwegian Institute for Air Research. Phosphate and silicate were analyzed by 
Standard molybdenum techniques (7.1.S. 1973b; T.1.S. 1976b). Nitrate plus 
nitrite was measured by the standard cadmium reduction technique (1.1.S. 
1972), while ammonia was analyzed using an Indopheno] Blue method (T.1.S. 
1973a). Chloride was analysed using the standard ferricyanide method modified 
for low levels (7.1.5. 1976a). In all of the wet chemical analyses, one blank 
and six standards were run between every 33 samples. 

Alkalinities were measured using a potentiometric method developed by 
Gran (1952). A 50 ml sample was titrated with 0.010 N HC] from pH 4.0 to pH 
3.3. The Gran's function of this titration curve was extrapolated to an 
endpoint to determine the sample alkalinity. A Fisher Accumet 420 pH meter 
and a Corning mode] 476182 pH electrode were used for this titration. The 
meter was calibrated with pH 4.0 and pH 7.0 buffers at the beginning and end 
of each set of titrations. Both samples and buffers used were at room 
temperature and air equilibrium. An initial air-equilibrium pH was recorded 





before each titration. 
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A summary of the analytical methods is presented in Table 2. 


Quality control. 











The study followed the following protocol on quality 








contro}. 

l. Precision: To determine precision, 5% of all water-quality samples 
collected were treated as replicates. The results demonstrated a 
precision for ali measurements of +10%. 

2. Accuracy: In the field, the accuracy of the field pH measurement 
was assured by calibration of the electrode with standard buffers 
before each measurement. The pH meter and electrodes were periodi- 
cally calibrated with cold, dilute acid standards. 

Table 2. Analytical techniques. 

Neasurenent Instrumentation Technique Summary 

parameter 

Mg Instrumentation Laboratory Atomic Lathanum added, aspirated in oxidizing 
Absorption Spectrophotometer, Mode! 751 flame and read at 285.2 nm using 

deuterium background correction 

Ca Instrumentation Laboratory Atomic Lathanum added, aspirated in oxidizing 
Absorption Spectrophotometer, Mode! 751 flame and read at 422.7 nm 

Na Instrumentation Laboratory Atomic Lithium added, aspirated in reducing 
Absorption Spectrophotometer, Model 751 flame and read at 589.0 nm 

K Instrumentation Laboratory Atomic Lithium added, aspirated in reducing 
Absorption Spectrophotometer, Mode! 751 flame and read at 766.5 nm 

50, Scientific Instruments Mode! A200 with Modification of the thorin technique 
custom-designed manifold developed by the Norwegian Inst. for 

Air Research (NILU) 

NO, Scientific Instruments Mode! A200 with Standard Cadmium Reduction technique 
stock manifold 

NH, Scientific Instruments Mode! A200 with Standard indopheno! blue technique 
stock manifold 

C) Scientific Instruments Mode! A200 with Standard ferricyanide method modified 
modified manifold for low levels 

$10, Scientific Instruments Mode! A200 with Standard molybdenum blue technique 
Auto-Analyzer with stock manifold 

PO, Scientific Instruments Mode! A200 with Standard molybdenum blue technique 
Auto-Analyzer with stock manifold 

Conductivity YSI mode! 3403 Cell, Beckman RC-16C meter Standard procedure 

pH Corning Mode! 3 Meter with Corning Mode! 476182 Standard two-point calibration with 
combination electrode pH 7.00 and pH 4.00 buffers 

Alkalinity Radiometer auto burette with Fisher Accumet Air-equilibrated Grans titration 


420 meter and Corning Mode! 476182 
combination electrode 
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The water-chemistry laboratory at the University of Virginia used two 
techniques to determine the accurary of analytical measurement. Standard 
additions were performed on representative samples. In addition, interlabora- 
tory calibration with EPA (1982) and USGS (1981) was conducted by the analysis 
of blind replicates for all ions. The results of these procedures showed that 
the accuracy of the analytical measurements was +10%. 

Soils and Surficial Materials 





William Locke (unpublished) has produced a general map of major soi! 
(landform-vegetation) groups in Rocky Mountain National Park which was 
available for our use. Within each basin six to sixteen sites, representative 
of the major soil groups and providing an altitudinal cross section, were 
selected for soil sampling (see Appendix A map). With the aid of Locke, soi! 
pits were prepared and a soil-sampling scheme devised. 

Soil samples were collected by excavating soil pits by hand to a maximum 
depth of 50 cm. Samples weighing 300 g to 700 g of each soil horizon exposed 
were collected and described. These samples were assumed to be representative 
of the soils encountered in the different landform types in the study area. 
Representative specimens of rock types were collected and identified. 
Rock-type names conformed to those used in the Park by Abbott (1974), O'Neil] 
(1981) and Braddock (unpublished). 

A complete sample description included colors, texture, structure, 
rooting depth, and estimates of coarse material. Slope, size of area, and 
Surrounding vegetation was also recorded along with photographs of the sites. 
Approximately 80-90 samples were collected during the 1981 field season. 

Air dried soil samples passed through a 2 mm sieve were used in all 
laboratory analyses. (The procedures used are described in detail in 
Appendix B). Briefly, they were as follows: pH was determined in a 1:1 soil 
water suspension. Exchangeable bases were extracted with 1 N NH, Ac and 
determined by Atomic Absorption Spectroscopy. Extractable acidity was 
determined using a BaCl,-Triethanolamine solution at pH 8.0, followed by a 
titration. Cation exchange capacity was determined by the sum of exchangeable 
bases and extractable acidity. Organic matter was estimated by loss on 
ignition at 500°C and by wet oxidation in KCr.0,-H,S0, solution (Mebius 
method). Percent sand, silt, and clay was determined by the pipette method 
after dispersion in sodium metaphosphate solution. 
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Selected samples were tested for sulfate absorption by equilibrating 10 
grams of soil with 25 ml] of dilute K,S0, solution at a pH of approximately 
4.1. The solutions were initially at 1, 5, and 10 mg of soil per liter. 
After shaking with soil and separation by filtration, solutions were analyzed 
for sulfate loss. 

Parent and surficial materials. Determination of the parent materials of 
the soils was accomplished by geologic reconnaissance of the area adjacent to 
each soil sample and by binocular microscop< examination of the 2- to 4-mm 
fraction of each soil sample. Glacial till is the parent material for most of 
the soil] samples analyzed and, in turn, is derived from the geologic units 
upvalley. Therefore, the composition of the 2- to 4-mm fraction of the soil 
samples is representative of both the glacial till and its parent material, 

_ the bedrock units in the area. The areal extent of mapped geologic units, 
surficial and bedrock, in each basin was determined using a planimetric 
digitizer and geologic maps. Weathering characteristics of the various 
bedrock and surficial materials were studied, since they directly affect the 
type and amount of soil present. 

Landform types (Table 3) are helpful in determining the parent material 
of the soils developed, the thickness of overburden and the existing vegeta- 
tion. They convey surface features and materials, along with associated 
vegetation that characterize mappable areas. This information can be obtained 








Table 3. Representative landform types (surface features and/or materials) 
observed in research area. 








Exposed bedrock Moraine veneer 
Talus slope Wet meadow 
Unglaciated regolith Dry meadow 
Moraine-ground Wet tundra 
latera! Dry tundra 
end 
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by analysis of aerial photographs. Stream gradients for al] rivers and 
Streams studied were determined with a linear digitizer. 

Seismic refraction surveys were conducted at seven sites to determine the 
thickness and types of surficial materials overlying bedrock. A Nimbus ES-125 
Single-channel signal-enhancement seismograph was used. The surveys were made 
in coniferous forests on moraine veneer, wet meadows and dry meadows. The 
seismic velocities obtained permitted grouping of geologic materials into four 
categories: soil, alluvial deposits, glacial till, and bedrock. Numerous 
factors influenced these values. They include thickness of organic material, 
water content, clay content, percent of boulders, lithologic makeup of glacial 
till and type and degree of weathering and/or jointing of bedrock. 


RESULTS 


Soils and Surficial Materials 





Characteristics. The seismic analysis revealed that generally, all soils 





had velocities corresponding to moist, loamy or silty soils, as reported by 
Redpath (1973), which indicate a low clay content. The glacial ti!! 
velocities correspond to values obtained by Redpath (1973) for dry glacial 
moraine deposits in the Sierra Nevada of California. These deposits had low 
clay content similar to those encountered in the RMNP project area. Most 
bedrock velocities were low, indicating a fairly large degree of weathering 
and/or jointing. This suggests that surface water penetrates to considerable 
depth. 

Seismic profiling revealed the thickness of the soils to be in the range 
of 0 to 1.9 m; alluvial materials, 3.0 to 3.7 m; and the glacial till, 0 to 
7m (Table 4). The high values represent materials in the center of the 
Stream valleys with decreasing thickness up the sideslopes. Although these 
depths are not large by comparison with nonalpine surficial materials, they 
could be quite adequate for huffering of acid precipitation if sufficient clay 
and organic material existed. Stream gradients are generally indicative of 
the landform type and are easily obtainable from topographic maps. The 
Steepness of gradient generally is indirectly proportional to the thickness of 
both glacial and alluvial overburden in the terrain encountered. 

The soils tend to be coarse, stony, and shallow with very little evidence 
of formation of secondary minerals or eluviation. (A detailed description of 
the soil characteristics of each sampling site may be obtained from the 
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Table 4. 


Seismic refraction data summary. 

















Glacial or 
Soil alluvial? Bedrock 
Profile # Landform type 
Velocity (m/s)- Velocity (m/s)- 
thickness (m) thickness (m) Velocity 
l 396-0. 7 1699-5.5 2417 Moraine veneer 
2 215-1.4 1390-6. 2 3135 Dry meadow 
3 224-0.8 1117-3.8 2391 Wet meadow 
4 374-1.9 None present 3941 Moraine veneer 
5 None present 955-3.2 3758 Wet meadow 
6 326-1.5 1594-5. 3 3208 Moraine veneer 
7 318-1.7 None present 2989 Moraine veneer 





*Unconsol idated parent material underlying soi}. 








° 
author). Slight development of a color B horizon was detectable in a few of 


the lower elevation soils developed in glacial till. The major variation in 
soils within the study area was caused by thickness of deposit, depositiona! 
mechanism (local alluvium, til] or colluvial deposit) and amount of accumu- 
lated organic matter. Differences in mineralogy of the parent materia! was 
not evidenced in soil chemical or physical properties. Elevation, especially 
as it relates to vegetative cover, appears to be important in determining the 
organic matter content. 

Coarse materials (>2 mm dia.) comprised much of the soil volume in most 
locations. Some of the meadows, however, were essentially stone-free. The 
talus slopes and most other locations had soils with 10 to 85% coarse 
fragments in the samples collected. The collection process itself was biased 
against large fragments therefore the percentages are conservative. 

In general the percent base saturation and pH of both the organic layers 
and mineral soil are low. All pH values are acid, in the range 4.2 to 5.6, 
with occasional values as low as 3.5 and as high as 6.8. The basic cation 
status of these soils is generally quite low due to the scarcity of basic 
rocks in the underlying material. The % base saturation of the soil exchange 
complex ranges from as low as 2% upwards to 90%, but very few samples exceed 
70%. The average for each of the basins ranges from 24 to 48%. 

The pH occurring most often in the mineral soil horizons was 5.2 with a 
strong central tendency around 5.2. In 15 of the 30 locations the mineral) 
soil pH in the top 50 cm was within the range 5.0-5.3 (Figure 3). The surface 
organic layers had a pH range similar to the mineral soils, but the distribu- 
tion was shifted slightly toward more acid values. 

The clay content of these soils was low. Only two sample locations 
contained horizons with more than 20% clay (% of <2 mm fraction). It is 
apparent that most of the exchange capacity is provided by the organic matter. 
Clay mineralogy was not determined in these samples, but the residual CEC 
which appears to be due to the clay after subtracting the CEC commonly 
associated with the organic fraction, indicates a material high in CEC such as 
vermiculite. Shroba (1977) reported mica alteration and formation of 10-18A 
mixed layer clays in Alpine soils and more extensive alteration and the 
presence of vermiculite in the spruce-fir region in the Rocky Mountains. 

Sampling of soils was limited to a few locations within each watershed. 
The variability within watersheds was great and masks many differences that 
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Figure 3. Distribution of pH values in the mineral soil and the 
surface organic layer of al! four watersheds. 


may exist between drainage basins (Tables 5 through 8). Fall River basin 
(Table 7) contained the only soi! samples with pH in excess of 6.0. These two 
mineral soil samples were from the lowest elevations sampled in tivis study. 
As a result, the average soi! pH in that basin was higher than that in the 
other three. Otherwise, the average pH in Fall River basin would have been 
5.3, essentially the same as those in Upper Colorado and East Inlet basins. 

The samples collected in Glacier Gorge Basin (Table 8) are the most acid. 
Only one mineral soil] sample exceeded a pH of 5.2. The mean of all 10 samples 
in that basin was 5.0, relatively acid for young soils. 

The pH of surface organic layers in each basin was usually slightly lower 
than the mineral soil] underneath it. This is a common occurrence in soils 
with a surface accumulation of organic material under coniferous forests. 

A cross section of soil samples were selected for determination of 
Sulfate absorption capacity which can contribute greatly to a system's 
resistance to the impact of acidic inputs. Sixteen soil samples were tested 
at 1, 5 and 10 yg S/m! solution for their ability to adsorb sulfate. None of 
the samples tested removed detectable quantities of sulfate from solution at 
any of the three concentrations. In many cases the soils released low amounts 
of sulfate to the solution. This was probably released from the organic 
matter as a result of the sample drying and rewetting and microbial activity 
during storage and processing. (Samples for sulfate retention determinations 
should be stored moist at 0° C until the test is performed.) Based on these 
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Table 5. 


Summary of soil properties at four locations in the East Inlet Valley. 





Location” Elevation 


Mineral soil (0-50 cm) 





Surface organic 











Vegetation 
100m org% =o pHi SSCS CClay%~—sOrg 8 pH type® 

EI 5 2.90 5 5.5 27 47 22 61 5.0 F 

EI 3 3.11 2 5.9 15 69 14 57 6.0 F 

EI 2 3.17 8 4.8 34 9 5 48 4.7 M 

El 1 3.35 6 5.1 16 23 8 -— = G 

x 5.3 37 5.2 

of = coniferous forest, M = wet meadow, G = grass, T = talus slopes with scattered grasses. 

or 


location see designation on maps Appendix A. 
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Table 6. Summary of soil properties at nine locations in the Upper Colorado drainage. 





Minera! soil (0-50 cm) Surface organic 











Location” aphid Negetat ies 

Org % pH CEC % BS % Clay % Org % pH type 
$G-3 3.05 3 4.8 20 12 14 38 4./] F 
0-1 3.14 8 5.6 39 50 — 31 5.4 F 
BX-1 3.17 4 5.1 22 21 16 38 5.3 F 
HG-2 3.20 2 5.1 20 4 13 52 aa F 
$G-2 3.23 2 5.2 16 25 13 31 5.2 M 
MN-2 3. 32 2 4.8 16 22 — =< - F 
HG-1 3. 32 2 5.8 17 45 =< -- -- T 
MN-1 3.35 4 5.5 16 33 7 -- == T 
$G-1 3.47 3 5.2 22 44 19 -- - T 
x 5.2 28 4.8 





af = coniferous forest, M = wet meadow, G = grass, T = talus slopes with scattered grasses. 
For location see designation on maps Appendix A. 
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Table 7. Summary of soil properties at sample locations in the Fall River drainage basin. 
‘acatiadl! Elevation Mineral soil (0-50 cm) Surface organic vegetation 
Org % pH CEC BS % Clay % Org % pH type 
(meq/100 g) 
L 81-18 2.62 6 6.1 16 56 6 = -- F 
L 80-6 2.65 <1 6.6 3 85 6 — -— G 
L 80-1 2.68 <1 5.3 2 52 “ 29 5.9 F 
yCS-1 3.08 <1 iF 6 30 4 = = F 
81-13 3.20 l 5.3 ll 25 7 = = F 
81-16 3.23 -— = - _ 60 4.8 M 
80-11 3.57 <1 5.9 8 68 4 -- -- T 
80-13 3. 66 4 5.0 14 23 7 21 4.8 G 
x 5.6 48 5.1 





2 = coniferous forest, M = wet meadow, G = grass, T = talus slopes with scattered grasses. 
For location see designation on maps Appendix A. 
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Table 8. Summary of soil properties at sample locations in Glacier Gorge Basin. 





Mineral soil (0-50 cm) Surface organic 











Location” — wegetation 
Org % pH CEC BS % Clay % Org % pH type 
(meq/100 g) 
GC-4 2.80 l 5.1 8 20 9 65 4.9 F 
1G-4 3.08 5 5.0 26 l 19 25 4.6 F 
GC-3 3.08 3 4.6 12 2 9 62 4.3 F 
TG-3 3.14 l ne 17 39 ] -- -- T 
AC-3 3.20 5 5.1 20 14 6 65 5.5 M 
1G-2 3.35 9 4.3 48 15 10 55 5.2 G 
GC-2 3.41 4 5.2 9 39 3 46 5.1 M 
GC-1 3.44 2 5.2 4 37 -- -- T 
AC-1 3.47 13 4.4 32 16 16 = -- G 
TG-1 3.51 7 5.8 19 52 7 -- -- T 
x 5.0 20 24 4.9 





HF = coniferous forest, M = wet meadow, G = grass, T = talus slopes with scattered grasses. 
For location see designation on maps Appendix A. 


I —_— 








tests, it is obvious that these soils have very little if any sulfate absorp- 
tion capacity at S concentrations up to 10 yg S/ml at a pH of 4.1. 

Soil buffering capacity. The soils of the Rocky Mountain National Park 
at the higher elevation (above 2800 m) are only slightly developed. They have 





surface organic horizons in some cases and usually have a darkened A horizon 
underlain by slightly weathered material derived from granites, diorites and 
other rocks low in bases. The soils are coarse, low in clay, low in basic 
cations and relatively acid (modal pH 5.2). The organic matter provides most 
of the cation exchange capacity and there is no evidence of sulfate adsorption 
Capacity. 

The physical properties of these soils--coarse, high in sand, low in clay 
and steeply sloping--promote rapid movement of water to the streams and lakes. 
Residence times of water will tend to be short especially at the higher 
altitudes where there is less vegetation, thinner soils and less accumulated 
organic matter. Much of the water flow will be rapid and through coarse 
channels such that only a small percentage wil! approach equilibrium with the 
Soi] materials. 

The chemical properties, low base saturation, low pH, and low cation 
exchange capacity provide little opportunity for neutralizing acidic inputs or 
providing significant alkalinity to the water system. If acidic deposition in 
the region reached pH values in the low 4.0s there would be a potential! for 
the soil to buffer the pH upward slightly. However, since most water moving 
through these coarse soils will not have a residence time sufficient for 
equilibration with the soil, there would be little effect on acidity. The 
current situation with the soils more acid than the streams indicates that the 
soils are not controlling the aquatic system pH. Since the dominant soi] pH 
is only slightly above 5 and the exchangeable basic cation supply is low, 
these soils would not be expected to provide strong buffering of the ecosystem 
against acidification. 

Even though aluminum was not determined in this study the present pH of 
these soils indicates that aluminum concentrations in the soil solution are 
not high. A lowering of soil pH to 4.5-4.6 would cause a signifi.« increase 
in aluminum in soil solution. Some additional aluminum would reach the 
aquatic system, but much more information would be needed to quantify the 
expected change in aluminum. 


26 


"Ayr? Perey A 


BEST COPY AVAILABLE 








The soil properties already mentioned plus the absence of sulfate absorp- 
tion capacity indicate that these soils should be considered in a “sensitive” 
category. They are among the soils which have low capacity to absorb and 
neutralize acidity and therefore the associated aquatic systems would not be 
protected from pH change by the soil ecosystem (Galloway et al. 1983). The 
soil itself is probably resistant to rapid acidification. The abundance of 
relatively unweathered minerals provides the soil's major buffering against 
lowering of the pH much below its present values. 

Chemistry of Surface Waters 

The lakes and streams in RMNP are characterized by low ionic strength. 
An average concentration for the inorganic constituents are presented in 
Table 9. The raw data for all of the measured constituents of the more than 
150 samples are presented in Appendix C, along with sampling dates for each 
subbasin. On an equivalent basis, the relative concentrations of the base 
cations are Ca>>Na>Mg>>K for most of the watersheds, although magnesium 
concentrations are greater than sodium in the Upper Colorado and Upper Fal] 





River Basins. 

The Ypsilon Creek, Tyndall Gorge, Loch Vale, and Glacier Gorge subbasins 
had similar stream chemistry. These watersheds are characterized as having 
lower alkalinity and base cation concentrations than the rest of the study 
area. Roaring River and East Inlet had siightly higher alkalinity, base 
cation, and silica concentrations. The Upper Fall River Basin had signifi- 
cantly higher concentrations of these constituents, while the Upper Colorado 
River Basin had the highest alkalinity, base cation, and sulfate concentra- 
tions of any of the watersheds. 


DISCUSSION 


Relationship of Surface Water Chemistry to Atmospheric Deposition, 
Elevation, Soils, Surficial Materials and Geology 
Calculation of “excess” cations. The streams and |lakes of RMNP have very 











low cation and anion concentrations. In the watersheds underlain by granitic 
and metamorphic rock, the average Cp (Na + K * Mg + Ca), is only ~85 peq/! 
(Table 9). In such dilute waters, the atmospherically deposited base cations 
may comprise a significant fraction of the stream water cation concentration. 
Previous researchers have corrected surface water chemistry data for the 
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Table 9. Mean ionic concentrations in Rocky Mountain National Park watersheds. ° 

Watershed pH Alk Na K Mg Ca NH, Cl NO, SO, $10, Cy Cp 
Roaring River 6.9 744.4 29.0 3.6 21.7 66.5 0.1 6.0 8.3 34.7 68.8 123.4 121.0 
Ypsilon Creek 6.6 48.0 229 41 16.2 43.6 0.1 3.4 9.8 30.2 61.3 91.4 86.2 
Upper Fall River 7.1 1805 40.2 7.5 67.3 106.9 0.1 7.9 4.9 46.3 110.8 239.6 222.0 
Andrews Creek 6.5 38.8 16.1 3.7 13.1 55.7 0.0 3.7 12.5 32.3 %.5 87.3 88.9 
Loch Vale 6.0 41.2 16.0 3.1 13.9 52.9 0.4 4.5 17.1 26.2 33.3 91.0 8.2 
Glacier Creek 6.6 40.3 14.2 3.0 10.0 4.3 0.1 3.1 11.3 «613.3 32.8 68.0 73.8 
Upper Colorado 

River 7.5 3318 46 7.6 80.0 234.3 0.3 7.9 6.3 84.3 81.5 430.3 357.1 
East Inlet 6.8 8.5 2.5 2.5 16.0 9.0 0.1 4.2 5.6 35.7. 67.9 131.0 125.5 





®a1) concentrations are in pea/l, except $10,, which is in pM/1. 





atmospherically deposited base cations by substracting the “seasalt" contribu- 
tion (Henriksen 1980; Almer et al. 1978). Although RMNP does not receive sea 
salt in its precipitation, the evaporite basins and arid regions upwind of the 
park result in the wet and dry deposition of base cations in the form of dust 
and salts (Gosz 1975). In an effort to separate “excess” cations (those 
released as a result of acid buffering reactions) from atmospherically 
deposited cations, we have made a correction for atmospheric deposition. 
Chloride, a conservative element with no terrestrial source (ij.e., not present 
in the bedrock) is used to make this correction. The calculation is made by 
subtracting from each cation concentration the product of the chloride 
conc7atration at that sample site and the ratio of that element to chloride in 
precipitation (NADP 1982). The result is the precipitation corrected or 
"excess" concentration--e.g., Na* = Na - (Cl x Na recip. / “precip. (3 
indicates “excess” concentration). Some base cations, especially potassium, 
are taken up by the ecosystem. As a result, our precipitation-corrected 
concentration of this element may turn out to be negative. If the calculation 
results in a negative value for an element, the precipitation-corrected 
concentration is taken as zero. The results of these calculations are 
presented in Table 10. This method of calculating the “excess” base-cation 
concentration (Na* * K* + Mg* + Ca* = Cy”) is a correction for the maximum 
contribution of C., (Na + K + Mg + Ca = Cy) by wet deposition. 

The stochiometry of primary dissolution reactions provides a check on 
this correction, since the release of base cations by mineral dissolution is 
accompanied by the concurrent release of HgSi0,, e.g.: 


Nay 7fap gAl, 4515 70, + 3.45 H,0 + 1.3 HCO, = 0.65 Al,S140,¢(OH) 4 


(oligoclase) (kaolinite) 
+ 1.4 HySi0, + 0.7 Nav? + 0.3 ca” + 1.3 HC0;! 


3 Ko (Mg,Fe,)Al,S1,0,,(0H), + 24 H40 + 18 HCO, = 


(biotite) 


2 a | 


, . +) +? ‘ 
2(Mg,Fe,)A1,51<059(0H)4-8 HA0 + 6K" + 3Mg + 3 Fe ; 


* 8 H,SI0, + 18 HCO 


(vermicullite) 
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Table 10. Comparison of mean of uncorrected (U) and 
“excess” (£) base cations. 








Watershed Ce Na K Mg Ca 
Roaring River (U) 121.0 29.0 3.6 21.7 66.5 
(E) 72.3 21.9 0.0 11.7 39.3 
Ypsilon Creek (U) 86.2 22.9 4 15.2 43.6 
(E) 56.2 18.8 0 9.4 28.2 
Upper Fall River (U) 222.0 40.2 7.5 67.3 106.9 
(E) 155.2 30.7 0.0 53.8 71.2 
Andrews Creek (U) 88.9 16.1 3.7 13.1 55.7 
(E) 57.3 11.9 0.0 7.3 38.1 
Loch Vale (U) 86.2 15.0 3.1 13.9 52.9 
(Et) 47.9 9.6 0.3 6.1 32.4 

Glacier Creek (U) 73.8 14.2 3.0 10.0 46.3 
(E) 47.6 10.5 0.2 4.8 32.5 

Upper Colorado River (U) 357.1 45.6 7 77.8 218.1 
(E) 280.1 35.1 1.4 63.0 181.0 
East Inlet (U) 125.5 26.5 2.5 16.0 90.0 
(E) 91.5 22.5 0.0 8.8 60.6 





San) concentrations in peq/!. 


We can assess the accuracy of our estimate of C,* by comparing the corrected 
Ch" to silica concentrations (as moles of $10,°). Graphs of CA" vs. $10,” for 
Loch Vale and Ypsilon Creek watersheds, with low intercepts show that our 
correction is appropriate for the watersheds with lowest alkalinity (Figure 
4). The variation in the slopes of the regression lines among the watersheds 
reflects a difference in the distribution of minerals. Watersheds with a 
slightly greater percentage of biotite will have larger slopes than those with 
more oligoclase in their drainage according to tne different stochiomatry of 
the weathering reactions. This check on our calculation of C,* is not 
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watersheds in RMNP. 


appropriate in watersheds where nonsilicious minerals and soi! cation-exchange 
reactions play a large part in the geochemistry (ij.e., the Upper Colorado 
River and parts of the Fal) River basins). 

Atmospheric deposition and water chemistry. Researchers have shown that 





the predominent direction of air mass movement over the Front Range of the 
Central Rocky Mountains is from west to east (Barry 1973), with periodic 
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upslope movement from the east (Kelley and Stedman 1980). Although most of 
the precipitation received by this area is dropped by the westerly air masses, 
the upsiope transport of pollutants from sources to the east may have a 
profound influence on the deposition chemistry of the Front Range (Kelley and 
Stedman 1980). The atmospheric deposition from the westerly air masses 
contains little sea salt (Junge and Werby 1958) but does result in the wet and 
dry deposition of airborne dust and salts from the dry, windy areas to the 
west (Gosz 1975). As a result, RMNP precipitation has a relatively high base 
cation concentration [22.0 peq/!] Ca in RMNP, compared with 13.0 peg/! in the 
Adirondacks (NADP 1982, Altwicker and Johannes 1980)). 

The nitrate concentrations in RMNP are relatively high, comprising 25.6% 
of the acid anions in precipitation, compared with 28.7% at Hubbard Brook 
(NADP 1982). These high nitrate concentrations are probably the result of 
upslope transport of NO. from the Denver metropolitan area (Kelley and Stedman 
1980). The bedrock in RMNP contains only trace amounts of chloride (Lovering 
and Goddard 1950), while sulfur-bearing minerals (e.g., chalcophyrite and 
gypsum) are found only in the Upper Fall River and Colorado River basins. As 
a result, atmospheric deposition is expected to be the primary source of 
chioride and nitrate 1 all of the streams in RMNP, as well as the primary 
source of sulfate in the waters of East Inlet, Glacier Gorge, Ypsilon Creek, 
and Roaring River. 

Stream water chloride and nitrate concentrations are fairly uniform over 
the park, although East Iniet, on the west side of the continental divide, has 
lower nitrate concentrations than the other four watersheds (Table 9). Two 
plausible explanations for the lower nitrate in East Inlet are (1) less 
deposition from upslope air masses bearing NO. from the Denver area and/or (2) 
more biological activity. Chloride and sulfate concentrations show little 
variability with elevation, indicating negligible variation in the atmospheric 
contribution of elements with elevation (Figures 5 and 6). Nitrate, the most 
biologically active anion, shows some variation over elevation, with the 
highest nitrate concentrations above the timberline, where biological activity 
is lowest (Figure 7). 

Streamwater chemistry and elevation. Elevation is related to several 





parameters that may effect streamwater chemistry. The thickness of surficial 
deposits (glacial till, alluvium, loess, and soils) and the length of flow 
path increase with decreasing elevation as we descend from headwater lakes 
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Figure 5. Graphs of Cl vs. elevation for subbasins in RMNP. 
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situated on scoured bedrock to lower elevation lakes surrounded by moraine 
veneer with thin soil cover. In addition, watershed slope decreases with 
decreasing elevation, further increasing the contact time between precipita- 
tion and neutralizing materials. These parameters are al] difficult to 
measure. Turk and Adams (1983) have shown that elevation used as a surrogate 
for these and other unmeasured watershed parameters (e.g., residence time, 
length of flow path, effective soi] and bedrock area), successfully predicts 
alkalinity. Our analysis shows that alkalinity, Cp” (excess Ce), and silicate 
are all inversely related to elevation in subbasins with homogeneous 
mineralogy and low alkalinities (Glacier Gorge, Loch Vale, Ypsilon Creek, 
Roaring River, and East Inlet) (Figures 8 through 10). 

This relationship between elevation and streamwater chemistry is obscured 
when the data from more than one subbasin are used in the analysis. This 
probably results from the differences in the distribution of glacial til] over 
elevation between subbasins in the same watershed. The relationship between 
elevation and streamwater chemistry is overwhelmed by the variability in 
mineralogy and subsequent weathering rates in the Upper Fall River and Upper 
Colorado River basins (Figure 8). Although elevation contributes to the 
variability of stream chemical composition in areas with homogeneous bedrock, 
variability in mineralogy and other interwatershed variables can overwhelm 
altitudinal effects. 

Relationships between bedrock geology, surficial materials, soils 





and stream water chemistry. The chemical composition of natural waters is 





primarily a function of the interactions between atmospheric deposition, 
bedrock geology, and surficial deposits. Previous researchers have found that 
an area's ability to neutralize acidic deposition is the result of mineral 
weathering and soil] ion-exchange reactions (Norton 1980; McFee et al. 1977; 
Johnson and Cole 1980). Under pristine conditions, the chemical dissolution 
of minerals by carbonic acid results in a loss of hydrogen ion (acidity) and 
the production of alkalinity (acid-neutralizing capacity), base cations, and 
Silica. The rates of these reactions vary with the mineralogy of the bedrock. 
Carbonaceous minerals, such as calcium carbonate, dissolve very rapidly, 
providing “infinite” buffering capacity, while silicious rocks, such as 
granite, dissolve very slowly, offering little buffering capacity. Surficial 
materials also play an important role in buffering acidic deposition, 
as soils neutralize acidic deposition through cation and anion exchange. 
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Figure 8. Graphs of alkalinity vs. elevation for subbasins in RMNP. 
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Surface water chemistry, particularly alkalinity is an integration of these 
mechanisms. As a result, alkalinity 1s frequently used as an index of sensi- 
tivity. 

The geology in the Sark ranges from Precambrian granite to tertiary 
intrusive and extrusive rocks. The difference in the weathering rates of the 
minerals in these formations is reflected by the streamwater chemistry. The 
highest concentrations of CB", alkalinity, and silicate occur in the Upper 
Colorado River Basin, an area underlain by highly weatherable ash flow tuff 
and andesite. The CB" and alkalinity concentrations in Glacier Creek, a 
watershed underlain by Silver Plume granite, are an order of magnitude lower 
than those in the Upper Colorado River Basin (Tables 9 through 10). 

R-mode factor analysis was used to investigate the relationships among 
the water chemistry variables and the bedrock and soil parameters. This 
method of analysis has been used by several authors (Dawdy and Feth 1967; 
Miller and Drever 1977a; Reeder et a). 1972) to irterpret natural water chem- 
istry in terms of its geochemical origin. In these analyses, factor analysis 
is used to reveal the geochemical processes responsible for the stream chemis- 
try composition. A brief review of factor analysis follows. 

Factor analysis is a statistical data reduction technique that rearranges 
the information contained in the correlation matrix for a set of variables 
into a smaller set of independent factors. These factors are linear combina- 
tions of the original variables. The first step in the analysis calculates 
the principal components, or factors that explain the greatest amount of the 
variances and co-variances in the correlation matrix. In the varimax solution 
used in this analysis, these principal components are rotated orthogonally to 
achieve a simple structure. This rotation produces a set of uncorrelated 
factors, so that the factor loadings (the extent to which each factor is 
associated with a particular variable) tend towards unity or zero. The re- 
Sulting factors are interpreted as source variables accounting for the rela- 
tionships between the origina! variables. 

The data set from RMNP was split into two geologic groups for this 
analysis. One group represents all samples from areas underlain by granite 
and gneiss; the other represents the samples collected from watersheds con- 
taining tertiary volcanics. Each group was analyzed with and without a bed- 
rock geology and soil parameter. The analysis of the stream chemistry alone 
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allowed a larger sample size. Variables that had many zero values (K* and 
NH) were not included in the analysis. All factors with e'tqen values 
greater than 0.51 are reported. 

The results of the factor analysis of the granite and gneiss group show 
that four factors account for nearly 93% of the variance in stream water 
chemistry data (Table 11). The first factor, accounting for almost 50% of the 


Table 11. R-mode varimax factor matrix of chemical data 
for 88 lake and stream samples underlain by granite and 
biotite gneiss and schist in Rocky Mountain National Park. 











Factor 
Variable 
1 2 3 4 

o! 
N3 934 0.243 -0. 185 -0.042 
Mg” 827 0. 346 0.092 -0.127 
Ca" 325 0.919 0.144 0. 068 
cl -0.031 0.082 0.112 0.975 
NO, -0. 150 -0.025 0.945 0.078 
S04 0.259 0.437 0.609 0.479 
Sid, 0.937 0.156 0.029 0. 180 
Blk 0.538 0.656 -0. 409 0.277 
% of variance 
explained by factor 49.5 23.6 12.0 7.7 
Cumulative % of 
variance 439.5 73.1 85.1 92.8 





variance, has high loadings for Na*, Mg", $10), and alkalinity. This factor 


is interpreted as an oligoclase and biotite weathering factor. These two 


minerals weather to kaolinite and vermicullite, respectively, releasing Na 7 


Mg’ 510,, HCO; ! and a smal! amount of ca’ 











Ca Si, 50, + 3.45 H,0 + 1.3 H,CO, = 0.65 Al $140, (OH), 


Nay 78 3A!) 3512 7% 2 2°93 2 
(oligoclase) (kaolinite) 
+ 1.4 HySi0, + 0.7 Nav) + 0.3 Ca’? + 1.3 HCO, 
3 K,(Mg,Fe,)Al,Si,0,(0H), + 24.H,0 + 18 H,C0, = 
(biotite) 
; +] +? +? ; -] 
2(Mg,Fe,)A1,Si<059(0H)4°8 H,0 + 6K") + 3.Mg’ + 3 Fe’ + 8 H,Si0, + 18 HCO,” 


(vermicul lite) 


The high percentage of the variance explained by this oligoclase and 
biotite weathering factor suggests that primary mineral weathering is one of 
the processes dominating stream water chemistry. 

The second factor, accounting for 23.6% of the variance, is less easily 
2, C03! and s0,° ca’ 


major bedrock element, although a smal! percentage is contained in the 


interpreted. This factor loads heavily on Ca’ is not a 
oligoclase. ca*é in precipitation, however, is quite high and is also con- 
tributed to by dry deposition of dust and salts. This factor can be inter- 
preted as a dry deposition factor, caused by the deposition of CaS0,, and 
possibly CaC0,. 

The third factor, accounting for 12% of the variance, is an acid deposi- 
tion factor. The high loadings of no,?, 504°, 
between the acid anions and alkalinity, represents the titration of lake and 
stream HCO," by anthropogenic acids. The fact that NO,” loads higher on this 
factor suggests that HNO, may be a more important component of acid precipita- 
tion than H4S0,. The fo:-th factor, explaining 7.7% of the variance repre- 
sents atmospheric deposition but does not appear to involve any buffering 


with an inverse relationship 


mechanisms. 

The addition of percent unreactive rock (quartz * microcline) and percent 
soil organic matter (% 0.M.) to the factor analysis results in a decrease in 
sample size and the loss of the acid deposition factor (Table 12). The first 


factor still represents oligoclase and biotite weathering, although ca’é loads 
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Table 12. R-mode varimax factor matrix of stream chemistry, % soi! 
organic matter, and &% unreactive rock for 40 lake and stream samples 
underlain by granite and biotite gneiss and schist in Rocky Mountain 
Nationa! Park. 











Factor 

Variable 

l 2 3 3 5 
Na” 0.958 —-0. 064 0.111 —--0. 066 0. 168 
Mg” 0.917 0. 208 0.084 -0.090 -0. 154 
Ca” 0.558 0.727 0.354 -0.042 -0. 048 
ci 0.009 0.273 0.212 0.923 -0. 082 
NO; 0. 226 0.851 -0. 212 0. 286 -0.093 
S0y 0.275 0.701 ~0.172 0.543 -0. 116 
Si0, 0.917 0.053 0.024 0.238 168 
Alk 0.741 -0. 066 0.591 0.140 0.049 
% 0.M. 0.101 -0.107 0.064 -0. 087 978 
% rock -0.118 0.085 -0. 953 -0.112 -0. 062 
% of variance 
explained by factor 41.7 25.1 12.4 9.0 6.1 
Cumulative % of 
variance 4].7 66.8 79.2 88.2 94.3 





somewhat higher than in the previous analysis. The second factor still repre- 
sents atmospheric deposition of CaSO, in dust and salts, although NO, also 
loads heavily on this axis. The third factor represents the inverse relation- 
ship between Heo," and unreactive rock, further evidence that primary mineral 
weathering plays a dominant role in stream chemistry. The fourth factor 
represents atmospheric deposition. The fifth factor represents %0.M., which 
is unrelated to any other variable, suggesting that soil ion-exchange 
processes play a very minor role in stream chemistry. The factor analysis of 
Stream samples from areas underlain by tertiary volcanics siows that four 
factors account for 86.4% of the variance (Table 13). 

The first factor, accounting for over 40% of the variance, indicates that 


the weathering processes in these watersheds are very different from those in 
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Table 13. R-mode varimax factor matrix of stream chemistry 
for 33 lake and stream samples underlain by tertiary 
intrusive bedrock in Rocky Mountain Nationa! Park. 











Factor 
Variable 
l 2 3 4 
1 
Na” 0.122 0.911 -0. 225 -0.076 
Mg” 0.835 0. 152 -0. 264 129 
Ca" 0.772 -0.410 0. 138 115 
C) 0.277 0.033 -0.025 955 
al 
NO, -0.03%6 -0.219 0. 946 -0.025 
a 0.853 0.121 0. 185 095 
SiO, 0. 006 0.884 -0. 066 121 
Alk 0.881 0.127 -0. 184 247 
% of variance 
explained by factor 40.3 26.5 10.0 9.6 
Cumulative % of 
variance 40.3 66.8 76.8 86.4 





the granite and gneiss areas. HC0,°, S0,°, ca’é 
3 4 


, and Mg’? all load heavily on 
this factor, representing the weathering of mafic minerals low in silica, as 
well as the dissolution of CaSO, in the Pierre shale of the Upper Colorado 
Basin. The second factor explains 26.5% of the variance and loads heavily on 
na’? and $10, This factor can be described as oligoclase weathering. The 
third and fourth factors, accounting for 19.6% of the variance, each explain a 
single variable and cannot be interpreted in terms of buffering mechanisms. 
The addition of soil and bedrock parameters in this analysis reduced the 
sample size to 16 observations, too few to be used with factor analysis. 

The results of the analyses of the two geologic groups suggest that 
primary mineral weathering is the major mechanism underlying the stream 
chemistry. The weathering of oligoclase and biotite, the dry deposition of 
CaSO, dust, and the titration of bicarbonate by anthropogenic acids appear to 
be the major geochemical processes in the watersheds underlain by granite and 
biotite gneiss and schist. An assemblage of low-silica, high-sulfate minerals 
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appears to be responsible for much of the stream chemistry in areas underlain 
by tertiary intrusives. Oligoclase weathering plays an importart role in the 
chemistry of these watersheds, while no evidence of acid titration of bicarbo- 
nate was found. 

Discriminant function analysis (DFA) was used to examine the differences 
in stream chemistry between the areas underlain by tertiary volcanic bedrock 
and those surrounded by granite and biotite gneiss and schist. DFA dis- 
tinguishes between two populations on the basis of observations of multiple 
variables. DFA may also be used to classify data on the basis of observed 
variables. In this two-group discriminant analysis, one discriminant function 
is formulated by the analysis. This function is a linear combination of the 
measured variables, i.e., 


L = BX, * BK, * . _ Bx : 

where L is the discriminant function, X, are the measured variables, and BW 
are the discriminant function coefficients. The discriminant function L is 
formulated to achieve the maximum discrimination between the two groups; ij.e., 
the variance in L within each group is much less than the variance in L 
between the two groups. The discriminant function analysis also indicates 
which variables differ most between the two groups. The correlation between 
the measured variables and the discriminant function (ranges between -1.0 and 
*1.0) gives an indication of the importance of that variable in differen- 
tiating between the two groups. 

In this application, DFA is used to differentiate between the two major 
geologic groups (tertiary volcanics vs. granite and biotite gneiss and schist) 
on the basis of the water chemistry parameters. The results of the DFA show 
that the two groups are significantly different (Table 14). The correlations 
between the water chemistry parameters and the discriminant function show that 
alkalinity, Mg’, Ca’, and so,? 
difference between the two geologic groups. These are the same variables that 


are the parameters responsible for the 


make up the first factor in the factor analysis of the tertiary volcanic 
group. This result is consistent with the belief that the differences in 
stream chemistry between the two geologic groups are caused by mineral 


weathering processes The reclassification of the data set on the basis of 
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geologic type. Areas with tertiary intrusive rocks present in their catch- 
ments can be expected to have significantly higher alkalinity, magnesium, 
calcium, and sulfate than areas underlain by granite and gneiss. 

Summary. The chemical! composition of the waters in RMNP is a product of 
the interactions between atmospheric deposition and bedrock geology and 
surficial materials. Soils play a minor role. The previous discussions show 
that bedrock mineralogy, atmospheric deposition, and elevation aii signifi- 
cantly affect the stream chemistry in RMNP. 

Primary mineral weathering appears to be the dominant mechanism deter- 
mining the concentrations of base cations, silica, and alkaiinity throughout 
the park. Factor analysis shows that the primary mineral weathering of 
oligoclase and biotite account for almost 50% of the variance in stream 
chemistry in areas underlain by granite and biotite gneiss and schist. In 
areas that contain tertiary volcanic bedrock, the mineral weathering of mafic 
minerals, sulfur-bearing minerals, and oligoclase account for more than 50% of 
the variance in stream chemistry. Soils appear to have little effect on 
Stream chemistry. The soils are highly permeable, low in clays, and very 
thin. Soil organic matter accounts for most of the soil CEC, but % 0.M. does 
not forr a factor with any of the stream chemistry variables. 

Alkalinity, Ce. and sulfate are inversely related to elevation in the 
Glacier Creek, Loch Vale, Ypsilon Creek, Roaring River, and East Inlet sub- 
basins. This is a result of the deeper glacial till (larger flow path) and 
gentile gradients (larger residence time) at lower elevations. 

Atmospheric deposition is the primary source of chloride and nitrate in 
the streams of RMNP. Atmospheric deposition is also the primary source of 
sulfate in RMNP waters, with the exception of the Upper Colorado and Upper 
Fall River basins, where the weathering of sulfur-bearing minerals is a source 
of sulfate. Chloride and sulfate are relatively constant with elevation, 
while nitrate concentrations are highest above the timberline, where biolog- 
ical activity is lowest. Atmospheric deposition is also a significant source 
of stream water Ce. a result of the deposition of airborne dust and salts from 
the dry, windy regions upwind of RMNP. 

Discriminant function analysis of the water chemistry data shows that 
bedrock mineralogy can be used to estimate the water chemistry of RMNP. 


Bedrock mineralogy will be used in later sections to assess the sensitivity 
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(as measured by alkalinity) of similar geologic areas in the central Rocky 
Mountain region. 


Alkalinity (Sensitivity) of RMNP Lakes and Streams 








An area's ability to neutralize acicic deposition through interactions 
with bedrock and surficial geology determines its vulnerability to acidic 
deposition. The lithological characteristics of a watershed, combined with 
‘ts water chemistry data, provide the basis for assessing the area's sensi- 
tivity to acid rain. Alkalinity, an integrator of watershed buffering mecha- 
nisms, 1s often used as an index of sensitivity. Hendrey et al. (1980) define 
sensitive waters as those with alkalinities lower than 200 yeg/1], a level low 
enough to be neutralized by acidic deposition and runoff. Using this defini- 
tion, we find that much of the RMNP is sensitive ta acidic deposition. 

As we have seen in the previous sections, the sensitivity of the water- 
sheds in RMNP is primarily determined by the bedrock geology of the water- 
sheds. Within each watershed, elevation may be used to further divide the 
basins into areas of different sensitivities. The classification scheme used 
for assessing the sensitivity of RMNP to acidic deposition defines sensitive 
waters as those with alkalinities less than 200 weq/!, while waters with 
alkalinities less than 100 peq/! will be considered very sensitive. Lakes and 
streams with alkalinities below 50 weg/! are classified as extremely sensi- 
tive. Following this classification scheme, we have evaluated the sensitivity 
above and below 3300 m in each subbasin (Figure 11). The results show the 
following classifications: 


Extremely sensitive Glacier Gorge 

(alkalinity < 50 peq/!) Ypsilon Creek 

Very Sensitive Roaring River 

(50 < alkalinity < 100 peq/1) Upper East Inlet 

Sensitive Lower East Inlet 

(100 yeq/! alkalinity < 200 peq/1) High elevations of 
Upper Fall River 

Nonsensitive Lower elevations of 

(alkalinity * 200 peq/!) Upper Fall River 


Upper Colorado River 


Current Acidification Status of RMNP 
Researchers have observed acidic deposition in the Rocky Mountains since 





1967 (Lewis and Grant 1979). These authors report a decrease in stream 
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Figure 11. Spatial distri 


bution of sensitivity im RMNP 
(see map Figure 1). 
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bicarbonate in the Como Creek watershed of north-central Colorado over a 
150-week period. This record, however, is too short to be considered evidence 
of increasing acidification. Lewis (1982) reports an average decrease in 
alkalinity from 22 mg/) to 18 mg/! as CO, [equivalent to a decrease from 1000 
yeq/!] to 818.2 wyeq/! (Lowenthal and Marias 1978)] for 104 lower-elevation 
lakes surveyed by Pennak from 1938-1942 and resurveyed by Lewis in 1979. 
Although Lewis attributes this loss to acidic deposition, the concurrent 
decrease in total residue suggests that hydrologic variability is responsible 
for this decrease in alkalinity. The high discharge in 1979 (~30% above 
average) and the similar decrease in total residue and alkalinity (22% for 
alkalinity, 21% for total residue) suggests that this decrease is mainly 
caused by dilution. In addition, a 180 yeq/! loss of alkalinity is much 
greater than any alkalinity loss observed in acidified regions of the north- 
eastern United States (Hendrey et al. 1980). Although several areas in the 
Rocky Mountains have been shown to be potentially sensitive to acidic deposi- 
tion (Harte et al., submitted; Dodson 1981; Turk and Adams 1983), no con- 
vincing evidence of acidification has been presented. 

At present, RMNP is receiving much less acidic deposition than acidified 
regions in the northeastern United States. A comparison of the (volume- 
weighted average) concentrations of the major ions in precipitation show that 
the precipitation at RMNP contains much less acid than that at Hubbard Brook 
(Table 16) (NAD? 1982). During 1981 the average precipitation pH at RMNP was 
5.07, compared with an average of 4.33 at Hubbard Brook. Sulfate concentra- 
tions at Hubbard Brook are also higher than those observed at RMNP (48.9 peq/! 
at Hubbard Brook, compared with 35,0 peq/! at RMNP), while the Cp concentra- 
tion at RMNP is more than double the concentration at Hubbard Brook. However, 
the nitrate concentrations observed at these two locations are very similar, 
21.9 yeq/| at RMNP, compared with 21.2 yeq/!| at Hubbard Brook. These results 
Suggest that much of the acidity observed in RMNP precipitation may be caused 
by nitrate, while the sulfate is probably a combination of anthrorogenic 
emissions of sulfur dioxide, along with sulfate associated with Cp from air- 
borne dust and salts. 

Since no historical water-chemistry data are available for RMNP, we base 
our assessment of its present acidification status on the current composition 


of its waters. According to electroneutrality conditions: 
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Table 16. Comparison of Rocky Mountain National Park 
and Hubbard Brook precipitation chemistry (Nationa! 
Atmospheric Deposition Program 1982). 

















RMNP Hubbard Brook 
Vol. wt. Yearly Vol. wt. Yearly 
average deposition average deposition 
lon (yeq/1) (meq/m- ) (yeq/1!) (meq/m ) 
HK 8.6 3.0 46.5 71.1 
so, 35.0 12.2 48.9 75.8 
NO,” 21.9 7.5 21.2 32.5 
cv? 4.9 1.7 5.3 8.3 
+2 
Ca 21.9 7.7 6.5 10.4 
Mg’? 8.1 2.9 3.6 5.6 
x*) 5.8 il 0.4 0.7 
+] 
Na 30 2.1 5.1 7.9 
7 | "2 aml -1_ tl | +n 
Cp = HCO, + 50, + NO, + Cl H metals 


Under acid rain conditions, anthropogenic acids (H,SO, and HNO, ) reduce the 
relative concentration of alkalinity in this equation by titration. We can 
use the “excess” concentration of base cations to provide an estimate of the 
alkalinity replaced by acid anions (i.e., amount of acidification). After 
correcting for the deposition of salts and excluding H” and metals, which are 
negligible at the pH of water in RMNP, the ion balance equation is reduced to: 
C,* - HCO, = C,* ; where C,* = NO, + S0,°* . 

B 3 A A 3 4 

The results of these calculations for the most sensitive watersheds 
(Clacier Creek, Loch Vale, East Inlet, and Ypsilon Creek) show that these 
waters have, at the most, suffered a smal! loss of alkalinity (Table 17). It 
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Table 17. Calculation of acidification status for Rocky Mountain 
National Park watersheds according to ion balance considerations. 








c _ Calculated 
Watershed ; “(eea/l) = Acidification 
(yeq/1) ” (yeq/1) 

Roaring River 12.3 74.4 Nore 
Ypsilon Creek 56.2 48.0 8.2 
Upper Fall River 155.2 180.5 None 
Andrews Creek 57.3 38.8 18.5 
Loch Vale 47.9 41.2 6.7 
Glacier Creek 47.6 40.3 7.3 
Upper Colorado River 280.1 319.0 None 
East Inlet 91.5 85.5 6.0 





is so small that the combined effect of analytical error and the error 
associated with the correction for atmospheric salts is probably as large as 
the calculated acidification. The Upper Colorado River, Upper Fal! River, and 
Roaring River show no evidence of acidification. 

A similar method for evaluating the current acidification status of low 
ionic strength lakes is the “predictor nomograph" developed by Henriksen 
(1980). This mode] was empirically derived to evaluate the impact of in- 
creased acid precipitation on lakes in Norway but may also be used to assess 
current acidification status. The model is based on electroneutrality condi- 
tions, assumes no increase in weathering of base cations, and considers 
atmospheric deposition to be the only source of $0;°. The nomograph is a plot 

e and Mg*? vs. “excess” so; 
means monmarine in origin. The lake 504° concentrations have been corre! ated 


of the “excess” Ca’ in lake waters. “Excess” 


with pH and $0, concentrations in rain by a regression analysis of data from 
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719 Norwegian lakes. The graph has been divided into three sections to 
represent three stages of acidification. The first stage represents waters 
that stil! contain enough alkalinity to buffer incoming acid precipitation, 
“bicarbonate lakes." The second stage of the nomograph is representative of 
“transition lakes." These waters experience rapid fluctuations in pH because 
of their low alkalinities. The third stage of acidification are the “acid 
lakes", characterized by low pH and increased aluminum concentrations. The 
nomograph is used to predict acidification status by plotting precipitation pH 


or excess jake s0,°* with “excess” ca’? ¢ Mg*¢ 


in peaq/l. This simple 
empirical mode] has been successfully applied to lakes in the Adirondack, 
Canada, Scandinavia, and Scotland (Wright et al. 1980). Although there are 
some differences between these areas and RMNP (notably the soils), the simi- 
larity of geochemical processes should allow our use of the nomograph with 
some changes. 

We have modified the Henriksen nomograph to use as a too! to evaluate the 
current and future effects of acid rain in the Rocky Mountains. ca’? and Mg"? 
are the major base cations in Henriksen's study area, while na”? + Mg’? are 
the dominant cations in RMNP. As a result, we have replaced the “excess” 
ca’? ¢ tg"? used by Henriksen as the y axis in the nomograph with Cp". As 
previously stated, RMNP receives minimal sea-salt deposition but does receive 
atmospheric deposition of salts and dust. Cp" is our best estimate of 
“excess” base cations. 

Since much of the acid deposition in this region is in the form of HNO, 
(NADP 1982; Lewis and Grant 1979; Kelley and Stedman 1980), we have further 
modifed the Henriksen nomograph by using C,* (s0"* + No; *) as the x axis in 
the nomograph. The waters in RMNP contain a large indeterminate “natural,” or 
background, concentration of 04°, a result of the atmospheric deposition of 
dust and salts. We expect that only part of the 50,” in RMNP lakes is the 
result of acid deposition. Almer et al. (1978) estimate the background 
concentration of 0, in Scandinavian lakes to be 20-60 peq/l. To estimate 
acidic deposition, we will use Almer's lowest estimate of background s0,° to 
calculate a maximum “excess” 504°": 


2 


04°" = $0, = 20.0 yeq/] 


; : Qs a "2s 
In this scenario, Cy = NO, + $0, . 
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Plotting this information on the nomograph shows that most of the lakes 
in RMNP can be classified as “bicarbonate”, while a few approach the 
“transition” stage (Figure 12). 
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Figure 12. Status of 23 lakes in RMNP. 


The high cation concentrations in RMNP precipitation indicate that the 
park is undoubtedly receiving “natural” deposition of sulfate in the form of 
salts and dust from the surrounding arid regions. However, the low pH of 
Rocky Mountain precipitation indicates that some of the sulfate and nitrate is 
being deposited as acid deposition. It is likely that the waters in RMNP have 
suffered some loss of alkalinity because of acidic deposition, but this loss 
(a maximum of ~10 yeqg/!) is minor compared with that experienced in lakes in 
the northeastern United States. 

Future Impacts of Acidic Deposition 

As we have seen in the previous sections, the Central Rocky Mountains are 
extremely sensitive to acidic deposition but have as yet suffered little or no 
acidification on a regional basis. An increase in acidic deposition could 
have some serious consequences. The development of major coal and oj! shale 
resources upwind of this region and the subsequent increase in acidic deposi- 
tion could have a serious impact on the surface water acidity in the region. 
Using the Henriksen predictor nomograph described in the previous section, we 
can estimate the best- and worst-case consequences of an increase in the 
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current rate of deposition, to that experienced in the northeastern United 
States. The scenari. for increased acidic deposition in RMNP assumes a 
decrease in precipitation pH to that experienced in the northeastern United 
States. This involves an tncrease in H” from pH 5.1 to 4.3, which wil! be 
accompanied by an increase in lake Cc." of approximately 80 yeq/!. 

In the worst-case scenario, the increase in acidic precipitation wi!! not 
increase <p” but WINE only or tah in decreased alkalinity. Using our estimate 
of current c,". (S0,°* ¢ NO, *"- 2( peq/!) a worst-case prediction using the 
nomograph indicates that most of the lakes in RMNP will reach “acid” status, 
while the remainder will be classified as “transition” lakes (Figure 13). 

In the best-case scenario, the increase in acidic deposition will be 
accompanied by an increase in Ce” of 0.4 wyea/! per 1.0 yeq/! c,". This 
increase in Ce" with increased acidic deposition was empirically calculated by 
Henriksen using data from low-alkalinity lakes in several areas The best- 
case scenario predicts that a few lakes wil! remain “bicarbonate,” the 
majority wil! become “transition” lakes, and several lakes wil! stil! reach 
"acid" status (Figure 13). 

Henriksen's evaluation of data sets from Norway, Sweden, Canada, and the 
United States shows that there is an increase in base cation weathering with 
acid rain for some, but not al] watersheds. Increases in Cp will probably 
occur in regions that have soils with high levels of exchangeable bases. When 
these areas receive acidic deposition, base cations on the soi! exchange sites 
will be replaced with H° by mass action, increasing surface water Ce”. In 
regions with poorly buffered soils, the increase in C" with acid deposition 
will be smaller, since fewer exchange sites contain base cations for ton 
exchange. Other researchers have shown that the rate of primary mineral! 
weathering in areas with granitic and metamorphic rock does not increase under 
acid rain conditons (Johnson et al. 1981). The low ion-exchange capacity of 
the soils in RMNP and the bedrock composition indicate that the increase in 
Ce" with increased acidic deposition in RMNP will be low and that with in- 
creasing acidic deposition many of the lakes will shift from bicarbonate 


towards acid status. 
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Figure 13. Predictor nomographs projecting the best and worst-case 
estimates of the effects of increased acidic deposition on lakes in RMNP. 











YELLOWSTONE NATIONAL PARK 


INTRODUCTION 


Located in northwestern Wyoming at the northern edge of the potentially 
sensitive Rocky Mountain region, Yellowstone National Park (YNP) is charac- 
terized by extreme variability in geology, geologic history, and water chem- 
istry. The 3,742 square mile park contains four large lakes (Yellowstone, 
Shoshone, Lewis, Heart), and numerous streams, rivers, and smal! back-country 
lakes, 4 number of which support trout populations. As in the case of Rocky 
Mountain National Park, the general objectives of the study in Yellowstone 
National Park were to determine the sensitivity of waters to acidification by 
current levels of acid precipitation and to determine whether this will affect 
fish populations. The accomplishment of these objectives is based primarily 
on the evaluation of a water chemistry data base developed over the last 30 
years. 


DESCRIPTION OF GEOLOGY 


Topographically, Yellowstone is dominated by a high plateau from about 
1800 to 2800m in altitude. The plateau is bounded on three sides by 
mountain ranges: the Gallatin and Beartooth mountains to the north, the 
Absaroka range to the east, and uplands from the Teton and Washakie ranges in 
the south (Cox 1973). The continental divide transects the park along the 
southwestern edge--approximately 20 percent of the park area drains west into 
the Snake River basin, and 80 percent drains east into the Missouri River 
drainage. 

Cox (1973) divided Yellowstone into seven hydrologic units following 
geologic and geographic boundaries. These are: Rhyolite plateau (RP), 
Gallatin (GT), Beartooth (BT), Absaroka (AB), Falls River (FR), Snake River 
(SR), and West Yellowstone (WY) (Figure 14). The plateau area (RP) is under- 
lain by Tertiary and Quaternary rhyolite flows, while north and south of the 
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Figure 14. Regional-geological map of Yellowstone National Park (adapted 
from Cox 1973). (AB) Absaroka region is predominantly andesitic lava 
flows and breccia, with basalt, and some occurrence of rhyolite, sandstone, 
and limestone; (B87) Beartooth region is a mix of Precambrian granites, 
Paleozoic and Mesozoic sandstones and shales, and Tertiary/Quaternary 
volcanics; (FR) Falls River region is Quaternary rhyolite and basalt, 
frequently overlain by alluvial and glacial deposits; (GT) Gallatin 
region is Precambrian granites, Paleozoic and Mesozoic limestones, 
sandstones, and shales, Tertiary/Quaternary volcanics; (RP) Rhyolite 
plateau region is predominantly Tertiary and Quaternary rhyolite flows; 
(SR) Snake River region is Paleozoic and Mesozoic limestones, sandstones, 
and shales, with some mg pe be of Tertiary rhyolite and andesite; (WY) 
West Yellowstone region is rhyolite overlain by -lluvial, glacial, and 
lacustrine deposits. 


plateau are a heterogeneous, heavily faulted mix of PreCambrian granites, 
Paleozoic and Mesozoic limestones, shales, and sandstones, and more recent 
(Tertiary, Quaternary) volcanic deposits (GT, BT, SR). Between the plateau 
region and the Gallatin range along the extreme western park border lies a 
series of deep (c.200 feet) alluvial and glacial deposits underlain by 
volcanic rocks (WY) East of the plateau is a mountainous area composed 





mostly of Tertiary andesitic lava flows and breccia (AB), while basalt flows 
occur In the extreme southwest (FR). 

Rhyolite, located in the southwest and central portion of the park as 
well as im outcroppings throughout the rest of the region, is clearly the 
dominant bedrock type. Extremely rich in silica, rhyolite is chemically 
nearly equivalent to granite (Bryan 1979). Although fairly uncommon through 
the world, nearly all geysers are associated with rhyolite formations, since 
the rock provides the pressure-tight “piping” necessary for geyser creation. 
The rhyolite in Yellowstone consists primarily of lava and welded tuff with 
assorted deposits of breccia, ash, and glass. All of the types are chemically 
similar, although tuff may be slightly more weatherable than lava due to 
differing extrusion and deposition procecses (Herzog 1982; Cox 1973). 
Rhyolite may be altered by hot water and gas near hot springs. 

Weatherability of the bedrock types may be roughly ranked by chemical 
Characteristics as follows (Herzog 1982 and personal communication; Loughnan 
1969): 


Ae granite, rhyolite 
ash flow tuff 


metamorphics--biotite gneisses and schists 





andes ite 
y basalt 
high limestones, shales, sandstones 


In addition to a heterogeneous geology, the park bears the ubiquitous remains 
of three major glaciations originating in “he Absaroka Range and Beartooth 
Mountains. Most of the current park deposits date to the last major glacia- 
tion, the Pinedale (10-30,000 years b.p.), which covered 90% of the park. 
Abundant till and kame deposits throughout Yellowstone attest to the influence 
of these geologic events 


DETERMINATION OF SURFACE WATER CHEMISTRY 


Surface water quality studies in Yellowstone date to as early as the 19th 
century, when Gooch and Whitfield (1688) published the first chemical data for 
Yellowstone Lake. Until! comparatively recently, however, surveys concentrated 
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on the four largest lakes, all of which are chemically insensitive to acidi- 
fication by acid precipitation (alkalinity »200 weq/!). Chemical surveys of 
the more vulnerable back country lakes were initiated in 1964 by the 
Yellowstone Nationa) Park Fishery Management Investigation, and continue to 
the present. 

Field pH and alkalinity were measured with a Hach Field engineers kit 
until! 1974, when the Hellige kit was introduced for pH measurements. During 
1969 and 19/0, ‘laboratory analysis of water samples were performed by the 
Bureau of Indian Affairs Soi! Laboratory in Gallup, New Mexico. (No labora- 
tory analyses were performed from 1971 to 1973.) Beginning in 1974 and con- 
t'nuing to the present, lab analysis of samples was performed by Orlando 
Laboratories, Orlando, Florida. Almost al! samples were refrigerated, mailed 
within two days, and processed within a week. Analytical methods used by both 
labs were in accordance with then-current Standard Methods for the Analysis of 
water and Wastewater. 

All pH values reported in this study were obtained colorimetrically in 
the field. Burns et al. (1981) foun’ that measurements from a Hellige kit 
agreed to within 0.15 pH unit with potentiometric values, although Pfeiffer 
and Festa (1980) report a systematic bias expressed by the followirg 
relationship: 


Old (Hellige) = 0.6639 (pH meter) + 2.534 


re = 0.9) 


All alkalinities reported were determined in the laboratory by cclori- 
metric titration to pH 4.6 (Standard Methods #403), except during the years 
1965-1966 and 1971-1973, when only field measurements were taken. Field 
alkalin‘ties were determined colorimetrically during this period with the Hach 
field engineers kit (Model DR-EL or AL-36-P). Colorimetric titrations do 
overestimate alkalinity by the amount of free hydrogen ion in solution at the 
endpoint; at titration to pH 4.6, this is 25 weg/! alkalinity for all samples. 
Analytical methods for major cations and anions are given in Table 18 for the 
years 1969-1970 and 1974-present The current (15th edition) Standard Methods 
is referenced in Table 18, although it its assumed that the most recently 
available version was followed for historical samples. Quality assurance was 
determined by calculating ton balances for al! those lakes with complete 











Table 18. Analytical methods, Yellowstone Nationa! Park. 





Parameter Method, 1969-70 Method 1974-present 





Alkalinity Standard method #403. Colorimetric titration using 
phenciphthalein (end point pH 8.3) and methy! orange (end 
point pH 4.6) indicators 


Calcium Standard method #306C Hach Method, p. 122 
EDTA titrimetric Titrimetric 

Magnes ium Standard method #313C; calculation from EDTA hardness 

Sodium Standard method #325B; flame photometric 

Potassium Standard method #322B; flame photometric 

Sulfate Standard method #427A EPA #375.4 
Gravimetric with ignition Turbidimetric 
residue 

Chloride Standard method #4074; argentometric 





(Standard methods 1980; EPA 1979; Hach 1978) 


chemical records (see Appendix 0D). Lakes with ion balances with absolute 
values less than or equal to 20% are considered usable for the purposes of 
this study. 

The 106 lakes discussed in this report represent a significant sample of 
all the lakes in YNP. All the major lakes have bee. surveyed, including those 
of special recreational or scientific significance. Every region, major 
geological formation and geochemical type of lake has been sampled in the 
survey. Because of its large size, Yellowstone Lake was surveyed at four 
different locations in the lake. The chemistry of the lakes not surveyed may 
be inferred by the chemistry of neighboring surveyed lakes. 


RESULTS AND DISCUSSION 


As stated in the introduction, the ®4AS project relied on available data 
from Yellowstone National Park. There was thus no opportunity to design a 
sampling program to test hypotheses relating lake and stream sensitivity to 


elevation, soils, and geology, as was the case for Rocky Mountain Nationa! 
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Park. YNP data were available for lakes only, and there had been no attempt 
to obtain data on complete watersheds or to relate sampling locations to 
elevation, soils, or geology In addition, water chemistry data were obtained 
over the years using different analysis methods and samples were analyzed by 
different laboratories. Al! of these factors made it impossible to carry out 
the extensive evaluations reported in the chapter for Rocky Mountain Nationa! 
Park. The evaluation of information from YNP focuses on sensitivity or 
alkalinity, which are the critical factors in assessing the health of fish 
populations. In addition, the general YNP findings can be compared with those 
from RMNP to verify the ability to extrapolate resuits from one region to 
another. 


Lake Alkalinities (Sensitivities) 





A spatial sensitivity map, using alkalinity as an index of vulnerability 
to acidification, was created for the park with recent alkalinity data 
(Figure 15). The selection of alkalinity as the best index of sensitivity to 
acid precipitation is based on its physical significance as the emergent sum 
of many acid-neutralizing processes occurring in the watershed, and its wel!- 
documented relation to pH (Henriksen 1979). Using Hendrey et al.'s (1980) 
convention, those lakes with alkalinity values <200 weq/! are considered 
potentially sensitive to long-term inputs of acid precipitation. For 
Yellowstone, lakes with alkalinity values reported as 230 or less are con- 
sidered in this category, due both to the underestimation implicit in colori- 
metric titrations, and to account for imprecision in the analysis. One fourth 
(30 of 106) the surveyed lakes may be considered “sensitive” by this crite- 
rion. Thirteen of these lakes show alkalinities 5100 peq/!, although at least 
Six are influenced by thermal springs or humic acids (YNP Fishery Management 
Investigations reports). 

Excluding those lakes that are located in the midst of major geotherma| 
areas, the lowest alkalinity measured was 40 yeq/!. Yellowstone iakes on the 
average apparently show considerably greater levels of alkalinity than those 
determined for Rocky Mountain National Park lakes, probably due primarily to 
differing geology between the two areas (most of the Rocky Mountain Nationa! 
Park samples are collected in areas underlain by Precambrian granites and 
metamorphic rocks). Of the four basins considered in the RMNP study, the 


Upper Colorado River watershed is geolugically the most similar to Yellowstone, 
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Figure 15. Alkalinity map, Yellowstone National Park. 


Yellowstone, as it shows a significant (ca. 8%) accumulation of volcanic 
rocks. This watershed shows the highest alkalinity and pH values of the RMNP 
Study sites [average alkalinity equals 248 peq/! (this study) ]. 

Although the low-alkalinity (<230-peq/1) lakes in Yellowstone show no 
clear geographical (spatial) patterns, a geological] pattern is evident. Most 
occur within the large rhyolite flow which rises from the southwest and 
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spreads along the central-west and centra! portions of the park. Lakes with 
slightly higher alkalinity values (ca. 300 peq/’} are found in the extreme 
Southwestern Fall River basalt formation and the northwestern part of the 
rhyolite plateau. Most of the northernmost lakes, dominated by andesites ana 
basalts, show very high alkalinity (>1000 yeq/]) and are not sensitive to acid 
deposition. The exception to this is a group of five high-altitude lakes of 
the Specimen Creek drainage in the northwestern corner of the park. These 
smal! headwater lakes have an andesite-basalt bedrock but show very low 
alkalinities (<200 peq/1l). The water supply of these lakes is primarily from 
snowmelt (YNP 1965-1981)) and lakewater chemistry appears to be influenced 
more by this dilute source than by reaction with the surrounding bedrock. 

field pH measurements were used to create a spatial pH map to complement 
alkalinity in determining geographical distributions of sensitivity 
(Figure 16). The pH map for Yellowstone shows most lakes with circumneutral 
pH's of 6.5-7.49, and most of the rest of the lakes siightly to strongly 
alkaline. Six of the 106 lakes have pH values less than 6.5. Of these six, 
two are unquestionably influenced by thermals and one probably so, and three 
are dystrophic lakes. As with the high alkalinity lakes, the greatest pro- 
portion of high pH lakes is found in the northern part of the park. 

A similar map denoting sulfate concentrations for the Yellowstone lakes 
(Figure 17), shows a cluster of high sulfate lakes in the north and central 
east areas; these lakes are generally characterized by high alkalinity (>1600 
yeq/!) and in some cases are found in low-lying marshy areas. Sulfate in 
these lakes probably has a biological and/or bedrock component § (YNP 
1965-1981). Sulfate concentration in lakes throughout the rest of the park 
Shows no particular trends; most lakes show sulfate values <200 yeq/!, while a 
majority of lakes in the rhyolite bedrock have sulfate concentrations <100 
yeq/1. 

The alkalinity, pH, and sulfate maps for Yellowstone must be viewed with 
a number of factors in mind. First, all lakes for which "recent" (1965 or 
later) data were available were used to create the maps, regardiess of ion 
balance. Some of the values, therefore, may be in question from a purely 
analytical view. Second, there exists a l6-year spread in chemical analyses 
for the maps: Lakes sampled once in 1965 are not distinguished from those 
sampled in 1980. In al] cases, the most recent chemical measurements are used 
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Figure 16. pH map, Yellowstone Nationa] Park. 


for each lake. Third, the water sources for the lakes vary, and reflect 
different residence times through different watersheds. Only two of the lakes 
with alkalinities <230 yweq/! have significant inlets, and most lack outflow. 
The primary source of water for these dilute lakes is groundwater and snowmelt 
(P Gresswel], personal communication). A number of the less sensitive 
la es have significant flow-through. Finally, humic lakes, geothermal ly 
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Figure 17. Sulfate map, Yellowstone National Park. 


influenced lakes, and those fed by subterranean seeps through calcareous 


deposits are not distinguished from those lakes which more readily reflect the 


influence of the surficial geology. Indeed, many lakes in this geologically 


active area may be influenced by unknown factors. Figure 18 shows the loca- 


tion of the major thermal springs and geysers in the park (Waring 1965): 
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Figure 18. Location of hot springs and geysers, Yellowstone National 
Park (from Waring 1965). 


springs and geysers in the park (Waring 1965): These are concentrated in 
geyser ybasins along the Firehole and Gibbon rivers, as wel] as north of 
Yellowstone Lake. Major hydrothermal activity in the park generally follows 
the southwest-northeast sweep of the rhyolite plateau. 

Examination of alkalinity and pH data indicates that most of the lakes 
are wel! buffered against potential increases in rainfall] acidity. One fourth 
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of the surveyed Yellowstone lakes do report alkalinity values <230 wpeq/! 
(defined as “sensitive” in this report), although a number of these (at least 
7 of 30) receive significant interna! acidity from humic acids or hydrothermal 
springs. All but four of the lakes in the < 230 peq/! class are naturally 
barren of fish. 

Factors Influencing Alkalinity (Sensitivity) 





Numerous investigators (e.g., Galloway and Cowling 1978) have considered 
surficial geology as the critical factor in determining water chemistry of 
many lakes and streams. Commonly, surficial geology is correlated with water 
chemistry when soils are derived from the underlying bedrock, surface water 
originates as flow of incident precipitation through a watershed, and chemical 
constituents are derived primarily through ion exchange/mineral weathering 
reactions occurring in the soil. 

As a whole, lakes in Yellowstone are influenced by far more complex 
factors. ‘hree major glaciers caused the deposition of extensive al locthonous 
rubble in some areas. The chemical characteristics of this transported rock 
material may, in some cases, dominate lake chemistry. The ash and lava flows 
covering iarge areas of the park are dotted with hydrothermals which provide 
internal sources of sulfate and other chemicals. Underground springs, 
especially in the north, may contribute large concentrations of dissolved ions 
to lakes in “unreactive” bedrocks. Finally, watersheds frequently show the 
effects of a number of different geologic events, making simple separation of 
most lakes into “dominant” bedrock type subject to error. 

Considering only lakes with good ion balances and no known humic or 
thermal influence, a separate variance t-test was used to determine if signif- 
icant differences existed in surface water chemistry between the volcanic and 
the mixed-geology regions. The results of the pairwise comparison indicate 
that the alkalinity and th. base cation sum are significantly different 
between the volcanic (FR, RP) and the older sedimentary-volcanic (GT, AB, BT) 
regions of the park, with p < 0.01 (Table 19). Both alkalinity and base 
cations are significantly lower in the volcanic regions. There appears to be 
no statistically significant difference in lake sulfate concentration between 
the two regions, indicating that local sulfur sources are not confined to a 
Single region, but may influence lakes in many areas of the park. 

The statistical results indicate that, with ihe best available separation 
of “poor” data (bad ionic balance, thermal or humic influence), lakes in the 
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Table 19. Regional alkalinity of Yellowstone Nationa! Park lakes. 





Mean alkalinity 





Region (yeq/1) No. of samples Range (yeq/!) 
rp? 182 3 80- 360 

FR 320 2 320 

GT 1468 8 640-2280 
AB 1493 3 240-3520 
BT 1907 3 160-4800 





[FR = Ralls River hydrologic unit, RP = Rhyolite Plateau hydrologic unit, GT = 
Gallatin hydrologic unit, AB = Absaroka hydrologic unit, BT = Beartooth 
hydrologic unit. 


Rhyolite Plateau and Fall River basalt region (i.e., those most influenced by 
volcanic bedrock) are potentially more sensitive to acidification than those 
in the non-volcanic northern and eastern ranges. The term “sensitive” is, as 
always, relative, for the more sensitive lakes in Yellowstone show similar 
geology to the least sensitive lakes in Rocky Mountain National Park (this 
Study). 

Alkalinity is commonly observed to decrease with increasing elevation, as 
soi! depth and development, and watershed residence time are decreased. A 
fairly strong alkalinity-eievation correlation exists for a number of water- 
sheds in Rocky Mountain National Park, with only a weak trend for the 
volcanic-andesite dominated Upper Colorado watershed (Figure 8). Results for 
the Yellowstone lakes dominated by rhyolitic becroc' are intermediate 
(Figure 19), showing a fairly strong correlation, except for two lakes located 
in deep depressions at the edge of the formation. The Yellowstone rhyolite is 
a somewhat more homogeneous bedrock than the Upper Colorado geology; this may 
explain the clearer gradient. 

An evaluation of the relationship of sensitivity to basin characteristics 
was attempted by statistical analysis of water quality data of lakes from 
differing hydrologic regions. Lakes in regions dominated by volcanic bedrocks 
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Figure 19. Alkalinity versus elevation, lakes draining rhyolite bedrock, 
Yellowstone Nationa! Park. 


show the lowest alkalinities and base cation concentrations. On the whole, 
these lakes may be considered the most sensitive in the park to acidification. 
A wide spread in chemistry values for the other lakes indicates that deep 
springs or local geologic deposits may greatly influence water chemistry. 
Historical Changes in Surface Water Chemistry 

Acidification trends may be determined in two major ways: Analysis of 
historical data, and use of (empirical or mechanistic) predictive models. The 
major difficulties with historical data are that (1) water chemistry may 
change daily and seascrally in relation to biological activity, (2) method- 
ology may change over the years, and (3) hydrological and meteorologica! 
conditions may not be comparable. Minimization of daily and seasonal fluctua- 
tions is best accomplished by holding these parameters as constant as 
possible, and by comparing mean values of a number of lakes. A! water 
Chemistry data used in this study are from the summer months of June through 
September, but daily fluctuations in the chemistry were not determined, since 
lakes were sampled at different times of the day from one year to the next 
Changes in methodology may be accounted for by application of “correction 
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factors” if earlier methods are systematically biased (e.g., Burns et a). 
1981; Pfeiffer and Festa 1980). When complete chemical analyses are avail- 
able, ion balances may serve as a quality check. For the Yellowstone data, 
all alkalinities were determined using colorimetric titrations, and ai! pH's 
with either the Hach or Hellige kit. Since a single lab performed analyses of 
Yellowstone Park water since 1974 and methods have mot changed considerably, 
data collected since that year are hig ly comparable 

The final difficulty with historical data is related to differing hydro- 
logical and meteorological conditions. Tonic concentrations may be substan- 
tially reduced by dilution during a wet period and increased by concentration 
during a dry period. While absolute concentrations have changed, relative 
concentrations (with respect to the sum of al! tonic constituents) have not. 
Figures 20 through 23 attempt to account for simple hydrologic variations in 
comparing historical with recent data. The dotted line in each figure 
connects the origin with the ion sum for each year. The position of each 
individual ion with respect to that line indicates whether its relative con- 
centration has increased (above the line) decreased (below the line) or 
remained constant over the period of time indicated (from Henriksen 1982). In 
Figures 20 through 23, a solid 45 degree line is included for comparison. By 
indicating relative changes in surface water concentration, the plots reduce 
interpretive error due to differences in hydrologic conditions. Figures 20 
through 23 represent available historical data for lakes with alkalinity 
values of <500 yeqg/!, including those with poor ion balances (ji.e., »°20% 
variations). 

Of the seven lakes shown in the historical comparison plots, (three show a 
relative decrease in alkalinity (High, Crescent, Grebe) and increase in 
sulfate over the 10-year period, two show the reverse trend (Wolf, Ice), and 
two show no change in alkalinity (Crag, Cascade). Thus, the overal! trend for 
the seven lakes appears random. Iwo of the lakes showing alkalinity weclines 
are located in the extreme northwestern corner of the park in a group of five 
small snowmelt-seepage lakes in the Specimen Creek watershed. These lakes 
show uniformly low alkalinities; however, ion balances are poor. More 
sampling would be necessary to make definitive statements about the sensi- 
tivity of these lakes. Of the three lakes which show an increase in lake 
sulfate, the relative sulfate concentration changes (averaging 80 weq/!) are 
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Figure 20. Historical chemistry comparisons, High and Crescent Lakes. 
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Figure 23. Historical chemistry comparisons, Cascade Lake. Al] 
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considerably greater than the sulfate load in precipitation (approximately 33 
yeq/|). These sulfate changes are probably too great to be caused hy acid 
precipitation and likely reflect the influence of local natural sources of 
sul fate. 

The relative concentrations of sulfate and alkalinity were also computed 
for the seven lakes and are shown in Table 20. Comparing data in terms of 
relative concentrations, or percent of each ion in reiation to the total ion 
sum, 7S another method to minimize the influence of sampling under different 
hydrologic conditions. t-tests on historical versus recent relative concen- 
trations of alkalinity and sulfate for the seven lakes show no significant 
changes in these parameters over time. 

Predictive models, such as that developed by Henriksen (1979), generally 
assume that bicarbonate lost in water acidified by atmospheric deposition is 
Stoichiometrically replaced by sulfate, as strong acids from anthropogenic 
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Table 20. Relative concentration® of alkalinity and 
sulfate, lakes with historical data. 














Alkalinity Sulfate 

Lake 

Historica] Recent Historical Recent 

(1963-70) (1978-80) (1963-70) . (1978-80) 
Crescent 47.2 27.8 11.3 25.2 
High 36.1 20.1 8.4 20.9 
Cascade 33.3 30.0 12.7 16.9 
Ice 18.6 36.2 28 8 9.1 
Crag 40.2 41.2 24.4 10.7 
Wolf 27.5 35.6 17.3 12.3 
Grebe 45.3 35.8 5.1 13.2 
x 35.5 32.4 15.4 15.5 





3Relative concentration = [ion of interest]/ion sum x 100. 


sources titrate existing alkalinity or replace bicarbonate as a major weather- 
ing ion. This assumption implies a negligible increase in base cation release 
from soils with acidification. A second assumption is that "internal" natural] 
sources of sulfate are negligible. However, Figures 20 through 23 indicate 
that historical changes in alkalinity are generally balanced by changes in 
both calcium (cation compensation for acidity increases) and sulfate, in 
violation of the first assumption. The second assumption is also invalid for 
Yellowstone, which is heavily influenced by numerous geothermal sources of 
Suifur. This may be illustrated by comparing the concentrations of chloride 
and sulfate among the Yellowstone lakes: chloride (Figure 24), assumed to be 
a conservative ion, remains relatively constant in nearly all the Park lakes. 
In this plot, four of the five outliers represent large, old lakes (three 
samples from Yellowstone Lake and one from Lewis Lake); these are probably 
influenced by long-term concentration by evaporation. In contrast, sulfate 
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CHLORIDE VS. ALKALINITY 
YELLOWSTONE NATIONAL PARK LAKES 





258 — 
~ 
- 
“ a 
‘ 
20e — 
S 1 
¢ : . 
_ - 4 
= 150- 
* 4 4 ° 
— 4 
= — 
oOo 
2 ‘ 
So 100- 
« 
aa, 4 
: a. 8s a4 Sa 44 4 a 4 
:. rN a 
_ 





eprrenrerryrrerperernr rrr 


50@ 75@ ieee 1250 1500 1780 


ALKALINITY IN vea/1 


Figure 24. Chloride versus alkalinity, Yellowstone National Park lakes. 


concentration throughout the Park lakes is highly variable (Figure 25). In 
Figure 25, a cluster of 17 low sulfate, low alkalinity lakes (mean sulfate 
concentration = 35.8 peq/1) is evident in the Park, along with a wide scat- 
tering of higher sulfate lakes. Thus, the unique characteristics of 
Yellowstone make use of a Henriksen-like predictive acidification mode] 
inappropriate. 

Detecting historical changes in surface water chemistry for Yellowstone 
is confounded by natural variability in water chemistry, internal sources of 
acidity, and a scarcity of reliable historical data. Yellowstone precipita- 
tion is currently not highly acid, nor do the lakes appear acidified by 
anthropogenic activity. No significant overall] change in water chemistry was 
noted from lakes sampled once, then resampled several years later. The dilute 
lakes in the rhyolite plateau and also in the far northwestern corner of the 
Park may be the most subject to future acidification. 
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Figure 25. Sulfate versus alkalinity, Yellowstone National Park lakes. 


Current Status and Future Trends in Surface Water Chemistry 

Current, historical, and future acidification scenarios depend upon 
changing composition of precipitation. Yellowstone receives an average of 
42 cm of rainfall annually (NOAA 1970), but, because of the Park's great size 
and varied terrain, precipitation distribution is not homogeneous. The 
Continental Divide intercepts eastward moving storms from the Pacific, causing 
the greatest amount of precipitation near the Divide in the southwest, and the 
lowest in the low altitudes. Up to 1.2 m of snow may accumulate in the Park 
on average winters. Figure 26 indicates how weekly averaged pH and sulfate 
concentrations over the NADP deposition-monitoring site vary with season from 
June 198) to January 1982. Although no overall] trends are evident, chemistry 
is quite variable, with a few sharp sulfate peaks over the year. These high 
peaks represent the lowest volume periods of the record, and are not corre- 





lated with especially low pH. 
In 1981, weighted average precipitation pH at the NADP collecting station 
was 5.2 and weighted average annual sulfate concentration was 33.5 peq/). 
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PRECIPITATION CHEMISTRY 
YELLOWSTONE NATIONAL PARK 
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Figure 26. Precipitation pH and sulfate concentration, Yellowstone National 
Park, WY, 1980-1981 (from NADP 1981). 


During the same year, over a “typically” acidified area of the eastern US., 
Hubbard Brook, New Hampshire, precipitation pH averaged 4.4 and precipitation 
sulfate averaged 48.9 yeg/!. A comparison of these two sites indicates that 
the northwestern Wyoming area is currently not threatened by highly acid 
precipitation. 

The major anionic contributor to strong acidity is the sulfate ion. 
Although sulfate mobility can be affected by soil adsorption and biological 
uptake (Johnson and Cole 1980), inputs and outputs of this ion have been found 
approximately equal in a number of lakes, such as those in the Canadian Shield 
region (NRCC 1981). Long-term increases in sulfate deposition to sensitive 
lakes in Yellowstone could result in increased sulfur concentration and 
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potential acidification. As discussed previously, the mean sulfate concentra- 
tion for the low sulfate-low alkalinity group of lakes was approximately 36 
yeq/!l. This value corresponds wel! to the estimated average annual! precipita- 
tion input of about 33 peq/! as derived from recent NADP data and may, as an 
extreme, represent a 100% contribution from rainfall (no sulfate derived from 
soils, bedrock, or nearby thermal sources). This simplification represents a 
maximum, and ignores any potential soil sulfate adsorption and residence time 
considerations. 

Yellowstone does not currently receive acid precipitation comparable to 
the eastern United States. The low sulfate concentration of 36 peq/! falls in 
the range of “background” sulfate for the Canadian Shield lakes (30-60 peq/!), 
and probably is also almost entirely of natural origin. This is further 
corroborated by the relatively high rainfall pH. In contrast, lakes and 
Streams in acidified regions throughout the U.S. and Canada show average 
sulfate ranges of 3 to 4 times that (Table 21). The higher values of sulfate 
found in some Yellowstone lakes is assumed due to groundwater, soil, bedrock, 
or geothermal! sources. Yellowstone probably experiences little current acidi- 
fication from acid precipitation. Considering the highly variable lake 
sulfate concentrations, however, local sulfate sources undoubtedly pervade the 
park. Given the existence of geothermal, loca] atmospheric, bedrock, and 
underground spring sources of sulfate, it is impossible to state how anthro- 
pogenic acidification will affect these lakes. 


Table 21. Sulfate levels in selected acidified waters. 








Area Sulfate range (yeq/!) Reference 
Southern Ontario (lakes) 160-22u Dillon et al. (1978) 
Nova Scotia (lakes) 100-140 NRCC (1981) 

New Hampshire (streams) 129-142 Likens et al. (1977) 





It is possible, however, to approximate how Yellowstone would be affected 
if rainfall over the park were to change composition to approximate the 
acidity of typically eastern precipitation. The monitoring station at Hubbard 
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Brook experimental! forest provides the longest continuous record of precipita- 
tion chemistry in the United States. Over the ten year period from 1963 to 
1974, the weighted annual mean concentration of sulfate in precipitation was 
60.3 weq/1, and for nitrate, 23.7 yeq/l. A far shorter record (1.5 years) is 
available for the evaluation of precipitation chemistry of Yellowstone 
National Park, however, a value of approximately 30 yeq/! sulfate and 10 yeq/! 
nitrate is reasonable (NADP 1981). Thus, a doubling of the sulfate and 
nitrate concentration of Yellowstone precipitation from approximately 40 to 80 
yeq/!] would roughly yield the mean precipitation chemistry for these constit- 
uents in the east. Two extremes exist in the response of watersheds to this 
chemical change. As a minimum, all of this new sulfate and nitrate would be 
incorporated into the watershed, i.e., adsorbed onto soil particles or tied in 
biotic cycling. As a maximum, al] of the new strong acid-derived anions would 
be introduced into the lake ecosystem. Thus, the range of acidification of 
Yellowstone lakes would be from 0-40 weq/]. The maximum loss of alkalinity 
expected, under these conditions, is 40 weqg/1l; the true loss is probably 
somewhat less. A loss of 40 pec/!] alkalinity from each lake in Yellowstone 
would shift the number of “sensitive” (alkalinity <230 peg/1; no significant 
hot springs) lakes from 24 to 33. Only one of these new lakes has a signifi- 
cant trout population, the rest are historically barren. 
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FISH POPULATIONS 


INTRODUCTION 


The loss of fish populations in acidified lakes and streams has been the 
most publicized effect of acidic deposition. Although it is believed that 
other aquatic biota are also unfavorably affected by increasing acidity, 
effects on fish have been the most widely studied and documented. The basis 
for the RMAS was to determine if, in fact, there was any evidence that the 
acidification of waters in the most sensitive areas of the Rocky Mountain 
region had progressed to a point that might endanger any of the fish species. 
The evaluation relied on the acidity and alkalinity data from RMNP and YNP and 
published and unpublished fisheries information for the Rocky Mountain region. 
Assessment Approach , 

Assessment of the potential responses of fish populations in the Rocky 
Mountain region to acidification, resulting from increased atmospheric strong 
acid deposition, was approached in the following manner: 





l. Areas potentially sensitive to acidification were defined on the 
basis of relationships between water quality, geology, and soils, 
developed in representative watersheds in RMNP and YNP. 


2. Information on predominant fish species, fisheries, and management 
policies prevailing in potentially sensitive areas were obtained 
from the published literature and unpublished fishery survey data. 


3. The possible effects of acidification on the species at risk were 
evaluated from published information on life history patterns, in 
relation to expected changes in water quality. Potentially critical 
life history stages were identified on the basis of documented 
effects on species with comparable life histories and habitat re- 
quirements, in regions where acidification has already occurred. 


4. The implications of possible acidificat‘on effects on fish popula- 
tions for State and Federal management policies were considered, 
critical information gaps were identified, and future research 
direction recommended. 
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Evaluations in this study of water chemistry data obtained from RMNP and 
YNP and data from other studies in the Rocky Mountains have not revealed any 
instances of chrenic acidification at levels that would be detrimental to fish 
survival. Chronic acidification of surface water occurs where rates of base 
supply, derived from weathering processes in the drainage basins, are exceeded 
by rates of strong acid input from atmospheric depositions. Such conditions 
have been observed in sensitive areas of eastern North America and Scandinavia 
(Wright et al. 1980), where chronically acidified lakes and streams exhibit an 
absence of bicarbonate buffering, low pH (<5.0), and increased levels of 
cytentially toxic metals. Fish populations are generally absent or extremely 
stressed at these chronic levels of acidification (Schofield 1976; Muniz and 
Leivestad 1980). However, it is important to recognize that perturbations 
leading to the observed decimation of fish populations in these chronically 
acidified waters probably were initiated at much earlier, transitional stages 
of acidification, prior to the complete loss of bicarbonate buffering in the 
systems (Henriksen 1980; Schofield 1982). 

As indicated in the previous sections on water chemistry and geology, a 
large number of headwater drainages in the geologically sensitive areas of the 
Rocky Mountains exhibit extremely low rates of base supply, as indicated by 
the low alkalinity and cation levels. These drainage systems border on a 
transitional stage of acidification, at current levels of atmospheric strong 
acid deposition. Even minor excursions of increased strong acid loading, as 
might occur during the early stages of snow-melt, could episodically acidify 
these systems. The presently available water chemistry data represent late 
| snow-melt or baseflow conditions, hence the contemporary occurrence of such 

events is unknown. The probable consequences of such perturbations for 
indigenous fish populations, should they occur in these sensitive drainage 
systems, are considered in subsequent discussions of fish responses to water 
quality change. 


RESULTS AND DISCUSSION 


Fish Resources in Potentially Sensitive Areas of the Rocky Mountains 

The general areas in the Rocky Mountains potentially sensitive to acidi- 
fication are situated in the alpine (>3,200 m) and upper montane (2,500-3,200 
m) altitudinal zones. Sensitivity of specific drainage basins is tempered by 
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local variations in bedrock and surficial geology, as discussed previously. 
Fish species richness and standing crops are relatively low in these headwater 
drainages and many were originally devoid of native fish populations, because 
of natural barriers to colonization in the form of falls and stream gradients 
impassible to fish. However, most of the orginally barren headwater lakes 
were subsequently stocked with either native or exotic salmonids, beginning in 
the early 1880s (Pennak 1963). The cutthroat trout (Salmo clarki) is the on’y 
native trout in the area and it is represented by several subspecies re- 
stricted to specific drainage systems. The only other native species locally 
present in these sensitive watersheds are the grayling (Thymallus signifer), 
several species of the mountain sucker (Pantosteus spp.), sculpin (Cottus 
spp.), and the mountain whitefish (Prosopium williamsoni). Introduced, non- 
native species that have become regionally established include the brook trout 
(Salvelinus fontinalis), rainbow trout (Salmo gairdneri), and brown trout 
(Salmo trutta). The intentional and inadvertent introduction of these 
species, particularly the brook trout and rainbow trout, has been most detri- 
mental to the native cutthroat trout stocks, both because of competitive 
exclusion and hybridization. Additionally, indiscriminate hatchery plantings 
of cutthroat have so blurred the genetic integrity of this species, that only 
in a few high lakes and streams can the original subspecies still be distin- 
guished morphologically (Pennmak 1963). The problems associated with the 
preservation and maintenance of these subspecies wil! be discussed further in 
the section on wild trout management, with specific reference to the water- 
sheds in RMNP. 

The extensive fisheries data base for YNP provides both a detailed chron- 
ology of the changes in fish distribution and management policies for this 
region since the 1800s. The comprehensive summary of fish stocking activities 
in YNP by Varley (1981) provides the most recent, detailed information on the 
Status of fish populations in the Park. An earlier publication by Fromm 
(1941) provides a more anecdotal! account of the early hstory of fish surveys. 
Although eighteen species of fish (12 native species) are currently recorded 
as being present in YNP waters, their distribution is quite limited. Many of 
the lakes and streams in YNP were and stil! are © shless, primarily because of 
the physical limitations to colonization described earlier. Lakes and streams 
in approximately 40% of YNP (Figure 27) were estimated to be devoid of fish in 
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Figure 27. Cross-hatching shows the area in Yellowstone National Park 
found barren of fishes by Jordan in 1889, with the exception of sculpins 
in the Gibbon River above Gibbon Falls (from Fromm 1941). 


the 1800's (Jordan 1889). Early fishery management. efforts sought to estab- 
lish fishes in these barren waters through extensive plantings of both native 
and exotic species. Varley (1981) provides a detailed chronology of these 
stocking efforts. Many of the previously fishless headwater streams and lakes 
received introductory stockings during the period 1920-1935. Although fish 
populations were probably established temporarily in most of these waters by 
stocking the majority of the initially fishless lakes are presently barren 
once again. However, this is not the case for streams. Most of these 
historically fishless streams are currently supporting viable fish populations 
(Table 22). The difference probably reflects the lack of suitable spawning 
habitat in many of the smal] headwater lakes. In the previous section 23 
potentially sensitive lakes were identified in YNP (alkalinity <200 weq/1) 
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Table 22. Current and historical fish population status of lakes 
and streams in Yellowstone National Park (from Varley 1981). 














Category Number 
Historical Current Lakes Streams 
Fishless Fishless 29 2 
Fishless Fish 18 | 38 
Fish Fish 1] 17 
Unknown 13 16 
Total 70 73 





from the water chemistry data available for 107 lakes. Information on histor- 
ical and current fish status obtained for 1] of the sensitive lakes is sum- 
marized in Table 23. Only one of the lakes supported cutthroat trout histor- 
ically and the only other species currently present is the introduced brook 
trout in three of these lakes. 
Fisheries Management Policies in Potentially Sensitive Areas 

State and federal policies which must be considered in assessing poten- 
tial acidification impacts on fish populations include state level wild trout 
management policies, National Park Service aquatic resource management policy, 
and U.S. Fish and Wildlife Service implementation of the Endangered Species 
Act, including cooperative agency efforts. 

The NPS policies governing the management of aquatic ecosystems in desig- 








nated Natural Zones are particularly relevant to the sensitive areas identi- 
fiec in RMNP and YNP. The primary goal of resource management programs in 
natural zones is the preservation and restoration of native aquatic eco- 
Systems, including those waters originally barren of fish. Park waters 
falling into this category are allowed (or in some cases “rehabilitated") to 
revert to their original fishless condition. This would include the majority 
of sensitive waters in YNP and a large number in RMNP as well. Strict inter- 
pretation and adherence to this policy would seem to make the question of 
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Table 23. Total alkalinity and fish population status for sensitive 
lakes in Yellowstone Nationa! Park (from this study, and Varley 1981). 





Fish Status 








Lake name Alkalinity Historical Current 
(yeq/!) 
wrangler 40 Fichless Fishless 
Summit 60 Fishless Fishless 
Shelf 80 Fishless Fishless 
Mt. Everts 160 Fishless Fishless 
Ice 200 Fishless Fishless 
Ranger 160 Fishless a 
Obsidian 80 Fishless brook trout 
High 170 Fishless cutthroat 
Forest 192 Fishless cutthroat® 
Trilobite 200 Fishless brook trout 
Robinson 100 cutthroat brook trout® 





asratus questionable 


potential acidification one of largely academic concern, with little relevance 
for fisheries programs. However, there are two important exceptions to this 
policy which should prioritize concern for the potential acidification of 
these waters. In some previously fishless waters, populations of either 
native or non-native fishes have become wel! established. If these estab- 
lished populations represent distinctive sub-species (particularly of endan- 
gered or threatened indigenous species) or valuable genotypes of non-native 
species, the populations may be designated as “naturalized” and managed as 
integral components of the ecosystem. Similarly, should research indicate the 
need to stock or re-locate threatened or endangered species, fishiess or 
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formerly fishless waters may be utilized as refugia for these species (Watson 
1980; Jones 1980). 

The relatively recent evoiution and implementation of “Wild trout manage- 
ment" programs (most notably, the Colorado Division of Wildlife Commission's 
Wild Trout Policy, adopted in 1982) in the Rocky Mountain states is indicative 
of the widespread concern for the depletion of wild trout stocks. The objec- 
tives of these programs are diverse and although they are us@ oriented, they 
share a common objective with NPS and FWS programs for the protection of 
native species. Cooperative efforts are particularly evident and vita! in 
programs developed for the protection and restoration of endangered or 
threatened species (Johnson and Rinne 1982; Behnke and Zorn 1976). Genetic 
swamping, habitat destruction, and competition with introduced species are 
primary reasons for the .widespread demise of native trout in the Rocky 
Mountain Area (Behnke 1979). The Endangered Species Act provides a vehicle 
for the restoration of these western salmonids and recovery action programs 
are currently in place throughout the region. Implementation of these 
programs first requires listing of the species as endangered or threatened, 
which involves extensive study and documentation (Johnson and Rinne 1982). In 
addition to affording the listed species protection, under the conditions of 
the Act, habitats necessary for their survival are also afforded protection 
from adverse alteration or destruction. How this condition might apply to 
federal programs that influence the potential for acidification of these 
habitats is unclear at this time. The general components of recovery plans 
entail definition of the genetic purity of the population(s), survey and 
monitoring of existing populations, evaluation of habitat in existing and 
candidate reintroduction waters, and reestablishment of the native 
populations. 

Endangered or threatened Salmo species present in potentially acid sensi- 
tive regions of the Rocky Mountains include several subspecies of the cut- 
throat trout. These fish and their listings are given in Table 24. The life 
history patterns of these species ind the potential problems that acidifica- 
tion might impose on restoration programs in sensitive watersheds are consid- 
ered in the following section. 
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Table 24. 


Rocky Mountain Region. 


Endangered and threatened Salmo genotypes in the 








Common name Scientific name States Listing status 

Greenback cutthroat® Salmo clarki stomias Colorado Endangered 

Colorado cutthroat Salmo clarki pleuriticus Colorado, Threatened 
Wyoming 

Westslope cutthroat Salmo clarki subsp. Montana Endangered 

Eastslope cutthroat Salmo clarki subsp. Montana Endangered 

Snake River cutthroat Salmo clarki subsp. Wyoming Rare 

Rio Grande cutthroat Salmo clarki virginalis Colorado, Endangered 
New Mexico 

Gila trout® Salmo gilae New Mexico Endangered 

Arizona trout® Salmo apache Arizona Threatened 



































Indicates nationally recognized as endangered. 


Potential Impacts of Acidification on Fish Populations 

The few remaining native trout populations in interior regions of the 
Rocky Mountains persist only in small, isolated headwater drainages. These 
native cutthroat populations prefer and function best at lower temperatures 
than other species. With few exceptions, the cutthroat coexist and dominate 
introduced species only in these cold, headwater situations (Behnke 1979). 
Unfortunately, the displacement of these rare and endangered genotypes to 
headwater drainages also makes them most susceptible to potential acidifica- 
Because of the current endangered status of 





tion in these sensitive habitats. 
these fish, the Salmo clarki complex must be considered the primary species at 
risk and priority is given here to a consideration of potential acidification 





impacts on these populations. 

Given the restricted distributions of the subspecific cutthroat popula- 
tions and the geologic heterogeneity in watershed sensitivities to acidifica- 
tion in the Rocky Mountains, current juxtaposition of specific populations 
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Significantly determines their relative sensitivity to any regional increase 
in acidic deposition. The risk of possible extinction, resulting from acidi- 
fication, could certainly be lessened if management recovery teams had prior 
knowledge of the distribution of the populations relative to specific water- 
shed sensitivities to acidification. In the process of selecting new sites 
for reestablishment of cutthroat populations, watershed acid neutralizing 
capacity should also be considered among the selection criteria. The current 
distributions of the Colorado River cutthroat (S. c. pleuriticus) and the 
greenback cutthroat (S. c. stomias), relative to the marked differences in 
watershed sensitivity observed in RMNP exemplify this point. The Colorado 
River trout is native to the upper Colorado River basin (Figure 28) and 
efforts are being made to reestablish this fish in headwater sections located 
in RMNP. The greenback occupies the headwaters of the Arkansas and South 
Platte drainages (Figure 29), on the eastern side of the Continental Divide. 
Comparable efforts have been made to reintroduce this form into these areas of 
the ~-*k. The higher alkalinities in the upper Colorado drainage would cer- 
tainly favor the maintenance of S. c. pleuriticus in the advent of increased 
acid deposition, whereas the low alkalinities on the other side of the Divide 
would probably not inhibit acidification of the greenback's prime habitat. 
Behnke's and Zorn's (1976) prophetic suggestion that the greenback trout may 
be the most vulnerable of al] western trouts to extinction, would likely be 
realized with acidification. 

In addition to distribution, there are species specific life history 
characteristics that must be considered in assessing potential sensitivity to 
acidification. All of the interior western trouts of the genus Salmo have 
basically similar life histories. They spawn in the spring when water temper- 
atures reach 5.5-9.0°C, which can be anywhere from early April to June or 
July, depending on latitude and elevation. All are obligatory stream spawners 
and fry emergence occurs in early to mid-summer. Growth, maturation, and 
fecundity are variable, depending on prevailing temperature regimes and pro- 
ductivity of the local habitats (Behnke and Zorn 1976). The basic life his- 
tory pattern outlined above is markedly different from that exhibited by 
salmonid populations inhabiting waters of eastern North America and Scandi- 











navia, where acidification impacts have been described (Schofield 1976; Muniz 
and Leivestad 1980). The predominant salmonid species in these areas 
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Figure 28. Indigenous distribution of Salmo clarki pleuriticus 
(from Behnke and Zorn 1976). 





(Salvelinus fontinalis, Salvelinus namaycush, Salmo trutta, and Salmo salar) 
are all fall spawners and either stream or lake spawning may be locally preva- 








lent. Fry emergence for these species is in early spring, often coinciding 
with snowmelt periods when water quality is very poor in acidified areas. The 
fry (particularly during and shortly after hatching) of most of these species 
have also been found to be physiologically more sensitive to acidification 
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Figure 29. Indigenous distribution of Salmo clarki stomias 
(from Behnke and Zorn 1976). 


than either embryos or older fish (Spry et al. 1981; Baker and Schofield 
1980). For these reasons, the early life history stages of development have 
been identified as critical periods for survival in acidified habitats. Given 
the different life history pattern of the western Salmo sp. and uncertain 
water quality conditions that might prevail during the early life history of 
these fish, under an acidification regime, it is difficult to extrapolate 
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these findings from eastern North America and Scandinavia to the Rocky 
Mountain region. Water quality conditions would be more favorable for cut- 
throat fry if episodic acidification were to occur before hatching. 

A more likely critical period for the headwater populations of this 
species might be during late winter and early spring, prior to or during 
Spawning. Studies by Johnson and Webster (1977) demonstrated a marked avoid- 
ance of acidic water by spawning brook trout and Flick et al. (1982) noted a 
tendency for emigration of brook trout populations from lakes during episodes 
of acidic snowmelt. These latter observations are particularly relevant in 
terms of the known propensity for emigration by cutthroat populations inhabit- 
ing headwater streams subject to unfavorable winter temperature extremes 
(Bjornn 1971). According to Behnke and Zorn (1976), this was also the cause 
of a failed transplant of greenback trout in the North Big Thompson River in 
RMNP, where all the fish migrated downstream over a barrier during the winter 
months. A dense brook trout population below the barrier made it doubtful 
that the cutthroat population could sustain itself there. These behavioral 
responses to adverse environmental conditions suggest a subtle, but potential- 
ly devastating mechanism whereby even relatively minor acidification excur- 
sions in headwater trout refugia could lead to population extermination. 
Behavioral studies would be needed to define thresholds of acidification that 
induce avoidance in the form of downstream migration. It is quite likely that 
these thresholds (e.g., in terms of pH change) would be much lower than those 
determined by classical bioassay for definition of dose-response functions, 
where response is death or acute physiological stress. 

As noted above, the relevance of any discussion of species specific 
dose-response functions, as usually defined, is somewhat questionable at this 
point, given the uncertainties in determining critical life history stages and 
population/community level responses. However, there are some potentially 
important, physiological level questions that need to be considered in defin- 
ing the sensitivity of this species to acidification. Most comparisons of 
relative tolerance to acidity among salmonid species indicate that the rainbow 
trout (which is closely related to the cutthroat) is the most sensitive to low 
pH (Haines 1981). However, no studies have yet been conducted with the cut- 
throat to define its relative tolerance. Intraspecific variation in acid 
tolerance to acidity has been observed in brook trout (Flick et al. 1982) and 


93 

















brown trout (Gjedrem 1976), and the marked polytypic character of S. clarki 





(Trojnar and Behnke 1974) suggests that it might also exhibit significant ~ 
variation in acid tolerance among the defined subspecies. Determination of 
the extent of this variation in acid tolerance among the extant populations in 
the Rocky Mountain region would be potentially useful for rehabitation pro- 
grams, in the event of acidification. 

The potential for aluminum mobilization by acidification in Rocky 
Mountain watersheds is uncertain, but probably not as great as in eastern 
forested systems where soils exhibit marked accumulations of amorphous forms 
of aluminum that are readily mobilized by acidic deposition (Cronan and 
Schofield 1979). The enhanced toxicity of acidified waters containing alum- 
inum is well documented (Schofield and Trojnar 1980) and if aluminum mobili- 
zation is not an integral facet of the acidification process in Rocky Mountain 
Soils, then comparisons of biological responses between the two areas would be 
further complicated by significant differences in solution chemistry. 

Another important difference in solution chemistry between Rocky Mountain 
waters and those of eastern North America are the levels of dissolved calcium. 
The low calcium levels in the headwaters of the Rockies are more similar to 
the alpine waters of Scandinavia, than those of eastern North America. The 
significance of calcium as a mediator of gill membrane permeability and acid 
Stress was noted by Brown (1981). Below 1 mg/! of calcium (typical for many 
of the Rocky Mountain headwaters) trout would be susceptible to acid induced 
osmoregulatory stress at much higher pH levels than populations inhabiting 
higher (>2 mg/1, typical for eastern waters) calcium waters. However, poten- 
tial adaptation (acclimation) of resident trout populations to low calcium 
environments might ameliorate expected acid stress responses (Guthrie 1981). 
Again, it is diff*cult to generalize dose-response function when dealing with 
fish populations that have evolved under unique environmental conditions. The 
intraspecific genetic diversity of Salmo clarki is quite remarkable, but the 





potential adaptability of the genotypes of acidified environments remains to 
be determined. 
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CONCLUSIONS 


INTRODUCTION 


This section summarizes the conclusions reached concerning the four major 
objectives outlined in the introduction. The first two sections are specific 
to the two areas studied, RMNP and YNP. The third section evaluates the 
results of the water chemistry studies in terms of effect on several species 
of fish important to this region. Finally, the assessment of sensitivity of 
lakes and streams to acidification in the Rocky Mountain Region in general is 
based on studies in each of the two Parks. 


ROCKY MOUNTAIN NATIONAL PARK 


Soils and Geology 





1. The soils of the Rocky Mour.ain National Park at the higher eleva- 
tions (above 9,000 feet) are only slightly developed. They have 
surface organic horizons in some cases and usually have a darkened 
A-horizon underlain by slightly weathered materia] derived from 
granite, diorites and other rocks low in bases. The soils are 
coarse, low in clay, low in base cations and relatively acid (modal 
pH 5.2). The organic matter provides most of the cation exchange 
capacity and there is no evidence of sulfate adsorption capacity. 


2. The physical properties of these soils--coarse, high in sand, low in 
clay, and steeply sloping--promote rapid movement of water to the 
Streams and lakes. 


3. The low base saturation, low pH and low cation-exchange capacity 
provide little opportunity for neutralizing the acidic inputs or 
providing significant alkalinity to the water system. 


4. The soil itself is probably resistant to rapid acidification due to 
the abundance of relatively unweathered minerals. 


Surface Water Chemistry 





1. Primary mineral weathering appears to be the dominant mechanism 
determining the concentrations of base cations, silica and alka- 
linity throughout the Park. Factor analysis shows that the primary 
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mineral weathering of oligoclase and biotite accounts for almost 50% 
of the variance in stream chemistry in areas underlain by granite 
and biotite gneiss and schist. In areas that contain tertiary 
intrusive bedrock, the mineral weathering of mafic materials, sulfur 
bearing minerals, and oligioclase account for more than 50% of the 
variance in stream chemistry. 


Atmospheric deposition is the primary source of chloride and nitrate 
in the streams of the Park. Atmospheric deposition is also the 
primary source of sulfate in the Park waters, with the exception of 
the Upper Colorado and Upper Fal! River basins where the weathering 
of sulfur bearing minerals is a source of sulfate. Chloride and 
sulfate are relatively constant with elevation, while nitrate con- 
centrations are highest above the timberline, where biological 
activity is lowest. Atmospheric deposition is also a significant 
source of streamwater C,, a result of deposition of airborne dust 
and salts from the dry, Windy, regions upwind of the Park. 


Most of the waters in the study watershed have alkalinities of <100 
yeq/|. The alkalinities are lowest at higher elevations. The 
waters of Glacier Gorge and Ypsilon Creek are extremely sensitive. 
Those of Roaring River and Upper East Inlet watersheds are very 
sensitive with alkalinities of <100 peq/!. The waters of Lower East 
Inlet and the higher elevations of Upper Fall River watersheds are 
slightly higher with alkalinities <200 peq/!. The waters in the 
lower elevations of Upper Fall River watershed are non-sensitive 
with alkalinities >200 peq/!. 


Based on the current concentrations of sulfate, nitrate, base 
cations and alkalinity we estimate that the waters of Glacier Gorge, 
Loch Vale, East Inlet and Ypsilon Creek watersheds may have suffered 
a smal] loss of alkalinity (<10 peq/]). We estimate that the waters 
in the Upper Colorado, Upper Fall River and Roaring River watersheds 
have suffered no loss of alkalinity. 


If precipitation in Rocky Mountain National Park becomes as acidic 
as in the eastern United States, we estimate that in the worst case 
(ji.e., if the increase in acidic deposition does not cause an 
increase in base cation loss from the watershed), most of the lakes 
in Rocky Mountain National Park will become acidified below pH 4.7. 
At the other extreme, if for every peq/! increase in acid sulfate in 
the waters there is a 0.4 wpea/! increase of base cation in the 
waters, we estimate using the Henriksen nomograph that the majority 
of lakes will become transitional and only a few will reach acid 
Status. The low ion-exchange capacity of the soils in the Park and 
the resistance of the bedrock to chemical weathering indicates that 
the increase in base cation concentrations with increased acidic 
deposition in the Park will be low and that with increasing acidic 
deposition many of the lakes will shift from a bicarbonate towards 
an acid status. 
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YELLOWSTONE NATIONAL PARK 


ie 


The region-chemistry t-tests and alkalinity-elevation plot showed 
that, after other influences are accounted for, surficial geology 
and elevation are correlated to water chemistry. Rhyolite, the 
dominant bedrock formation in the park, appears to be the most 
sensitive (conferring the least buffering capacity), although the 
subalpine (mean elevation = 2651 m) lakes located in the andesitic- 
volcaniclastic rocks of the northwestern corner also may be poten- 
tially sensitive. 


Lakes found in other regions of the park--the andesitic Absaroka 
Mountains, the mixed metamorphic Gallatin and Beartooth ranges-- 
appear not immediately sensitive to acidification by acid precipita- 
tion (alkalinities generally above 500 weq/!). These findings are 
corroborated by a streamwater chemistry study in the Absaroka 
Mountains (Miller and Drever 1977b), where mean alkalinity for 14 
stream samples in the Shoshone River Basin was found to be 650 
yeg/). 


Finally, no clear trends in surface water chemistry appear from 
examination of historical data. It is suggested that if any chem- 
istry changes have occurred at all, these are limited to the north- 
west corner lakes. Since estimated “baseline” (precipitation- 
derived) sulfate levels in Yellowstone lakes are still fairly low in 
comparison to those levels in known acidified lakes throughout the 
continent, any acidification that has occurred is minimal. Most of 
the park lakes appear well protected from acidification in the 
future. For the present, the dilute lakes of the rhyolite bedrock 
and Fall River basalt region, as wel! as the Specimen Creek drainage 
basin lakes, are most vulnerable to changes in chemistry by acid 
rain. 


EFFECTS OF CURRENT AND FUTURE ACIDIFICATION LEVELS IN FISH POPULATIONS 


IN YELLOWSTONE AND ROCKY MOUNTAIN NATIONAL PARKS 


l. 


Currently, there is no evidence of chronic acidification in Rocky 
Mountain waters that would be detrimental to fish survival. 


Tne very low base cation concentrations observed in the headwater 
drainages of RMNP suggests extreme sensitivity to acidification. 
Fish populations present in these low calcium waters may be particu- 
larly susceptible to osmoregulatory stress from episodic acidifica- 
tion. 


The few remaining native trout (Salmo clarki) populations located in 
interior regions of the Rocky Mountains persist only in small, 
isolated headwater drainages. The cutthroat coexists with and 
dominates introduced species only in cold, headwater situations. 
Displacement of these rare and endangered genotypes to headwater 
drainages also makes them most susceptible to potential acidifica- 
tion in these sensitive habitats. 
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4. Although a large proportion of the headwater lakes in YNP and other 
areas are historically fishless, and many are still in that condi- 
tion, these waters represent potentially invaluable retugia for the 
reestablishment of endangered species. 


5. Comparisons of alkalinity and bedrock geology in the upper Colorado 
and eastern Divide drainages in RMNP, suggests that the greenback 
trout (Salmo clarki stomias), which occupies the latter basins, is 
at potentially greater risk than the Colorado River cutthroat (Salmo 
clarki pleuriticus). 








6. A comparison of the life history patterns of the western Salmo with 
eastern salmonids, tn relation to seasonal changes in acidity, 
indicates that different life history stages may be impacted by 
acidification. Hatching and fry development have been identified as 
critical periods for the fall spawning eastern salmoids. However, 
early-late summer emergence of fry in the western Salmo populations 
(spring spawners) indicates this may be a less critical life history 


stage. 


7. Avoidance response, in the form of downstream emigration, to epi- 
sodic acidification in headwaters occupied by S. clarki is suggested 
as a more subtle, but potentially devasting impact of acidification 
in the Rocky Mountain watersheds. 


8. The relative sensitivities of the cutthroat genotypes to acidifica- 
tion stress, capacities for adaption in low calcium water, and 
dose-response functions are currently unavailable. 


SENSITIVITY EVALUATION OF THE CENTRAL ROCKY MOUNTAIN REGION 


The examination of the geochemistry of Rocky Mountain National Park has 
shown that many areas in RMNP are sensilive to acidic deposition and that this 
sensitivity is primarily determined by bedrock geology. In addition, sensi- 
tivity varies inversely with elevation in watersheds with consistent geology. 
The results of the evaluation of the existing data in YNP, though not so 
definitive, corroborate the general findings in the RMNP studies. An evalua- 
tion of the sensitivity of the Central Rocky Mountains (Colorado and Wyoming), 
using geologic maps and elevation, can be provided based on this information. 

The analyses in RMNP show that watersheds underlain by granite and bio- 
tite gneiss and schist are equally sensitive to acidic deposition. The lakes 
and streams in these watersheds had alkalinities <200 peq/!, while the waters 
at higher elevations (>3300 m) were very sensitive (alkalinity < 100 yeq/!). 
The Upper Colorado River Basin and the Upper Fal! River Basin contain tertiary 
intrusive rocks in their drainage, resulting in low (alkalinity »200 peq/!) 
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sensitivity. In YNP, low-alkalinity or sensitive lakes were found in regions 
underlain by rhyolite flows or basalt and andesitic-basaltic flows (USGS 1977; 
Tweto 1979). 

The analysis of sensitive aquatic systems in the Central Rocky Mountain 
Region has been accomplished by extrapolating the results from RMNP and thus 
delineates areas underlain by granite biotite gneiss and schist and similar 
gneisses and schists (Figures 30 and 31). Areas underlain by these formations 
are classified as sensitive (alkalinity «200 weq/!), lakes and streams located 
at higher elevations (>3300 m) can be classified as very sensitive (alkalinity 
<100 peq/!). Although areas of YNP have a limited number of moderately sensi- 
tive lakes (>200 peq/!), areas underlain by tertiary intrusive rocks are 
generally classified as nonsensitive (alkalinity »>200 yeq/1). 
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Figure 30. Spatial distribution of sensitivity in the central Rocky 
Mountain Region: Colorado. 
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Figure 31. Spatial distribution of sensitivity in the central Rocky 
Mountain Region: Wyoming. 


Mountain ranges underlain by rock formations similar to those in RMNP 
include the Sawatch Range, the Front Range (Lovering and Goddard 1950), and 
the Sangre de Cristo Mountains (Litsey 1958) in Colorado (Figure 30), and the 
Medicine Bow Mountains (Hill et al. 1968) and Wind River Range (Bassett and 
Giletti 1963) in Wyoming (Figure 31). All of these mountains contain granite, 
gneisses, schists, and tertiary intrusive rocks. 

Although this method of evaluation is adequate for a regional-scale 
sensitivity essessment, sensitivity evaluations for individual lakes and 
streams must be carried out on a drainage by drainage basis. Smal! outcrops 
of highly weatherable rocks, too smal! to be seen on regional-scale maps, may 
have an overwhelming influence on the sensitivity of a watershed. This is the 
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case Im tive Upper Fall River Basin in RMNP The bedrock in this basin is 
primarily granite and biotite gneiss and schist. However, a smal] deposit of 
tertiary intrusive rock at the head of the watershed gives the fall River a 
relatively high alkalinity (>200 peq/1) Differences in hydrologic flow path 


and sot! development may also dominate sensitivity on a local scale 
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RECOMMENDATIONS FOR RESEARCH AND MANAGEMENT 


INTRODUCTION 


While these studies and most others would indicate that areas of the 
Rocky Mountain West were not experienceing significant impacts from acid 
deposition, they also demonstrate that many headwater lakes and streams are 
very sensitive (alkalinities < 200 peg/e). It is therefore believed prudent 
to undertake some long term programs to more fully assess the current status 
of both deposition and surface water chemistries and to develop some long term 
measurement programs. In addition strategies should be developed to protect 
indigenous fish populations in the advent that this area does exper ence 
increased sulfate and nitrate deposition. The following recommendations are 
formulated to address these issues. 


ROCKY MOUNTAIN NATIONAL PARK AND ROCKY MOUNTAIN REGION 


l. Long term watershed experiments should be established at several 
points in the Rocky Mountain Region. At minimum, routine sampling 
(on at least a bi-weekly basis and in the spring, every 2 to 3 days) 
should be performed. The samples should be analyzed for all] major 
inorganic constituents, DOC, aluminum, alkalinity and total acidity. 
Based on the research of this project, we recommend that the Glacier 
Gorge watershed in RMNP be considered for long-term monitoring. (As 
a result of this recommendation the National Park Service is now 
conducting a long-term study in the Loch Vale subbasin of the 
Glacier Gorge watershed. ) 


2. Additional surveys of water chemistry should be performed in 
mountainous areas not only in the Rocky Mountain region but also in 
other mountainous areas of the western United States. It is most 
probable that the acidity of precipitation will increase in future 
years and it will be to our benefit to obtain background data at 
this time. 


3. The results from the watershed studies should be made available to 
research groups having models on watershed response to acidic depo: 
Sition. This will ensure that if there are pecularities about 
western watersheds, that the modeis wil! be developed with those 
pecularities taken into account. 
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Since one of the differences between the Rocky Mountain region and 
other regions experiencing acid rain is the greater importance of 
nitric acid relative to sulfuric acid, watershed studies should be 
designed to take this into account. 


A greater effort must be made to determine the rate of dry deposi- 
tion of neutral salts (e.g., CaSO,). 


YELLOWSTONE NATIONAL PARK 


F 


Future studies in Yellowstone National Park should concentrate on 
those regions which this report has determined as potentially sensi- 
tive to acidification. In particular, reliable water chemistry from 
the Specimen Creek drainage basin is necessary to determine if these 
lakes are indeed as dilute as indicated. 


A more intensive sampling of the low alkalinity lakes (i.e., several 
alkalinity measurements over the course of a year) would indicate 
seasonal fluctuations in alkalinity and provide a more complete 
estimate of sensitivity. It is suggested that the analytical tech- 
nique used to determine alkalinity of these lakes be one designed 
specifically for low alkalinity water, i.e., Gran's plot or double 
endpoint potentiometric titration. 


Finally, a survey of the headwater streams of the Park, particularly 
those in sensitive regions or with important fisheries, is necessary 
for a complete understanding of the response of Yellowstone to 
potential acidification. 


FISHERY RESEARCH AND MANAGEMENT 


Ae 


The potential sensitivity to acidification, of watersheds currently 
occupied by endangered or threatened populations of S. clarki or of 
watersheds considered as potential candidate sites for reintroduc- 
tion, should be determined and given consideration in recovery plans 
for the species. 


Any further experiments designed to develop dose-response functions 
for interior western Salmo clarki populations, should consider the 
possibility of significant variations in tolerance of this polytypic 
species. Additionally, these responses to acidification need to be 
determined in very low calcium media, typical of the headwater 
habitats of this species. 





The behavioral responses and emigration tendencies of S. clarki 
populations exposed to episodes of acidification should be deter- 
mined experimental ly. 


The potential for aluminum mobilization in Rocky Mountain watersheds 
exposed to increased acid deposition should be determined. 
103 


CT ¢ FAVA Ley 
3eST COPY AVAILABLE 














5. Future studies of winter and spring lake and stream chemistry in 
headwater catchments should also evaluate movement of trout popula- 
tions, in response to chemical change. 
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APPENDIX A 
SOIL AND WATER CHEMISTRY SAMPLING SITE MAPS, 
ROCKY MOUNTAIN NATIONAL PARK 
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— A~1. Soil and chemistry sampling sites in Upper Colorado basin, 
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Figure A-2. Soil and chemistry sampling sites in East Inlet basin, RMNP. 
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Figure A-3. Soil and chemistry sampling sites in Glacier Gorge basin, RMNP. 
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Figure A-4. Soil and chemistry sampling sites in Fall River basin, RMNP. 
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APPENDIX B 
- SOIL LABORATORY PROCEDURES, PROFILE DESCRIPTIONS AND CHEMICAL ANALYSES 


SOIL LABORATORY ANALYSIS PROCEDURES 


Cation Exchange Capacity 





Extractable acidity. To 10 g of air dry soil add 100 ml of extraction 
solution (BaCl., triethanolamine). Stopper, shake well, and let stand over- 
night. Run blank. Filter, wash with extracting solution. Bring up to 250 
ml. Add indicator (bromcreso] green and methy! red) and titrate with 0.2 N 
HC] to pink orange end-point. Calculate. (Black 1965) 





Base extraction. Use macerated paper in 60 m] fritted glass funnels. 
Extract with a total of 100 ml 1 N ammonium acetate pH 7.0 (allow to sit 
overnight in 30 ml extracting solution). Run Ca, Mg, K, and Na on extract. 
Calculate C.E.C. and report in meq. Calculate percent base saturation. 
(Black 1965) 





Potassium: 1.5 g soil in 15 ml of 1 N ammonium acetate, pH 7.0. Shake 5 
min at 200 opm. Filter. Read on AA. Compare with standards. (Carson 
1975) 








Calcium and magnesium: To 1.5 g soil add 15 ml 1 N ammonium acetate pH 
7.0. Shake 5 min at 200 opm a! filter. Dilute with Lanthanum Chloride 
(1500 ppm). Read on AA Spectrophotometer. Compare with standards. 
(Walsh 1971). 


Sodium: To 5 g soil add 15 ml 1 N ammonium acetate pH 7.0. Shake 5 min 
at 200 opm. Do not filter. Read supernatant liquid on AA Spectropho- 
tometer. Compare with standards. (Walsh 1971). 


Organic Matter 





Oxidize with potassium dichromate and conc. sulfuric acid by standing for 
30 min; add water, phosphoric acid, sodium flouride, and diphenylamine indi- 
cator. Titrate with 0.5 N ferrous ammonium sulfate. (Black 1965). 
pH 


1:1 soil to water volume. Mix 5 seconds. Let stand 10 min. Stir, read 
pH with pH meter. (McLean 1975). 
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Particle Size 





Samples are sieved to remove particles larger than 2 mm. The remaining 
sample is dispersed with sodium hexametaphosphate, placed in settling con- 
tainers and subsampled at appropriate intervals with a pipette. The subsample 
is dried at 105 C and weighed. (Franzmeier et al. 1977). 


Loss on Ignition (LOI) 





Ovendried (105 L) samples are ignited in a muffle furnace at 550 C and 
the loss in weight is determined. In soils with low clay contents, the loss 
is a valid estimate of organic matter. (Black 1965). 


Sulfate Absorption 





Soi] samples (10 g) were equilibrated with 25 ml] of dilute K2SO, solution 
at a pH of 4.1. The solutions were initially at 1, 5 and 10 mg sulfur per 
liter. After filtration the solutions were analyzed for remaining sulfate 
utilizing an indirect Atomic Absorption method. (@ien 1979). 
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Table B-l. Particle size distribution of soil samples. 





Percent of less than 2 mm 








3 Depth »>2 ™ Texture! 
Location (cm) (%) Sand Silt Clay class 
1-1 0-11 55 68 27 5 SL 
TG-1 11-25 60 82 13 5 S 
1G-2 O-17¢ 0 35 55 10 SIL 
16-4 0-5 1 21 66 13 SIL 
1G-4 5-11 2 21 60 19 SIL 
16-4 11-24 8 35 46 19 L 
AC] 0-8 - 52 39 9 L 
GC-1 Talus 66 91 6 3 LS 
GC-2 24-50¢ 56 83 14 3 $ 
GC-3 0-6 3 23 70 7 SIL 
GC-3 9.5-30 19 7$ 15 10 SL 
tI-] Talus 56 65 27 8 Si 
EI-2 0-9 0 20 75 5 SIL 
EI-3 3-13 58 39 45 16 L 
EI-3 13-28¢ 76 18 70 12 SIL 
EI-5 0-16 8 43 35 22 L 
EI-5 16-31¢ 10 43 34 23 L 
EI-6 0-4 23 51 4) 8 L 
Bk-] 0-Si+ 53 50 35 15 L 
WG-2 0-5 ? 34 54 12 SIL 
WG-2 5-35¢ 70 64 23 13 SL 
MN- 1 Talus 43 70 23 7 SL 
IL-2 6-0 -- 43 54 3 SIL 
180-1 5-0 =< 23 73 4 SIL 
L80-1 0-19 - 7 19 a $ 
180-1 19-54 -- 93 4 3 LS 
L60-1 54-80 -- 84 12 4 $ 
180-1 80- 8S¢ : 80 16 4 S 
L80-6 0-18 -< 63 28 9 S 
L80-6 18-50 -- 49 46 5 SL 
L80-6 50-80¢ =< 85 ll 4 $ 
L801) 0- 35+ -- 85 ll 4 $ 
(80-13 Muck -- 7 56 ? SIL 
L80-13 Grave! - 77 16 7 S 
(81-13 0-1 55 37 8 St 
(81-13 1-18 : 61 32 7 
(81-13 18-46 62 
(81-13 46-60+ 73 
181-16 31-42¢ 32 
L81-18 0-24 68 
(81-18 24-% 78 
L81-18 36-59+ 68 
vCS-1] 4-11 62 
vCS-1 11-21 83 
vCS-1 21-36 
vCS-] 45-60+ 82 





"See maps Appendix A 
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Table B-2. Soil chemical analyses®. 





Site Depth Base Each 





10° (cm) 3 LOI + Oo pH Exch. Ca Exch. Mg Exch. « Exch. Na Exch. # cEC sat bases 
6-1 0-11 13.8 10.2 6.0 14.28 2.12 0. 366 0.040 11.42 28.23 69.5 16.5 
6-1 11-25 4.20 3.40 5.7 3.03 0.586 0.119 6.035 4.35 8.09 46.23 37 
G-2 30-0 , 47.2 55. 20 5.2 6.59 0.630 0. 305 0.131 41.10 48. 7% 18.7 1? 
6-2 0-17 14.2 9.03 4.3 6.00 0. 880 0.181 0. 136 41.27 48.47 14.9 12 
iG-3 Talus 2.16 1. 30 5.2 5 25 i.01 0.194 0.043 10. 28 16. 78 387 6.6 
16-4 3-0 32.6 25.0 46 4.20 0. 980 0.875 0.078 41.88 48.01 12.8 6.2 
1G-4 0-5 13.0 10.4 4.8 0.9% 0. 346 0.334 0.077 38.31 40. 00 42 1? 
1G-4 $-11 11.4 

1G-4 11-24 9.21 

16-4 24-%6 4.92 3.19 5.0 0.208 0.034 0.075 0.031 25.95 26.29 1.3 ‘3 
AC-2 

AC-3 5-0 60.9 65.5 5.5 9.15 3.29 3.52 0. 732 33. 86 50.55 33.0 18.2 
AC-3 0-9.5 8.57 5 45 5.1 1.84 0.425 0. 209 0. 342 17.44 20. 26 13.9 2.8 
AC-3 9.5-12.5 11.2 

AC-3 12. 5-50¢ 9.18 1.62 4.9 1.84 0.731 0. 060 0.279 19.93 22.84 12.7 2.9 
GC-1 Talus 2.01 1.62 5.2 1.10 0.20 0.095 0.030 2.44 3.87 37.0 1.4 
GC-2 4-0 32.5 46.1 51 31.27 6.71 3.58 0.147 39.47 81.18 $1.4 41.6 
GC-2 0-24, 5 85 4.33 5.2 2.93 0.538 0.169 0.035 5.75 9.42 39.6 37 
GC-2 24-50" 3.16 3.25 5.0 1. 70 0.217 0.072 0.052 4.48 6.52 31.3 2.0 
GC-3 5-0 ~ 75.8 62.7 4.3 0.963 0.535 0.672 0.107 45.29 47.57 4.8 2.3 
GC-3 0-6 28.0 28.5 4.2 0.315 0.438 0. 467 0.187 41.03 42.44 3.3 1.3 
GC-3 6-9.5 18.0 

GC-3 9. 5-30 3.52 2.67 4.6 0.101 0.063 0.058 0.039 12.13 12. 40 2.2 0.2 
GC-4 4-0 75.5 65.4 4.9 14.58 2.07 1.25 0.131 34.07 $2.10 4.6 18.0 
GC-4 0-4 9.32 6.14 4.9 4.12 0.927 0.633 0.134 21.80 27.61 21.0 5.9 
GC-4 4-10, 2.21 0.911 5.1 1.11 0. 286 0.101 0.035 6.17 7.70 19.9 1.5 
GC-4 10-46 2.7% 

E1-0 5-0 

1-0 0-3 

1-0 3-13 

£1-0 13-24, 

EI-0 24-35 

EI-1 Talus 6.99 6.44 5.1 2.63 0.556 0.311 0.039 12.13 15.67 22.6 5 
EI-2 23-0 $1.9 47.7 47 7.53 1.29 0.937 0. 369 42.06 62.19 19.4 10.1 
EI-2 0-9 19.7 17.9 4.8 3.15 0.315 0. 086 0.132 35. 04 38.72 95 3? 
E1-2 9-27" 11.2 7.74 4.8 2.27 0.224 0.040 0 080 27.% 30.07 8.7 2.6 


/2/ 








Table B-2. (continued) 








Site, Depth Base Each 
1.0. (om) % LO! . Oo pe Exch. Ca Exch. Mg Exch. « Exch. Na Exch. # CEC sat bases 
EI-3 12-7 75.0 55.3 5.9 4.44 5.40 1.11 0.126 26.37 79.45 66.8 $3.1 
EI-3 7-0 60.3 59.5 6.1 72.% 7.93 0.581 0.272 26.81 107.9% $.2 81.2 
EI-3 0-3 29.7 27.6 6.3 50.83 5.80 0.412 0.214 22.43 79.69 71.9 $7.3 
EI-3 3-13 5. 26 
EI-3 13-28 2.67 2.17 5.9 8.99 1. 30 0.129 0.083 4.73 15.23 68.9 10.5 
EI-4 
EI-5 19-0 68.8 60.9 5.0 24.53 4.09 1.65 0.246 %. 82 67.¥ 45.3 30.5 
EI-5 0-16, 6.70 5.3 5.5 10. 82 1.88 0. 726 0.11) 15.09 28.63 47.3 13.5 
EI-5 16-31 5. % 4.93 5.6 8. % 1.63 0.647 0.121 14.88 26.24 43.3 11.4 
EI-6 0-4 19.5 
EI-6 4-10, 4.53 
El-6 10-3) 5.02 
= EI-7 20-0 18.8 17.1 5.1 8.76 1.88 0. 383 0.4é 23.48 4.4 32.8 11.5 
nN El-? 0-22 ll.4 12.4 5.3 13. 00 3.4% 0.1% 0.452 18. 20 35.31 48.5 17.1 
EI-7 22-37 3.84 1.98 5.5 10.71 2.98 0.042 0. 288 6.23 20.25 69.2 14.0 
EI-? 37-39 3.82 
SG-1 Talus 6.63 3.14 5.2 7. 2.04 0.458 0.044 12.45 22.33 “42 9.7? 
4 SG-2 17-0 “4.3 31.5 $.2 10.49 1.59 0.259 0.101 27.40 39. 84 31.2 12.4 
as SG-2 0-17 3. 66 2.15 $.2 3.40 0.717 0.138 0.052 8.95 13.26 32.5 4.3 
a“ 56-2 17-30, 4.15 
62 $G-2 30-41 7.09 2.57 $.1 2.53 0. 556 0.124 0.070 15.42 18. 70 17.5 3.3 
=" SG-3 2-0 69.0 38.5 4.7 13.39 4.99 1.87 0.115 33.33 53.70 37.9 20.0 
>~ SG-3 0-5 P 4.48 2. 6 4.5 1.92 1.60 0.285 0. 066 18.29 22.16 17.5 4.0 
«- S63 5-24 5.75 2.60 4.9 0. 909 0. 286 0.223 0. 062 17.47 18.95 7.8 1.5 
c BK] 2-0 P 40.9 38.1 5.3 15.41 5.4] 1.82 0.121 29.55 $2.31 43.5 22.6 
a Bx) 0-$1 5.60 3.58 5. 3.64 0. 740 0.311 0.101 17.47 22.26 21.5 4.8 
_ Bx) 14-16 6.62 
rv 
BK-2 Talus 
4G-1 Talus 5.42 2.6 5.8 6.57 0.84) 0.376 0.048 9.6 17.50 “48 8.1 
WG-2 4-0 50.8 52.3 3.5 1.28 2.42 1.50 0. 185 47.20 58.59 19.4 11.4 
WG-2 0-5 P 19.9 18.5 4.3 1.02 0.425 0. 393 0.137 4.87 36. 85 5.4 2.0 
WG-2 5-35 5.50 2.41 5.1 0.505 0. 185 0.093 0. 062 19.5 20. 38 42 0.8 
Di 3-0 48.0 31.1 5.4 23.32 5.4% 1.17 0. 106 31.2 61.53 49.2 30.1 
D1 0-1) 11.2 1% 5.6 14.79 4.08 0.602 0.115 19.3 38. 89 50.4 19.5 
D1 li-2i, 4.52 
D1 21-48 4 68 1.26 5.7 7.47 2.27 0.31] 0.070 10. 89 21.01 46.2 10.1 
MN) Talus 4.8) 3.55 5.5 4.60 0. 700 0.174 0.030 10.99 16.49 33.4 5.3 
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Table 8-2. (concluded) 








Site Depth Base Each 











1.0 . (cm) % Ll ‘oOo pe Exch Ca Each Mg Exch. & tach. Wa tach. # crc sat bases 
2 0-4 4.64 2.51 a4 2 60 0.730 0. 308 0. 082 13.07  ™ 22.3 7 
MN- 2 4-40, 3.45 LM $2 2.2 0.639 0.31% 01M 12.13 ss 21.9 15 
~-o 40-50 27.19 
li-l Talus 6.72 5.0 5.0 324 0.690 0.285 0 048 15. 42 20. 28 24 48 
lke 6-0 5.8 27.7 5.0 25.2 5. 92 1.07 0.125 22.85 $$. 2) 58 6 i? 2 
Ii-2 0-15 10.9 721 5.0 6 + 2.02 0 2% 0.305 271.38 0.2) 29.2 . 
L80-) 5-0 4? 3 78.79 5.9 a 48 46) 0.64} 0. O84 20.13 49 9) t 79.8 
80-1 0-19 1% 0.87 5.3 1.40 0. 483 0 108 0.017 2.65 46 $3] 2.0 
L80- 1 19-54 0.92 0.31 $3 0.80 0.183 0.079 0.030 1.0? 2.1] $i 1.0 
180-1 54-80, 0.71 0.22 $3 0. 80 0.180 0.077 0.017 0.61) 1 & 63.2 1 
L80- 1 80-85 0. 0 $8 0 0 ’ 1.37 1.0 
180-6 0-18 1% 
80-6 18-50, 1“ 
(80-6 50- 80 0.83 
80-1) o- 35° 2.39 
L80-13 Muck” “9 
L80- 13 Grave! 5 3 
181-13 0-] 5 34 
181-13 l-18 2.85 
181-13 18-46 2.60 
(81-13 46-60" 2.02 
81-16 0-10 76.8 
L81-16 10-31, 63.1 
L8l-l6 3) 42 27.0 
181-18 0-246 | 
81-18 26-3, 1M 
L8i-18 59 1.15 
vC$~] 4-1) 3% 
vCS*1 11-21 1.28 
vCS*1 21-369 1% 
vCS+] 45-60 1.6 
yCS=e 0-4 3 
CS} %-60" 0.72 
from Mcfee (unpub! ished) 


c nee maps Appendi« A 
Film container samp!es 











APPENDIX C 
WATER CHEMISTRY DATA, ROCKY MOUNTAIN NATIONAL PAR) 


The raw water chemical data collected in Rocky Mountain National Park are 
presented with the elevation of each sample site. Missing data are denoted by 
an asterisk. Sample location maps are included in Appendix A. The sample 
identification code fo! lows: 


Roaring River 

Ypsilon Creek Fall River Basin 
Upper Fall River 

Andrews Creek 

Loch Vale Glacier Gorge 
Glacier Gorge 

Upper Colorado River Basin Mw East Iniet 


“S84SE58S 


The sample code indicates whether the sample was taken from a stream (5S), 
lake inlet (1), lake outlet (0), or lake surface (L). Replicate samples are 
denoted by an ‘'R' after the sample number. 


The lake inlet and stream samples were used to calculate the mean concen- 
trations of the major anions and cations for each subbasin. 





Table C-l. Lake and stream pH and alkalinity data 
for Rocky Mountain National Park. 














Samp le 

1D Elev F pH? Temp° Cond? L pH® ark! 

Roaring River 
GAOL! 3511. 30 6.23 2.0 12.7 6.45 26.0 
GA02S 3499.10 6.95 8.0 19.2 6.98 96.0 
GAO4SL 3511. 30 6.63 13.0 9.6 6.81 47.0 
GA050 3511. 30 6.50 11.0 10.2 6.66 50.0 
GA080 3511. 30 6.50 11.5 10.6 6.66 44.0 
GA10I 3364. 99 7.05 11.0 19.2 7.09 98.0 
GA12L 3352.80 6.03 14.0 15.8 6.70 81.0 
GA160 3352.80 6.19 14.0 14.3 6.39 83.0 
GA18S 3291. 84 6.08 14.0 10.0 6.98 83.0 
GA20S 3194. 30 6.65 13.0 17.1 6.96 78.0 
GA24S 3169. 92 6.47 13.0 16.0 6.78 75.0 
GA26S 3017.52 6. 30 13.5 15.1 6. 96 95.0 
GA28S 2901. 70 6.40 13.0 13.5 6.87 77.0 
GA30S 2804. 16 6.09 13.0 13.2 7.01 72.0 
GA32S 2621.28 6.14 8.0 15.4 6.97 92.0 
GA32SR 2621.28 6.14 8.0 15.4 6.87 97.0 
GA46S 2926.08 6.00 14.0 12.3 6.83 57.0 
Ypsilon Creek 

GB021 3413.76 5.90 5.0 -- 6.43 31.0 
GB060 $413. 76 5.63 8.0 4 6. 56 31.0 
GB08 1 3352. 80 5. 96 9.0 20.9 6.56 35.0 
GBO9I 3352. 80 6.20 6.0 =< 6.78 66.0 
GB101 3352. 80 6.05 13.0 19.9 6.65 42.0 
GB120 3352.80 5.94 11.0 19.9 6. 66 42.0 
GB14] 3279.65 6.27 13.0 21.9 6.68 60.0 
GB160 3279.65 6.50 11.0 21.4 6.75 54.0 
GB18S 3108. 96 =< =< 11.8 6.69 63.0 
GB19] 3462.53 6.41 11.0 16.5 6.09 16.0 
GB220 3462.53 6.41 8.0 19.1 6.05 19.0 
GB260 3462.53 6.58 12.0 20.4 6.53 26.0 
GB30L 3462.53 6.83 14.0 22.2 6.58 34.0 
GB320 3462.53 6.78 13.0 19.0 6.61 38.0 
GB34S 3352.80 6.12 8.0 22.7 6.56 33.0 
GB35S 3352. 80 6.17 8.0 20.4 6. 33 28.0 
GB38! 3218.69 6.45 8.0 15.9 6.65 33.0 
GB40L 3218.69 7.00 13.0 3 6.59 42.0 
GB40LR 3218.69 7.00 13.0 15.1 6.65 36.0 
GB420 3218.69 6.83 13.0 15.5 6. 56 37.0 
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Table C-l. (continued) 

















ID Elev? F pd? Temp Cond? L pH® ae 
Upper Fal! River 
GC02S 3413.76 0.00 7.0 42.2 7.32 377.0 
GCO2SR 3413.76 0.00 7.0 42.2 7.61 358.0 
GC03S 3413.76 0.00 6.0 26.5 7.12 172.0 
GC04S 3535.68 0.00 6.0 10.0 6. 88 69.0 
GC06S 3523.49 6.99 5.0 15.0 7.13 138.0 
GCO08s 3523.49 6.89 5.0 39.6 7.47 357.0 
GC10S 3511. 30 6.90 5.0 24.8 7.07 216.0 
GC11S 3401.57 0.00 10.0 9.6 6.65 51.0 
GC12S 3401.57 0.00 10.0 22.5 7.07 184.0 
GC135S 2926.08 0.00 8.0 23.6 7.38 170.0 
GC14S 2816. 35 0.00 9.0 14.6 7.05 119.0 
GC16S 2767.58 0.00 9.0 21.9 7.21 162.0 
GC18s 2657. 86 6.13 7.0 22.1 7.17 167.0 
GC20S 2657. 86 6.23 8.0 14.3 6.42 87.0 
GC24S 2657. 86 6.20 8.0 20.4 6.93 146.0 
Tyndal! Gorge 
MA021 3084.58 5.66 8.0 27.4 6. 36 29.0 
MAO2ZIR 3084. 58 5. 66 8.0 27.4 6. 36 24.0 
MA060 3084. 58 5.61 12.0 23.4 6. 56 31.0 
MA081 3023.62 5.81 11.0 27.0 6.46 39.0 
MA120 3023.62 5.77 15.0 26.6 6.77 38.0 
MA16S 2804. 16 - -- 12.0 6.60 62.0 
Loch Vale 
MCO021 3474.72 6. 66 1.5 53.3 6. 46 46.0 
MCO4L 3474.72 6. 66 9.0 30.9 6. 46 28.0 
MC060 3474.72 6.48 8.0 31.3 5.98 36.0 
MC12S 3267.46 6.55 9.0 34.) 6. 38 27.0 
MC14S 3169.92 5.92 6.0 29.5 6.43 32.0 
MC14SR 3169.92 5.92 6.9 29.5 6.47 29.0 
MC 181 3328.42 6.21 2.0 36.3 6.04 17.0 
MC20L 3328.42 0.00 -- -- 6.53 22.0 
MC220 3328.42 7.05 11.0 31.9 6.42 32.0 
MC24] 3316.22 6.24 12.0 35.7 6.44 24.0 
MC280 3316.22 6.04 13.0 34.5 -- 32.0 
MC 301 3108. 96 6.51 13.0 38.8 6. 46 40.0 
MC32L 3108. 96 6.81 19.0 37.1 6.74 43.0 
MC 340 3108. 96 5.93 15.0 33.4 6. 66 47.0 
MC 36S 2987.04 5.94 16.0 10.0 6.57 43.0 
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Table C-l. (continued) 

















Sample 

ID Elev® F pH? Temp* Cond? L pH® alk’ 

Glacier Creek 
MDO21 3535.68 6.10 8.0 28.0 5.66 10.0 
MD060 3535.68 6.15 12.0 18.5 5.89 11.0 
MD081 3230.88 6.22 11.0 26.7 -- -- 
MD101 3230.88 6.21 6.0 20.8 6.09 30.0 
MD 140 3230.88 6.34 11.0 24.4 6. 56 52.0 
MD24] 3474.72 6.17 3.0 33.7 6.52 33.0 
MD26L 3474.72 6.82 12.0 26.6 6. 66 38.0 
MD280 3468.62 6.49 12.0 28.2 6.66 35.0 
MD32L 3413.76 6.89 14.0 29.3 6.65 43.0 
MD 340 3413.76 6.90 13.0 23.9 6.52 64.0 
MD38S 3145.54 6.47 14.0 25.5 6.63 53.0 
MD40S 3035.81 6.44 14.0 26.5 6.52 54.0 
MD420 3035.81 6.51 16.0 25.7 6.79 47.0 
MD44. 3035.81 6.48 18.0 26.5 6.66 50.0 
M0460 3035.81 6.61 15.0 26.3 6.75 45.0 
MD485S 2974.85 5.97 16.0 9.9 6.73 65.0 
MD50S 2865. 12 5.88 12.0 10.1 6.73 65.0 
MD525S 2755. 39 0.00 0.0 10.1 0.00 53.0 
Upper Colorado River 
Voll 3511. 30 0.00 11.0 10.5 6.91 85.0 
v02S 3230.88 0.00 0.0 76.0 7.93 682.0 
v030 3499.10 0.00 0.0 30.2 7.38 280.0 
v04S 3255. 26 0.00 0.0 29.6 7.62 242.0 
V06S 3182.11 0.00 0.0 53.0 7.61 427.0 
voss 3169. 92 0.00 0.0 56.4 7.69 607.0 
V10S 3108. 96 0.00 0.0 46.6 7.58 375.0 
V10SR 3108. 96 0.00 0.0 0.0 7.61 366.0 
v12S 3115. 06 6.95 10.0 72.9 8.03 580.0 
V145$ 3121.15 6.71 9.0 64.8 7.56 374.0 
V16S 3121.15 6.70 10.0 52.0 7.44 361.0 
v181 3486.91 6. 82 9.0 18.7 7.29 146.0 
V20L 3486.91 9.00 11.0 41.6 7.72 342.0 
¥220 3486.91 6.85 110.0 41.6 7.65 341.0 
v245S 3230. 88 7.19 7.0 45.8 7.58 357.0 
v25S 3121.15 0.00 0.0 0.0 7.69 367.0 
V28S 3121.15 6.70 8.0 64.0 7.76 480.0 
V28SR 3121.15 0.00 0.0 0.0 7.76 480.0 
V30S 3133.34 0.00 0.0 0.0 7.33 202.0 
V32S 3133.34 0.00 0.0 26.1 6. 56 50.0 
V34S 3133.34 0.00 0.0 0.0 0.00 0.0 
V345R 3133.34 0.00 0.0 0.0 0.00 0.0 
V36S 3133.34 0.00 0.0 48.9 7.78 443.0 
V38S 3108. 96 6.40 9.0 22.9 7.20 189.0 
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Table C-l. (concluded) 





Samp le b 








ID Elev? F pH Temp* Cond@ L pH® al! 
East Inlet 

MwO021 3316. 22 5.78 7.0 13.0 0.00 80.0 
Mw060 3316.22 8.13 11.0 12.0 6.92 80.0 
Mw08 I 3169. 92 0.00 12.0 11.4 6.79 72.0 
Mw120 3169. °2 0.00 16.0 10.9 6.84 93.0 
Mwi5! 3145.54 0.00 8.0 8.8 7.00 94.0 
Mwl6] 3145.54 0.00 15.0 10.2 6.80 68.0 
Mw200 3145.54 0.00 16.0 9.8 7.10 76.0 
MW221 3108. 96 0.00 13.0 9.3 6.73 60.0 
MwW260 3108. 96 0.00 15.5 10.0 6.84 83.0 
Mw30S 3084. 58 0.00 0.0 11.6 6. 96 77.0 
Mw32I 3023.62 0.00 0.0 12.6 6.85 73.0 
Mw 360 3011.42 0.00 17.0 10.2 6.81 91.0 
MwW38S 2865.12 0.00 0.0 12.3 6.76 80.0 
Mwe0S 2889.50 0.00 0.0 15.4 6.45 121.0 
MW40SR 2889.50 0.00 0.0 15.4 6. 86 120.0 
Mw425S 2865.12 0.00 0.0 14.4 6.64 94.0 
MwWw44S 2791.97 0.00 0.0 13.4 7.10 96.0 





pElev = elevation in meters. 

cf pH = field pH in SU. 

qlemp = temperature in a 

econd = conductivity in ymohs/cm. 
¢l pH = lab pH in SU. 


Talk = alkalinity in peq/!. 
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Table C-2. Lake and stream chemistry data for Rocky Mountain National Park. 














Samp le P b 

1D Elev Na . My Ca NH, Cl NO , +O 50, PO, 510, 

Roaring River 
GAOLI 3511. 30 9.70 3.45 13.4] 56.64 0.39 9. % 19. 39 31.94 0.33 40. 90 
GAO2S 3499.10 37.67 3.48 30.20 86.83 0.00 5.50 13.60 44.3) 0.14 144.32 
GA0SL 3511. 30 10.79 2.66 11.44 59.43 0.00 4.37 3. 32 24.2 0.00 35. 70 
GA0SO 3511. 30 10.05 2.71 10.45 55.14 0.00 5.05 3.31 27.11 0.00 2.33 
GA080 3511. 30 11.27 2. 8 11.93 55.79 0.00 8.10 9.93 30. 92 0.00 29.00 
GALO! 3365. 00 41.06 4.30 29.54 93.76 0.00 6.35 11.97 38.83 0. 00 142.13 
GAL2ZL 3352. 80 29.75 4.42 23.21 82.19 0.00 6.29 6.05 34. 50 0.09 97.43 
GA160 3352.80 29.41 3.50 22.% 79.69 0.00 5.98 5.79 33.52 0.09 95.22 
GA18S 3291. 80 30.45 3.91 18.51 59.53 0.00 5.78 0.35 37.52 0.00 101.34 
GA20S 3194. 30 31.62 3.45 18.93 56.54 0.00 4.65 6.84 %.19 0.09 105. 80 
GA24S 3169. 90 30.62 3.22 26.83 73.30 0.00 4.09 6.47 31.92 0.14 112. 86 
GA26S 3017.50 31. 32 2.58 21.64 64.07 1.05 6.26 3.11 29.13 0.14 94. 16 
GA28S 2901. 70 29.97 4.07 21.48 65. 32 0.00 4.80 7.00 32.71 0.13 112. 86 
GAI0S 2804. 20 30.54 4.07 19.83 57.98 0.00 7.50 6.95 32.13 0.14 116. 5 
GA325 2621. #0 39.89 4.02 28.64 77.3 0.00 5.58 5.77 38.37 0.09 147.73 
GA32SR 2621. 30 40.50 4.35 29.29 78.14 0.00 5.98 5.79 37.81 0.13 145. 98 
GASES 2926. 10 23.97 3.94 19.01 50.85 0.00 4.60 8.76 32.71 0.17 97 % 
Ypsilon Creek 

GBO21 3413. 80 17.62 5.04 11.68 38.07 1.11 3.27 18. 06 27.94 0.27 66.79 
GB060 3413. 80 19.97 3.58 10.78 42.61 0.00 2.28 20.21 3%. 3) 0.27 75.42 
GB08] 3352. 80 19.44 2.74 11.11 45 4 0.00 2.62 16. 00 35.73 0.22 75.69 
GBO9! 3352. 80 33.84 4.17 13.08 43.4% 0.00 3.64 0.26 24. 86 0.00 119. 86 
GB10L 3352. 80 24.71 4.48 14.40 52.59 0.00 4.88 14.89 37.02 0.22 “4M 
GB120 3352. 80 24.23 3.94 12.51 47.95 0.00 2.48 14.45 3%. 10 0.14 93.50 
GB14I 3279.60 32.02 4. 86 16.95 53.14 0.00 2.00 12.64 40.% 0.20 116.63 
GB 160 3279.60 31.45 4.60 16.54 50.85 0.00 1.35 4.98 37.92 0.00 114.05 
GB18S 3109. 00 29. 88 4.35 16.62 53.49 0.00 2.71 9.40 38.81 0.11 113.78 
GB191 3462. 50 3.52 2.12 5.60 15.12 0.00 0.14 8.53 18.01 0.20 =e 
GB220 3462. 50 8.92 3.45 }.44 27.15 0.00 0.20 12.11 25.92 0.14 31.47 
GB260 3462. 50 8.00 2.92 9.63 33.43 1.22 4.15 11.11 21.49 0.00 26.71 
GB 301 3462. 50 18. 40 6.55 21.31 %.13 0.33 3.81 9.74 29.50 0.00 54.40 
GB320 3462. 50 17.70 5.47 18. 76 34. 88 0.00 4.80 2.44 30.15 0.25 46.19 
GB 345 3352. 80 16. 66 4.12 11.52 31.54 0.00 2.26 13.82 21.26 0.27 78.95 
GB35S 3352. 80 16.31 3.76 12.59 30.34 0.00 2.51 13.24 24.26 0.27 64 35 
GB 38] 3218.70 16.18 4.4 12.59 35.58 0.00 2.82 11.61 27.17 0.19 64.95 
GBS0L 3218. 70 14. 88 3.53 11.93 4. 48 0.00 3.95 7.10 20.24 0.00 50.35 
GB40LR 3218. 70 18. 88 4 68 13.50 37.03 0.00 5.92 9.18 23. 86 0.08 62. 30 
GB420 3218. 70 15. 66 3.89 13.25 3%. 93 0.00 3.58 8.85 23.86 0.13 62.07 
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Table C-2. (continued) 














1D Elev* nae Mg Ca WH, C1 NO ,*N0., $0, PO, i0, 
Upper Fal! River 
GC02s 3413 80 31.75 9.23 162.11 198.15 0.00 5 8] 3.06 87.62 0.00 137.62 
GCOZSR 3413.80 2.84 8. 80 161.12 198.00 0.00 6.52 3.18 87.89 0.00 129.34 
GCO3S 3413. 80 29.49 8.31 80.97 120. 56 0.00 6.43 8 66 %.21 0.11 144.12 
GCO4S 3535. 70 23.10 7. 80 13.25 44.66 0.00 4.20 2.27 -- 0% 146. 22 
GC06S 3523. 50 %. 93 6.29 37.52 73.95 0.00 5 84 3.71 19.43 0.09 167.00 
Gcoe8s 3523. 50 44.98 10.54 137.75 200. 20 0.00 14.27 3.00 71.40 0.20 162.97 
GC 10S 3511. 30 40. 06 7.95 78.75 123.80 0.00 8. 80 3.37 44.25 0.00 164.23 
GC1IS 3401.60 30.58 2.84 12.59 30.84 0.00 3.55 5.35 16.68 0. 30 138.95 
GC12$ 3401. 60 3 76 6.68 68. 30 107.53 1.1] 9.93 1.74 37.83 0.02 167.98 
GC135S 2926.10 42.85 11.33 67. 48 102. 30 0.55 11.31 6.27 44.83 0.00 122. 48 
GC14S 2816. 40 46.54 2.33 35.63 66. 82 0.00 5 S6 1.00 %. 98 0.09 216.43 
GC16S 2767.60 46.4) 10.28 56.04 98. 00 0.00 8.55 4.63 49.70 0.00 185 88 
GC 18s 2657. 90 45. 46 9.72 57.93 99. 70 0.00 R 2% 7.02 6) 82 0.00 180. 12 
GC2Z05 2657.90 40.45 3.09 26.74 70.% 0.00 6.68 9.29 20.11 0.09 189.23 
GC245 2657. 90 47.11 8.2] 53.16 %.61 0.00 8.2] 6.81 46.04 0.14 194.93 
Tyndal! Gorge 
MAC? 3084. 60 14. 40 3.79 9.79 37.67 0.00 -- 1.6 . 0.20 48.29 
MAO2ZI& 3084 60 14.35 3.8] $9.7] 38. 32 0.00 2.88 6.53 30.09 0.27 46. 28 
MAO 3084 60 16.0] 425 12.84 54.14 0.00 3.33 19.56 32.92 0.08 51.% 
MALE! 3023. 60 15. %6 3.% 12.67 55.79 0.00 3.41 12. % 35.19 0.08 62.95 
MA120 3023.60 16.44 3.48 13.4] 57.63 0.00 3.75 11.19 31.65 0.09 54.5% 
MALES 2804. 20 18.27 3.20 17.12 62.43 0.00 3.92 12.48 31.65 0.08 70.31 
Loch Vale 
C02! 3474.70 27.97 ).65 33.08 138.72 0.67 12.72 57.90 109.02 1.19 119.11 
MCOSL 3474.70 1.22 174 6.42 22.% 0.00 1.97 -- -- 0.20 33.78 
MC 060 3474.70 7.39 1.89 6.42 22.06 0.00 1.07 4.98 -- 0.20 30.10 
mC 125 3267.50 12.18 2.12 7.24 29.09 0.00 0.59 10. 42 16. 70 0.20 51.29 
mC 145 3169. 90 12.48 2.10 10 45 %.93 0.39 4.06 14.37 21.45 0.14 45.74 
mC 1458 3169. 90 12.48 2.02 9.71 39.72 0.00 1.30 14.74 26.71 0.19 57.19 
mC 18! 3328.40 10. 16 1.74 10. 20 %.08 0.33 3.33 17.81 18.51 0.25 26.99 
MC 201 3328. 40 12.09 2.10 9.13 39. 22 0.00 3.13 10.35 30.27 0.19 M26 
mC 220 3328. 40 9.87 1.79 11.19 37.43 1.55 448 11.53 23.19 0.22 28.85 
MC 24] 3316. 20 9.92 1.71 11. % 37.77 1.44 4.3) 12.93 24.90 0.20 30.84 
mC 290 3316. 20 10. 18 1.97 11.44 40.4) 0.61 4.23 11.52 25.% 0.22 30.8) 
mC 30] 3109. 00 11.70 2.02 12.18 42.42 0.67 40 11.27 24.88 0.20 38. % 
mC 321 3109. 00 19. % 3.17 13.91 4471 1.55 931 10.69 27.49 0.31 % 33 
mC 349 3109.00 14.70 2.35 14.32 46 66 1.16 5.02 10.45 27.69 0.14 5.64 
MC 36S 2987.00 13.14 2.05 14.32 44.9) 0.67 5 64 11.02 28.25 0.17 %. 04 
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Table C-2. (continued) 














Ca wm Ci WO ,°M0, Oy m0, 510, 
biacter Crees 

18.8) 1.33 2.09 13.72 2.31 0.42 

23.80 1.6] 3.50 11.53 14.49 0.09 15.23 

1.97 0.00 3.75 0.18 18. 8 0.00 %.04 
4 0.00 1.4] 14.4) 15.6 0M 40. % 
4) 4) 0.00 2.31 3% 20. 0.00 46 06 
47.0% 0.00 2.79 22.99 15.57 0.3] 11 
48.50 0.00 3. 92 13.8] 17.03 0.19 % 2/7 
45 %& 1.05 3.55 13.19 16.03 0.27 55.75 
51.60 0.00 14) 12.32 18.47 0.23 40 75 
4) 4) 0 00 3.50 11.98 16.0] 0.16 $7.71 
53.59 0.00 3.4) 8.13 13.9) 0.11 52. 66 
5344 0.00 3.4) 764 493 0.08 4. % 
50.45 0.00 6 % 7.63 13.33 0.09 53.25 
51.55 0.00 4.0) 6.42 1]. 66 0.00 $2.22 
$1.75 0.00 4.4) 6 Ww 24.53 0.00 42.62 
51.9 0.00 40) 8 8) 13.26 0.17 48.77 
4 O04 0.00 4.9%) 8.32 14.85 0.09 47.99 
57.14 0.00 1.69 17.72 10. 39 0.3) 60. 38 
Upper colorado River 

77.00 0.00 4 68 2.13 20.47 0.23 9.82 
1? 6 0.00 9 % 1.63 133. #0 0.1) 182. 16 
197.85 0.72 0.08 8.32 28.63 0.17 30. 88 
118. % 0.00 ) 42 7.03 48.22 0.27 480.9) 
257.14 0.00 7.33 3.60 101. 88 0.13 203.9% 
404.09 0.00 1044 4.8) 153.12 0.13 118.62 
215.12 0.00 10.24 1.71 86.60 0.08 276. 39 
215.82 0.00 10.52 1.82 86.0) 0.09 274.00 
370. 76 0.00 10.0) 144 154.24 0.13 124. 85 
381.84 0.00 7.93 7.9% 159.93 0.00 105. 68 
337.62 1.72 7.64 2.40 144.94 0.00 107. 90 
116.87 0.00 1» 9. 45 -< 0.27 17. 
311.38 1.16 44 2.73 55.11 0M 51.24 
277.25 0.00 6.15 9.06 58.15 0.33 59.73 
290.77 0.00 8.83 16. 92 78.79 0.17 65.89 
262.72 0.00 15.09 0.92 79.9) 0.14 14.69 
4% 63 1.33 27.13 11.M 113.89 0.00 59.13 
47) 0.00 7” 17.13 107. % 0.00 70.62 
19]. 47 0.00 5 34 3.29 93.55 0.00 107. 98 
119.4] 0.00 2.20 12.37 119.8) 0.00 174. % 
101.60 0.00 2.82 0.79 44 4) 0.00 33.78 














Table C-2. (concluded) 





Ss 
4 ‘ 5 
a ee a i Se, SO °? 
Last Inlet 
Mw? 316 20 6. 83 2.25 6.79 119. % 0. © 5 6] ». 61.19 0.39 72.4 
Me ve 3306 15. % 2.33 14 % 14.74 0.0 10 95 i 2 2.37 0.4 58.65 
Me & 169 & 20. 62 } 48 15. 4% % 86 0 Ow 5.2 9 03 Ww 3 0M 71.69 
Mw 120 1}69 9 0.84 1s 4. 1.6] 1.22 4 KF 2.92 22.82 0.4 $5.1 
Mw lS! 3}45. 5 1? 06 2.20 15.63 Ss 40 0.33 423 1.73 | 0.1) 172.0] 
Mls! 3345. 277.38 2.33 13.74 73. ww 0.33 40 1.15 » 4 0 08 95.19 
Me 200 3345. 50 s% 4 1.6] 2. % WwW 0.44 3 92 0 $$ 5.67 0. 80.42 
Ma? ?! 3,09 Ox Me 1.64 2.5] 65. 92 0 00 4.20 mm. % 0. 84 92 
Ma 260 3109 OX 22.™ 1” 13.99 68. % 02 499 1 28.73 0. 71.4] 
_ MOS 3084 60 23.4 1” 2.99 63.67 0.0 4 20 1.6] 79.09 0.09 69.79 
w Mio 3023. 60 26.77 2.2 13.9) 68 4) 0 00 4.0) 2.8 28 67 0.0 68 49 
™ Me de 01). 40 23.19 2% 4 ee 0. 0 00 4% 1.3 27.82 0.09 64 9) 
Mw ie 27865. 10 24.23 ?.33 15. 22 70 % 0. 3.67 2. 27.80 0 68. SO 
Md 2ae89. & 97 2.89 22.95 % 4) 0 3.53 , #9 ? 40 0. 100. $7 
Mod 8 2889. 50 24 1.22 23.04 $2.22 0 +5 3.89 5.67 0.00 10) 4 
Med oS 2865. 10 4? 02 2.9? 16.05 48 0. 00 319 0. 4 5.2 0 ®& 211.87 
er Mods 2792.00 28.19 2.48 17.03 'S. 80 0. 2.93 2 28. 82 0. 93.57 
—_ 
i eee ~ ” ee — ee ee eee OO 
— M, lew in meters 
m4 Al) ions 1A preg } 50, A im } 
4 missing values . 
— 
= 
fn 
— 
IS. 9 
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Table D-l. Recent lake chemistry data for Yellowstone Nationa! Park Lakes. 











vet 


Sum of lon 

Lake Alk year Elev ph Ca Mg Na . $0, Ci cations balance 
Turbdid 0 1977 2389 4.20 419.2 31.9 522.0 4.7 3185.5 138.2 1377.8 82.8 
Beaver | 40 1976 2249 5.10 698.6 279.7 3654.0 95.6 2290.2 2482.5 4727.9 -1.8 
Buffalo 40 1976 2M) 6 9 #9.9 8.2 26.1 5.9 41.6 28.2 80.1 -31.3 
wrangler 40 1977 23% /.20 9.9 8.2 30.5 7.9 20.8 25.4 86.5 0.3 
Summit 60 1977 2607 6. 30 39.9 16.5 30.5 9.2 20.8 %.7 %.1 ~20.0 
Dewdrop 80 1979 2481 6.50 79.8 123.4 60.9 12.8 124.9 %.7 276.9 13.6 
Dryad 80 1977 2529 6.70 59.9 16.5 39.2 13.3 20.8 39.5 128.9 6.5 
Mary 80 1977 2519 7.00 99.8 16.5 60.9 4.0 20.8 42.3 211.2 38.4 
Obs idian &0 1980 2365 6.50 119.8 32.9 4.8 15.3 114.5 42.3 202.8 15.5 
Shelf 80 1979 2733 6.50 73.8 $i.i 26.1 i0.2 4i.6 i4.i is/.2 is 
Mirror 100 1975 2705 6.80 39.9 16.5 $2.2 12.3 20.8 28.2 120.9 -20.8 
Robinson 100 1976 1981 5.90 79.8 41.1 17.4 5.6 20.8 28.2 143.9 -3.5 
Scaup 104 1980 2420 6.75 109.8 32.9 60.9 20.5 $2.1 42.3 224.1 12.2 
High 120 1979 2675 7.00 159.7 41.1 78.3 51.1 124.9 22.6 330.2 21.0 
wWyodaho 120 1976 2071 6.75 79.8 41.1 47.9 12.5 41.6 28.2 181.3 4.6 
Sedge 130 1970 2484 7.10 68.9 19.7 30.0 24.0 73.3 4.2 142.6 442 
Beach 146 1970 2484 6.90 
orag 160 1979 2683 6.80 79.6 41.1 30.5 15.3 41.6 19.8 166.7 28.2 
Crescent 160 1979 2622 7.20 119.8 57.6 39.2 23.0 145.7 31.0 239.6 -33.7 
mt. Everts A 160 1976 2211 6. 90 159.7 41.1 21.8 4.6 62.5 0.0 227.2 2.1 
Ranger 160 1976 2128 7.00 99.8 57.6 4.8 11.0 41.6 28.2 203.2 “12.3 
Sheridan 160 1978 2249 8.75 79.8 123.4 213.2 20.5 229.0 141.1 4%.9 -19.3 
Forest 192 1980 2262 7.00 174.7 a4) 60.9 10.2 43.7 42.3 270.5 “2.7 
Cygnet A 200 1978 2530 6.90 79.8 41.1 174.0 15.3 20.8 39.5 310.2 17.5 
Ice 200 1980 2402 7.10 14.7 32.9 78.3 0.7 50.0 4 276.6 0.4 
Mallard 200 1973 2447 6.80 
Trilobite 8 200 1978 2676 7.80 548.9 181.0 26.1 5.1 41.6 28.2 761.1 95.3 
Delacy E 220 1972 25% 6.50 
Delacy Ww 220 1972 25% 6.60 
Harlequin 220 1979 2100 5.30 379.2 123.4 365.4 89.5 374.8 2.5 957.5 4%.3 
Cascade 240 1978 2432 7.20 159.7 41.1 156.6 12.8 135.3 53 6 370.2 -4.7 
Cygnet 8 240 1978 2530 6.60 49.9 131.6 230.6 12.8 20.8 48.0 424.9 31.6 
Cygnet C 240 1978 2530 6.60 119.8 41.1 200.1 23.0 41.6 33.9 384.0 19.6 
Cygnet £ 240 1978 2530 6.75 119.8 41.1 178.4 23.0 41.6 39.5 %2.3 12.1 
Duck 240 1979 2374 6.80 254.5 82.3 91.4 33.2 20.8 28.2 461.4 45.9 
Lake of the Woods 240 1974 2348 6.00 159.7 0.0 74.0 40.9 187.4 0.0 274.6 43.5 
Virginia 240 1974 2499 6.80 119.8 0.0 95.7 23.0 41.6 0.0 238.5 ~ 16.6 
Wapiti 240 1975 2573 7.25 119.8 57.6 104.4 33.8 20.8 28.2 315.6 8.8 
Glade 246 1979 2952 7.00 119.8 57.6 78.3 17.9 208.2 28.2 273.6 -$5.2 
Unnamed 8 280 1976 1981 6.85 159.7 41.1 108.8 28.9 41.6 286.2 338.5 “3.3 
Lilypad 320 1976 1957 7.10 199.6 123.4 82.7 4.1 20.8 28.2 419.8 12.9 
Phone! ine 320 1976 1945 6.40 279.4 197.4 104.4 17 124.9 28.2 588.9 21.8 
Unnamed A 320 1976 1981 7.10 159.7 62.3 204.5 16.9 104.1 28.2 463.4 2.4 
Winegar 320 1976 1969 7.10 199.6 82.3 82.7 8.7 41.6 28.2 373.3 43 
Delusion 340 1971 2385 7.20 

2) 
134 











Table D-1. (continued) 











set 


Lake Ale rear t lew pe Ca Mg Na . 0, C catior balance 
fern 360 1975 2513 5.60 299.4 57.6 870.0 155.6 624.6 338.5 1382 4.4 
Pocket 360 1972 2470 6. 80 
wolf 360 1979 2438 7.40 239.5 $1.1 169.7 48.6 124 28.2 $98.9 2 
Cygnet 0 380 1979 2530 6.60 2593.5 139.8 204.5 6.8 232.0 13.9 60.¢ 4 
Hering 400 1973 2250 7.10 
Trilobite C 400 1978 2713 8. 0 279.4 4).1 8.7 5.1 41 6 3] 334.3 43 
Ribbon 440 1972 2384 7.00 
Grebe 480 1979 2448 7.20 239.5 156. 3 208.8 51.1 177.0 28.2 655 4.4 
Tanager 500 1973 2125 8.00 
YS w Thum 524 1980 2357 7 40 319.4 82.3 408.9 51.1 374.8 172.1 ih | 71.¢ 
Riddle 540 1979 2413 6. 80 239.5 238.6 130.5 0.0 20.8 28.8 608 . 6 3.3 
YS St Are 55? 1980 2357 7.40 279.4 123.4 356.7 4 0 208.2 143.9 605.5 i2.1 
YS S Are 557 1980 255? 7. 80 259.5 181.0 600.2 51.1 208.2 158.0 891.8 3.5 
Sho s hone 574 1980 2375 /. 249.5 98.7 826.5 61.4 4.3 297.6 1236.1 25.6 
YS Steve Island 574 1980 2357 7.40 279.4 197.4 391.5 48.6 229.0 172.1 916.9 6.2 
by ivan 606 197 2563 7.60 319.4 62.3 9.7 77 63.3 31.0 $05.2 4.5 
Lewis 639 1980 2372 7. 40 188.0 87.7 695.7 66.1 125.6 217.2 1057.5 4 
Panther Creek 640 1979 2591 6.50 479.0 156.3 30.5 ,.7 41.6 19.8 673.5 4) 
Heart 716 1980 2271 ?. 60 2794.4 151.0 1609.7 83.5 333.0 912.0 2138.6 8.7 
Alder 720 1976 2363 8.20 399.2 370.8 130.5 41.4 124.9 28.2 891.9 2.1 
Goose 760 1979 2198 6.00 239.5 2%. 1 6%. 0 }4.2 87.4 56.4 1305.8 %.4 
m. Everts 8 760 1976 2217 8.10 439.1 320.8 78.3 56.8 20.8 28.2 895. ( 10.1 
Beula 800 1973 2249 7.20 
Gooseneck 84U 1979 224) 7.20 469.1 41.1 609.0 61.4 104.1 28.2 1180.6 19.4 
Pass 840 1974 2793 9.00 $98.8 238.6 56.6 23.5 104.1 0.0 937.5 2.9 
Trilobite A 880 1978 2650 7.80 479.0 156.3 174.0 40.9 41.6 135.4 850.2 21.7 
Upper Goosenect 900 1979 2240 7.80 44.3 41.1 565.5 63.9 41.6 28.2 1014.8 45 
Lost 90 1980 2034 6.50 596.8 238.6 274.1 18.4 333.1 0.0 1129.9 13.5 
Sportsman 1040 1974 2357 8.50 648.3 238.6 14.9 47.1 41.6 0.0 1268.9 15.9 
Grizzly 1080 1979 2289 7. 80 698.6 361.9 $22.0 25.6 249.9 28.2 1608. | 16.9 
Rainbow ™ 1340 1974 1793 9.00 598.8 238.6 1131.0 33.8 333.) 0.0 2002.2 17.9 
Divide 1400 1970 2207 9.20 91.) 389.9 30.0 19.9 63.5 42 1369.9 6.9 
Gallatin 1440 1979 2693 8.40 548.9 904.9 30.5 17 104.1 28.2 1492.0 $.2 
Out let 1480 1978 2362 9.40 499.0 436.0 $13.3 7 62.5 42.3 1356.0 15.6 
Basin Creek 1580 1979 2253 8.10 1297.4 378.4 169.7 17.9 416 4 28.2 1863.4 “8.1 
McBride 1600 1975 i999 9.20 1247.5 238.6 104.4 44.5 20.8 28.2 1635.0 -0.9 
Squaw 1680 1979 2372 7,40 598.8 301.9 8917.5 227.6 541.3 7898.8 10106.0 0.1 
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Table D-1. (concluded) 
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Table 0-2. Chemistry data for Yellowstone Nationa! Park lakes 19/70 and earlier. 
lon 

Lake Year Alk Ca Mg Na % 50, Cl balance 
Cascade 1969 340 299.9 50.2 180.1 19.9 129.9 : 
Crag 1970 100 39.9 9.9 10.0 21.0 60.6 7.1 69.9 
Crescent 1970 270 150.2 29.6 20.0 27.1 64.5 11.2 -41.6 
Grebe 1963 440 159.7 79.8 159.6 68.3 50.0 14.1 -7.5 
High 1970 280 289.9 29.6 60.0 45.0 65.2 7.1 18.6 
Ice 1969 110 99.8 50.2 150.1 10.0 170.1 - 
Wolf 1969 270 250.0 50.2 180.1 29.9 170.1 31.0 7.9 
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